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If due attention is given in formulating the basic equations for the Gamma-Ray Burst (GRB) 
phenomenon and in performing the corresponding quantitative analysis, GRBs open a main avenue 
of inquiring on totally new physical and astrophysical regimes. This program is very likely one of the 
greatest computational efforts in physics and astrophysics and cannot be actuated using shortcuts. A 
systematic approach is needed which has been highlighted in t hree basic new paradigms: the relative 
space-time transformation (RSTT) paradigm (RufBni et al. Il44ll '). the interpretation of the burst 
structure (IBS ) par adigm (RufHni et al. 145]), the GRB-supernova time sequence (GSTS) paradigm 
(RufEni et al. |146| '). From the point of view of fundamental physics new regimes are explored: (1) 
the process of energy extraction from black holes; (2) the quantum and general relativistic effects 
of matter-antimatter creation near the black hole horizon; (3) the physics of ultrarelativisitc shock 
waves with Lorentz gamma factor 7 > 100. From the point of view of astronomy and astrophysics 
also new regimes are explored: (i) the occurrence of gravitational collapse to a black hole from a 
critical mass core of mass M > lOM©, which clearly differs from the values of the critical mass 
encountered in the study of stars "catalyzed at the endpoint of thermonuclear evolution" (white 
dwarfs and neutron stars); (ii) the extremely high efficiency of the spherical collapse to a black hole, 
where almost 99.99% of the core mass collapses leaving negligible remnant; (iii) the necessity of 
developing a fine tuning in the final phases of thermonuclear evolution of the stars, both for the star 
collapsing to the black hole and the surrounding ones, in order to explain the possible occurrence of 
the "induced gravitational collapse" . New regimes are as well encountered from the point of view of 
nature of GRBs: (I) the basic structure of GRBs is uniquely composed by a proper-GRB (P-GRB) 
and the afterglow; (II) the long bursts are then simply explained as the peak of the afterglow (the E- 
APE) and their observed time variability is explained in terms of inhomogeneities in the interstellar 
medium (ISM); (III) the short bursts are identified with the P-GRBs and the crucial information on 
general relativistic and vacuum polarization effects are encoded in their spectra and intensity time 
variability. A new class of space missions to acquire information on such extreme new regimes are 
urgently needed. 
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I. INTRODUCTION 



In understanding new astrophysical phenomena, the solution has been found as soon as the energy source of the 
phenomena has been identified. This has been the case for pulsars (see Hewish et al. [t^) where the rotational energy 
of the neutron star was identified as the energy source (see e.g. Gold |6l,|63). Similarly, in binary X-ray sources 
the accretion process from a normal companion star in the deep potential well of a neutron star or a black hole 
has clearly pointed to the gravitational energy of the accreting matter as the basic energy source and all the main 
features of the light curves of the sources have been clearly understood (Giacconi & Ruffini |6^). In this spirit, our 
work in the field of Gamma-Ray Bursts (GRBs) has focused to identify the energy extraction process from the black 
hole (Christodoulou & Ruffini 30]) as the basic energy sources for the GRB phenomenon: a distinguishing feature 
of this process is a theoretically predicted energetics of the source all the way up to 1.8 x 10^^ {Mbh /Mq) ergs for 
3.2Mq < AfsH < 7.2 X IO^Mq (Damour & Rufhni Ull). In particular, the very specific process of the formation of 
a "dyadosphere" , during the process of gravitational collapse leading to a black hole endowed with electromagnetic 
structure (EMBH), has been in dicat ed as ori gina ting and giving the initial boundary conditions of the onset of the 
GRB process (Preparata et al. |l24j . Ruffini 113711. Our model has been referred as "the EMBH model for GRBs", 
although the EMBH physics only determines the initial boundary conditions of the GRB process by specifying the 
physical parameters and spatial extension of the neutral electron_ positron plasma originating the phenomenon. 

Traditionally, following the observations of the Vela (Strong 175]) and CGRO^ satellites, GRBs have been charac- 
terized by few parameters such as the fluence, the characteristic duration (Tgg or T50) and the global time averaged 
spectral distribution (Band et al. 0). With the observations of BeppoSAX^ and the discovery of the afterglow, and 
the consequent optical identification, the distance of the GRB source has been determined and consequently the total 
energetics of the source has been added as a crucial parameter. 

The observed energetics of GRBs, coinciding for spherically symmetric explosions with the ones theoretically pre- 
dicted in (Damour & Ruffini 33]), has convinced us to develop in full details the EMBH model. For simplicity, 
we have consider ed th e vacuu m pol arization process occurring in an already formed Riessner-Nordstrom black hole 
(Preparata et al. |124| . Ruffini jl37| 1. whose dyadosphere has an energy E^ya- It is clear, however, that this is only 
an approximation to the real dynamical description of the process of gravitational collapse to an EMBH. In order 
to prepare the background for attacking this extremely complex dynamical process, we have clarified some basic 
theoretical issues, necessary to be implemented prior to the descripti on of the f ully dynamical process of gravitational 
collapse to an EMBH (Cherubini et al. |2S|, Ruffini & VitagUano [Tfi^ ITsTj . see section IXXVIl)l . We have then 
described the following five eras in our model. Era I: the e+e~ pairs plasma, initially at 7 = 1, expands away from 
the dya dosphere as a sharp pulse (the PEM pulse), reaching Lorentz gamma factor of the order of 100 (Ruffini et 
al. |l42j ). Era II: the PEM pulse, still optically thick, engulfs the remnant left over in the process of gravitational 
collapse of the progenitor star with a drastic reduction of the gamma factor; the m ass M b of this engulfed baryonic 
material is expressed by the dimensionless parameter B — Mbc^ / Edya (Ruffini et al. |l43l ]'). Era III. the newly formed 
pair-electromagnetic-baryonic (PEMB) pulse, composed of e+e" pair and of the electrons and baryons of the engulfed 



^ see http://cossc.gsfc.nasa.gov/i3atse7 
^ see ,http: / /www. asdc.asi.it/bepposax/ 



material, self-propels itself outward reaching in some sources Lorentz gamma factors of lO'^-lO^; this era stops when 
the transparency condition is reached and the emission of the proper-GRB (P-GRB) occurs (Bianco et al. jT^I)- Era 
IV: the resulting accelerated baryonic matter (ABM) pulse, ballistically expanding after the transparency condition 
has been reached, collides at ultrarelativistic velocities with the baryons and electrons of the interstellar matter (ISM) 
which is assumed to have a average constant number density, giving origin to the afterglow. Era V: this era represents 
the transition from the ultrarelativistic regime to the relativistic and then to the non relativistic ones (Ruffini et al. 

Ii42|). 

Our approach differs in many respect from the ones in the current literature. The major difference consists in the 
appropriate theoretical description of all the above five eras, as well as in the evaluation of the process of vacuum 
polarization originating the dyadosphere. The dynamical equations as well as the description of the phenomenon 
in the laboratory time and the time sequence carried by light signals recorded at the detector have been explicitly 
integrated (see e.g. Tab. and Ruffini et al. |l49l Il53| |'). In doing so we have also corrected a basic conceptual 
mistake, common to all the current works on GRBs, which led to the wrong spacetime parametrization of the GRB 
phenomenon, preempting all these theoretical works from their predictive power. The description of the inner engine 
originating the GRBs has never been addressed in the necessary details in the literature. In this sense neither the 
specific boundary conditions originating in the dyadosphere nor the needed solutions of the relativistic hydrodynamic 
and pair equations for the first three eras described above have been considered. Only the treatment of the afterglow 
has b een widely conside red in the literature by the so-called "fireball model" (see e.g. Meszaros & Rees IS^jIsSli Piran 
[llfif . Rees & Meszaros |l28j and references therein). 

However, also in the description of the afterglow, which is represented by the two conceptually and technically 
simplest eras in our model, there are major differences between the works in the literature and our approach: 

a) Processes of synchrotron radiation and inverse Compton as well as an adiabatic expansion in the source generating 
the afterglow are usually adopted in the current literature. On the contrary, in our approach a "fully radiative" 
condition is systematically adopted in the description of the X-ray and 7-ray emission of the afterglow. The basic 
microphysical emission process is traced back to the physics of shock waves as considered by Zel'dovich & Rayzer 
|l9.l| . A special attention is given to identify such p roce sses in the comoving frame of the shock front generating the 
observed spectra of the afterglow (see Ruffini et al. |l5Clj| '). 

b) In the literature the variation of the gamma Lorentz factor during the afterglow is expressed by a unique power-law 
of the radial co-ordinate of the source and a similar power-law relation is assumed also between the radial coordinate 
of the source and the asymptotic observer frame time. Such simple approximations appear to be quite inadequate and 
do contrast with the almost hundred pages summarizing the needed computations which we recall in the rest of this 
article. In our approach the dynamical equations of the source are integrated self-consistently with the constitutive 
equations relating the observer frame time to the laboratory time and the bound ary cond i tions are adopted and 
uniquely determined by each previous era of the GRB source (see e.g. Ruffini et al. EM EH IHS 1113 ) • 

c) At variance with the many power-laws for the observed afterglow flux found in the literature, our treatment 
naturally leads to a "golden value" for the power-law index n = —1.6. The fit of the EMBH model to the observed 
afterglow data fixes the only two free parameters of our theory: the E^ya and the B parameter, measu ring the remnant 
mass left over by the gravitational collapse of the progenitor star (Ruffini et al. |l48i Il49l Il5(l HH^ 'l. 

It is not surprising that such large differences in the theoretical treatment have led to a different interpretation of 
the GRB phenomenon as well as to the identification of new fundamental physical regimes. The introduction of new 
interpretative paradigms has been necessary and the theory has been confirmed by the observation to extremely high 
accuracy. 

In particular from the definition of the complete space-time coordinates of the GRB phenomenon as a function of 
the radial coordinate, the comoving time, the laboratory time, the arrival time and the arrival time at the detector, 
expressed in Tab. it has been concluded that in no way a description of a given era is possible in the GRB 
phenomena without the knowledge of the previous ones. Therefore the afterglow as such cannot be interpreted unless 
all the previous eras have been correctly computed and estimated. It has also become clear that a great accuracy in the 
analysis of each era is necessary in order to identify the theoretically predicted features with the observed ones. If this is 
done, the GRB phenomena presents an extraordinary and extremely precise correspondence between the theoretically 
predicted features and the observations leading to the exploration of totally new physical and astrophysical process 
with unprecedented accuracy. This has been expressed in the relative space-time transformation (RSTT) paradigm: 
"the necessary condition in order to interpret the GRB data, given in terms of the arrival time at the detector, is 
the knowledge of the entire worldlinc of the source from the gravitational collapse. In order to meet this condition, 
given a proper theoretical description and the correct constitutive equations, it is sufficient to know the energy of the 
dyadosphere and the mass of the remnant of the progenitor star" (Ruffini et al. (144] ). 

Having determined the two independent parameters of the EMBH model, namely Edya and B, by the fit of the 
after glow we have introduced a new interpretative paradigm for the burst structure: the IBS paradigm (Ruffini et al. 
|l45l| ). In it we reconsider the relative roles of the afterglow and the burst in the GRBs by defining in this complex 



phenomenon two new phases: 

1) the injector phase starting with the process of gravitational collapse, encompassing the above Eras I, II, III and 
ending with the emission of the Proper-GRB (P-GRB); 

2) the beam-target phase encompassing the above Eras IV and V giving rise to the afterglow. In particular in the 
afterglow three different regimes are present for the average bolometric intensity : one increasing with arrival time, 
a second one with an Extended Afterglow Peak Emission (E-APE) and finally one decreasing as a fu nctio n of the 
arrival time. Only this last one appears to have been considered in the current literature (Ruffini et al. |14It| |'1. 

The EMBH model allows, in the case of GRB 991216, to compute the intensity ratio of the afterglow to the P-GRB 
(1.45 • 10^^), and the arrival time of the P-GRB (8.413 • 10~^s) as well as the arrival time of the peak of the afterglow 
(19.87s) (see Figs. I12I6I11II . The fact that the theoretically predicted intensities coincide within a few percent with 
the observed ones and that the arrival time of the P-GRB and the peak of the afterglow also do coincide within a 
tenth of millisecond with the observed one can be certainly considered a clear success of the predictive power of the 
EMBH model. 

As a by-product of this successful analysis, we have reached the following conclusions: 

a) The most general GRB is composed by a P-GRB, an E-APE and the rest of the afterglow. The ratio between the 
P-GRB and the E-APE intensities is a function of the B parameter. 

b) In the limit B=0 all the energy is emitted in the P-GRB. These events represent the "short burst" class, for which 
no afterglows has been observed. 

c) The "long bursts" do not exist, they are just part of the afterglow, the E-APEs. 

We are currently verifying these theoretical predictions on the following GRBs: GRB 991216, GRB 980425, 
GRB 970228, GRB 980519. It is very remarkable that, although the energetics of GRB 980425 (see Fig. ^ dif- 
fers from the one of GRB 991216 by roughly five orders of magnitude, the model applies also to this case with success. 
Furthermore from these analysis we can claim that both in the case of GRB 991216 and in the case of GRB 980425 
there is not significant departure from spherical symmetry. 

While this analysis of the average bolometric intensity of GRB was going on in the radial approximation, we have 
proceeded to the full non-radial approximation, taking into account all the relativi stic corre ctions for the off-axis 
emission from the spherically symmetric expansion of the ABM pulse (Ruffini et al. |l48l Il5 We have so defined 
the temporal evolution of the ABM pulse visible area (see Fig. I13|l . as well as the equitemporal surfaces (see Fig. I13|) 
(Ruffini et al. jl4Slil53i| V 

We have then addressed the issue whether the fast temporal variations observed in the so-called long bursts, on time 
scales as short as fraction of a second (Ruffini et al. |l48j '). can indeed be explained as an effect of inhomogeneities in 
the interstellar medium. 

We are making further progress in identifying the basic mechanisms of energy release in the afterglow by presenting 
a new theoretical formalism which as a function of only one parameter fits the entire spectral distribution of the X-ray 
and 7-ray radiation in GRB 991216 (Ruffini et al. [l50l |'l. 

Finally the GRB-supernova time sequence (GS TS) paradigm introduces the concept of induced supernova explosion 
in the supernovae-GRB association (Ruffini et al. |146| |) leading to the very novel possibility of a process of gravitational 
collapse induced on a companion star in a very special evolution phase by the GRB explosion. 

Before concluding, we also present some theoretical developments which have been motivated by preparing the 
analysis of the general relativistic effects during the process of gravitational collapse itself and we also show how such 
results motivated by GRB studies have already generated new results in the fundamental understanding of black hole 
physics. 

In the next section we briefly summarize the main results and we will then give the summary of the treatment in 
the following sections. For the complete details we refer to the quoted papers. 

II. SUMMARY OF THE MAIN RESULTS 
A. The physical and astrophysical background 

Gamma-ray bursts (GRBs) are rapidly fueling one of the broadest scientific pursuit in the entire field of science, 
both in the observational and theoretical domains. Following the discovery of GRBs by the Vela satellites (Strong 
|175| ). the observations from the Compton satellite and BATSE had shown the isotropic distribution of the GRBs 
strongly suggesting a cosmological nature for their origin. It was still through the data of BATSE that the existence of 
two families of bursts, the "short bursts" and the "long bursts" was presented, opening an intense scientific dialogue 
on their origin still active today, see e.g. Schmidt 170] and section Ixlll 

An enormous momentum was gained in this field by the discovery of the afterglow phenomena by the BeppoSAX 
satellite and the optical identification of GRBs which have allowed the unequivocal identification of their sources at 
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Figure 1: Selected events in the EMBH theory are represented. For each one the values of the energy density of the medium 
and the distances from the EMBH, in the laboratory frame and in logarithmic scale, are given. 

cosmological distances (see e.g. Costa H^)- It has become apparent that fluxes of 10^'' erg/s are reached: during the 
peak emission the energy of a single GRB equals the energy emitted by all the stars of the Universe (see e.g. Ruffini 
158]). 

From an observational point of view, an unprecedented campaign of observations is at work using the largest 
deployment of observational techniques from space with the satellites CGRO-BATSE, Beppo-SAX, Chandra^, R- 
XTE'', XMM-Newton^, HETE-2^ as well as the HST^, and from the ground with optical (KECK^, VLT^) and radio 
(VLA^") observatories. The further possibility of examining correlations with the detection of ultra high energy cosmic 
rays, UHECR for short, and in coincidence neutrinos should be reachable in the near future thanks to developments 
of AUGER" and AMANDA^^ (see also Halzen |7|). 

From a theoretical point of view, GRBs offer comparable opportunities to develop entire new domains in yet 
untested directions of fundamental science. For the first time within the theory based on the vacuum polarization 
process occurring in an electromagnetic black hole, the EMBH theory, see Fig.^ the opportunity exists to theoretically 
approach the following fundamental issues: 

1. The extremely relativistic hydrodynamic phenomena of an electron-positron plasma expanding with sharply 
varying gamma factors in the range 10^ to 10** and the analysis of the very high energy collision of such an 
expanding plasma with baryonic matter reaching intensities 10'^* larger than the ones usually obtained in Earth- 
based accelerators. 

2. The bulk process of vacuum polarization created by overcritical electromagnetic fields, in the sense of Heisenberg, 
Euler (Heisenberg & Euler ,2SJ) Schwinger (Schwinger ^72]). This longly sought quantum ultrarelativistic 
effect has not been yet unequivocally observed in heavy ion collision on the Earth (see e.g. Ganz et al. Heinz 
et al. [i^Ij Leinberger et al. The difficulty of the heavy ion collision experiments appears to be that the 
overcritical field is reached only for time scales of the order H/mpC^, which is much shorter than the characteristic 
time for the e~^e~ pair creation process which is of the order of h/nieC^, where nip and nie are respectively the 
proton and the electron mass. It is therefore very possible that the first appearance of such an effect occurs in 
the present general relativistic context: in the strong electromagnetic fields developed in astrophysical conditions 
during the process of gravitational collapse to an EMBH, where no problem of confinement exists. 



^ see|http: //chandra. harvard. edu/| 

* see http: / /heasarc. gsfc.nasa.gov/do cs/xte/] 
^ see http://xmm.vilspa.esa.es/ 

^ see http: //space. mit.edu/ HETE/'l 
see http://www.stsci.edu/l 

* see http://www2.keck.hawaii.edu:3636/ 
® see http://www.eso.org/projects/vlt/ 

see http://www.aoc.nrao.edu/vla/html/VLAhome.shtml 
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see http://amanda.berkeley.edu/amanda/amanda.html 
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Figure 2; This qualitative diagram illustrates the relation between the laboratory time interval At and the arrival time interval 
Ata for a pulse moving with velocity v in the laboratory time (solid line). We have indicated here the case where the motion 
of the source has a nonzero acceleration. The arrival time is measured using light signals emitted by the pulse (dotted lines). 
Ro is the distance of the observer from the EMBH, to is the laboratory time corresponding to the o nset of the gravitational 
collapse, and r is the radius of the expanding pulse at a time t = to + At. See also Rufhni et al. (144| |. 

3. A novel form of energy source: the extractable energy of a black hole. The enormous energies released almost 
instantly in the observed GRBs, points to the possibility that for the first time we are witnessing the release of 
the extractable energy of an EMBH, during the process of gravitational collapse itself. This problem presents 
still some outstanding theoretical issues in black hole physics. Havi ng p rogressed in some of these issues (see 
Cherubini et al. |2^, Ruffini & Vitagliano jl5fl 115*^ . Rufhni et al. |159| |) we can now compute and have the 
opportunity to study all general relativistic as well as the associated ultrahigh energy quantum phenomena as 
the horizon of the EMBH is approached and is being formed (see section |XXVH|) . 

It is clear that in approaching such a vast new field of research, implying previously unobserved relativistic and 
quantum regimes, it is not possible to proceed as usual with an uncritical comparison of observational data to 
theoretical models within the classical schemes of astronomy and astrophysics. Some insight to the new approach 
needed can be gained from past experience in the interpretation of relativistic effects in high energy particle physics 
as well as from the explanation of some observed relativistic effects in the astrophysical domain. Those relativistic 
regimes, both in physics and astrophysics, are however much less extreme than those encountered now in GRBs. 

There are three major new features in relativistic systems which have to be properly taken into account: 

1. Practically all data on astronomical and astrophysical systems is acquired by using photon arrival times. It was 
Einstein |43| at the very initial steps of special relativity who cautioned about the use of such an arrival time 
analysis and stated that when dealing with objects in motion proper care should be taken in defining the time 
synchronization procedure in order to construct the correct space-time coordinate grid (see Fig. . It is not 
surprising that as soon as the first relativistic bulk motion effects were observed their interpretations within the 
classical framework of astrophysics led to the concept of "superluminal" motion. These were observations of 
extragalactic radio sources, with gamma facto rs ^ 10 (Biretta et al. |l6l |) and of mi croq uasars in our own galaxy 
with gamma factor ^ 5 (Mirabel & Rodriguez ?102]). It has been recognized (Rees |127| ) that no "superluminal" 
motion exists if the prescriptions indicated by Einstein are used in order to establish the correct space-time grid 
for the astrophysical systems. In the present context of GRBs, where the gamma factor can e asily surp ass 10^ 
and is very highly varying, this approximation breaks down (Bianco et al. T^, Ruffini et al. |144L Il53| ). The 
direct application of classical concepts in this context would lead to enormous "superluminal" behaviors (see e.g. 



Tab.QJ. An approach based on classical arrival time considerations as sometimes done in the current literature 
completely subverts the causal relation in the observed astrophysical phenomenon. 

2. One of the clear successes of relativistic field theories has been the understanding of the role of four-momentum 
conservation laws in multiparticle collisions and decays such as in the reaction: n —^ p + e~ + u^^. From the 
works of Pauli and Fermi it became clear how in such a process, contrary to the case of classical mechanics, it is 
impossible to analyze a single term of the decay, the electron or the proton or the neutrino or the neutron, out of 
the context of the global point of view of the relativistic conservation of the total four momentum of the system. 
This in turn involves the knowledge of the system during the entire decay process. These rules are routinely 
used by workers in high energy particle physics and have become part of their cultural background. If we apply 
these same rules to the case of the relativistic system of a GRB it is clear that it is just impossible to consider 
a part of the system, e.g. the afterglow, without taking into account the general conservation laws and whole 
relativistic history of the entire system. Especially since in astrophysics the "somewhat pathological" arrival 
time coordinate is basically used (see Fig.O. The description of the afterglow alone, as has been given at times 
in the literature, indeed possible within the framework of classical astronomy and astrophysics, is not viable 
in a relativistic astrophysics context where the space-time grid necessary for the description of the afterglow 
depends on the entire previous relativistic part of the worldline of the system (see also section IXVjl . 

3. The lifetime of a process has not an absolute meaning as special and general relativity have shown. It depends 
both on the inertial reference frame of the laboratory and of the observer and on their relative motion. Such 
a phenomenon, generally expressed in the "twin paradox" , has been extensively checked and confirmed to 
extre mely high accuracy as a byproduct of the elementary particle physics (g-2) experiment (see e.g. van Dick 
|l78j '). This situation is much more extreme in GRBs due to the very large (in the range 10^-10^) and time 
varying (on time scales ranging from fractions of seconds to months) gamma factors between the comoving frame 
and the far away observer (see Fig. IHJ. Moreover in the GRB context such an observer is also affected by the 
cosmological recession velocities of its local Lorentz frame. 

B. The Relative Space-Time Transformations: the RSTT paradigm and current scientific literature 

Here are so me of the reasons why we have presented a basic relative space-time transformation (RSTT) paradigm 
(Ruffini et al. |l44j ') to be applied prior to the interpretation of GRB data. 
The first step is the establishment of the governing equations relating: 

a) The comoving time of the pulse (r) 

b) The laboratory time {t) 

c) The arrival time at the detector (ta) 

d) The arrival time at the detector corrected for cosmological expansion (tf ) 

The book-keeping of the four different times and corresponding space variables must be done carefully in order to 
keep the correct causal relation in the time sequence of the events involved. 

As formulated the RSTT paradigm contains two parts: the first one is a necessary condition, the second one a 
sufficient condition. The first part reads: "the necessary condition in order to interpret the GRB data, given in 
terms of the arrival time at the detector, is the knowledge of the entire worldline of the source from the gravitational 
collapse" . 

Clearly such an approach is in contrast with articles in the current literature which emphasize either some too 
qualitative description of the sources and the quantitative description of the sole afterglow era. In this quantitative 
description they oversimplify the relations between the radial coordinate of the source and its gamma Lorentz factor 
as well as the relation between the radial coordinate and the arrival time using power-law relations which do not 
correctly take into account the complexity of the problem. 

In the current literature several attempts have addressed the issue of the sources of GRBs. They incl ude s cenarios 
of binary neutron stars mergers (see e.g. Eichler et al. Meszaros & Rees [13, 1131' Narayan et al. |104| ') blac k 
hole / white dwarf (Fryer et al. pg ) and black hole / neutron star binaries (Meszaros & Rees lOol. Paczyhski llOSl), 
hypernovae (see Paczyhski jllflj), fa i led s upernovae or collapsars (see MacFadyen & Wooslev |9f|. Woosley |l9fl| |). 
supranovae (see Vietri & Stella |l8lL Il83 |\ Only those based on Ijinary neutron stars have reached the stage of a 
definite model and detailed quantitative estimates have been made. In this case, however, various problems have 
surfaced: in the gener al en ergetics which cannot be greater than ^ 3 x 10^^ erg, in the explanation of "long bursts" 
(see Salmonson et al. '1621, Wilson et al. |187| |V and in the observed location of the GRB sources in star forming 
regions (see Bloom et al. ^J). In the remaining cases attention was directed to a qualitative analysis of the sources 
without addressing the overall problem from the source to the observations. Also generally missing are the necessary 



details to formulate the equations of the dynamical evolution of the system and to develop a complete theory to be 
compared with the observations. 

Other models in the literature have addressed the problem of only fitting the data of the afterglow observations by 
simple power-laws. They are separated into two major classes: 

The "internal shock model" , introduced by Rees & Meszaros Il28ll. by far the most popular one, has been developed 
in many different aspects, e.g. by Fenimore Fenimore et al. |5n||. Paczyhski & Xu ^109], Sari & Piran 165]. The 
underlying assumption is that all the variabilities of GRBs in the range At ~ 1 ms up to the overall duration T of 
the order of 50 s are determined by a yet undetermined "inner engine" . The difficulties of explaining the long time 
scale bursts by a single explosive mod el has evolved into a subclass of approaches assuming an "inner engine" with 
extended activity (see e.g. Piran |ll7j | and references therein). 

The "external shock model", see e.g. Cavallo & Rees fS^, Meszaros & Rees Shemi & Piran |17,'^ . is less 
popular today. Paradoxically, some of the authors who have qualitatively high light ed distinctive featur es of this 
model have later disclaimed its validity (see e.g. Meszaros & Rees fo^, Piran |116| . Rees & Meszaros '128] and 
references therein). Possibly they were carried to this extreme conclusion by an impressive sequence of mistakes they 
made in implementing the basic physical processes of the model. This model relates the GRB light curves and time 
variabilities to interactions of a single thin blast wave with clouds in the external medium. The interesting possibility 
has been also recognized within this model, that GRB light curves "are tomographic images of the density distribution 
of the medium surrounding the sources of GRBs" (Dermer & Mitman ^3|)i see also Dermer et al. ji^, Dermer |43| 
and references t here in. In this case, the structure of the burst is assumed not to depend directly on the "inner engine" 
(see e.g. Piran |117| and references therein). 

All these works encounter the above mentioned difficulty: they present cither a purely qualitative or phenomeno- 
logical or a piecewise description of the GRB phenomenon. By neglecting the earlier phases, the relation of the 
space-time grid to the photon arrival time is not properly estimated. To tell more explicitly, their clocks are out of 
the proper synchronization and the theory is emptied of any predictive power! 

We will explicitly show in the following how an unified description naturally leads to the identification of new 
characteristic features both in the burst and afterglow of GRBs. Our theory, in respect to the afterglow description, 
can be generally considered an "external shock model" and fits most satisfactorily all the observations. 

C. The EMBH Theory 

In a series of papers, we have developed the EMBH theory (Ruffini |137| ) which has the advantage, despite its 
simplicity, that all eras following the process of gravitational collapse are described by precise field equations which 
can then be numerically integrated. 

St arting from the vacuum polarization process a la Heisenberg-Euler-Schwinger (Heisenberg & Euler |78| , Schwinger 
|l72j| ) in the overcritical fi eld of an EMBH first computed in Damour & Ruffini 33] , we have developed the dyadosphere 
concept (Preparata et al. |l24j ). 

The dynamics of the e+e^-pairs and electromagnetic radiation of the plasma generated in the dyadosphere prop- 
agating away from the EMBH in a sharp pulse (PEM pulse) h as be en studied by the Rome group and validated by 
the numerical codes developed at Livermore Lab (Ruffini et al. |l42j ). 

The collision of the still optically thick e^e~-pairs and electromagnetic radiation plasma with the baryonic matter 
of the remnant of th e pro genitor star has been again studied by the Rome group and validated by the Livermore Lab 
codes (Ruffini et al. |l43j ). The further evolution of the sharp pulse of pairs, electromagnetic radiation and baryons 
(PEMB pulse) has been followed for increasing values of the gamma factor until the condition of transparency is 
reached (Bianco et al. 

As this PEMB pulse reaches transparency the proper GRB (P-GRB) is emitted (Ruffini et al. il45l] ) and a pulse 
of accelerated baryonic matter (the ABM pulse) is injected into the interstellar medium (ISM) giving rise to the 
afterglow. 

D. The GRB 991216 as a prototypical source 

In the early phases of development of our m odel, the EMBH theory was developed from first principles by the 
EMBH uniqueness theorem (Ruffini & Wheeler |l4lj |). the energetics of black hole (Christodoulou & Ruffini [20|) as 
well as the quantum description of the vacuum polarization process in overcritical electromagnetic fields (Damour & 
Ruffini j3^). Turning now to the afterglow, the variety of physical situations that can possibly be encountered are 
very large and far from unique: the description from first principles is just impossible. We have therefore proceeded 




to properly identify what we consider a prototypical GRB source and to develop a theoretical framework in close 
correspondence with the observational data. 

The criteria which have guided us in the selection of the GRB source to be used as a prototype before proceeding 
to an uncritical comparison with the theory are expressed in the following. It is now clear, since the observations of 
GRB 980425, GRB 991216, GRB 970514 and GRB 980326 that the afterglow phenomena can present, especially in 
the optical and radio wavelengths, features originating from phenomena spatially and causally distinct from the GRB 
phenomena. There is also the distinct possibility that phenomena related to a supernova can be erroneously attributed 
to a GRB. This problem has been clearly addressed by the GRB supernova time sequence (GSTS) paradigm in which 
the time sequence of the events in the GRB supernova pheno mena has been outlined (RufBni et al. |l46j| '). This has led 
to the novel concept of a n induced supernova (Rufhni et al. |146| |'). This problem will be addressed in a forthcoming 
paper (Ruffini et al. |152| V 

In view of these considerations we have selected GRB 991216 as a prototypical case (see Fig. OJ for the following 
reasons: 

1. GRB 991216 is one of the strongest GRBs in X-rays and is also quite general in the sense that it shows relevant 
cosmological effects. It radiates mainly in X-rays and in 7-rays and less than 3% is emitted in the optical and 
radio bands (see Halpern et al. |Zj|). 

2. The excellent data obtained by BATSE on the burst (BATSE Rap id Burst Response @) is complemented by 
the data on the afterglow acquired by Chandra (Piro et al., |l20l |') and RXTE (Co rbet & Smith |^). Also 
superb data have been obtained from spectroscopy of the iron lines (Piro et al., Il20l |'). 

3. A value for the slope of the ener gy emission during the afterglow as a function of time has been obtained: 
n = -1.64 (Takeshima et al. ^2^) and n = -1.616 ± 0.067 (Halpern et al. 17^1). 



E. The interpretation of the burst structure: the IBS paradigm and the different eras of the EMBH theory 

The comparison of the EMBH theory with the data of the GRB 991216 and its afterglow has natur ally l ed to a 
new paradigm for the interpretation of the burst structures (IBS paradigm)) of GRBs (Ruffini et al. |l45l |). The 
IBS paradigm reads: "In GRBs we can distinguish an injector phase and a beam-target phase. The injector phase 
includes the process of gravitational collapse, the formation of the dy ado sphere, as well as Era I (the PEM pulse). 
Era II (the engulfment of the baryonic matter of the remnant) and Era III (the PEMB pulse). The injector phase 
terminates with the P-GRB emission. The beam-target phase addresses the interaction of the ABM pulse, namely the 
beam generated during the injection phase, with the ISM as the target. It gives rise to the E-APE and the decaying 
part of the afterglow" . The detailed presentations of these results are a major topic in this article. 
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Figure 4: The dyadosphere of a Reissner-Nordstrom black hole can be represented as constituted by a concentric set of shells 
of capacitors, each one of thickness h/rrieC and producing a number of e'*'e~ pairs of the order of ~ Q/e on a time scale of 
10~^^ s, where Q is the EMBH charge. The shells extend in a region Ar, from the horizon r+ to the dyadosphere outer radius 
Vds (see text). The system evolves to a thermalised plasma configuration. 

We recall that the injector phase starts from the moment of gravitational collapse and encompasses the following 
eras: 

The zeroth Era: the formation of the dyadosphere. In section IIIII we review the basic scientific results which lie at 
the basis of the EMBH theory: the black hole uniqueness theorem, the mass formula of an EMBH, the process of 
vacuum polarization in the field of an EMBH. We also point out how after the discovery of the GRB afterglow the 
reexamination of these results has led to the novel concept of the dyadosphere of an EMBH. We have investigated 
this concept in the simplest possible case of an EMBH depending only on two parameters: the mass and charge, 
corresponding to the Reissner-Nordstrom spacetime. We recall the definition of the energy E^ya of the dyadosphere 
as well as the spatial distribution and energetics of the e+e" pairs. See Fig. ^ We return in section IXXVIII to the 
theoretical development of the time varying process lasting less than a second in the process of a realistic gravitational 
collapse. In reality the vacuum polarization process will lead to a final uncharged black hole, but the analysis based 
on a Re issner-Nordstrom black hole is an excellent approximation to the description of this phenomenon (Ruffini et 

ai. jie^). 

In order to analyse the time evolution of the dyadosphere we give in the three following sections the theoretical 
background for the needed equations. 

In section Hvl we give the general relativistic equations governing the hydrodynamics and the rate equations for the 
plasma of e"*^e~-pairs. 

In section we give the governing equations relating the comoving time r to the laboratory time t corresponding 
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Figure 5: Comparison of gamma factor for the one-dimensional (1-D) hydrodynamic calculations (Livermore code) and slab 
calculations (Rome code) as a function of the radial coordinate (in units of dyadosphere radius) in the laboratory frame. The 
calculations show an excellent agreement. 



to an inertial reference frame in which the EMBH is at rest and finally to the time measured at the detector ta which, 
to finally get if, must be corrected to take into account the cosmological expansion. 

In section fVll we describe the numerical integration of the hydrodynamical equations and the rate equation developed 
by the Rome and Livermore groups. This entire research program could never have materialized without the fortunate 
interaction between the complementary computational techniques developed by these two groups. The validation of 
the results of the Rome group by the fully general relativistic Livermore codes has been essential both from the point 
of view of the validity of the numerical results and the interpretation of the scientific content of the results. 

The Era I: the PEM pulse. In section IIVI by the direct comparison of the integrations performed with the Rome 
and Livermore codes we show that among all possible geometries the e"'"e~ plasma moves outward from the EMBH 
reaching a very unique relativistic configuration: the plasma self-organizes in a sharp pulse which expands in the 
comoving frame exactly by the amount which compensates for the Lorentz contraction in the laboratory frame. The 
sharp pulse remains of constant thickness in the laboratory frame and self-propels outwards reaching ultrarelativistic 
regimes, with gamma factors larger than 10^, in a few dyadosphere crossing times. We recall that, in analogy with the 
electromagnetic (EM) pulse observed in a thermonuclear explosion on the Earth, we have defined this more energetic 
pulse formed of electron-positron pairs and electromagnetic radiation a pair-electromagnetic-pulse or PEM pulse. 

The Era IT. We describe the interaction of the PEM pulse with the baryonic remnant of mass AIb left over from 
the gravitational collapse of the progenitor star. We give the details of the decrease of the gamma factor and the 
corresponding increase in the internal energy during the collision. The dimensionless parameter B = Mbc^ / Edya 
which measures the baryonic mass of the remnant in units of the Edya is introduced. This is the second fundamental 
free parameter of the EMBH theory. 

The Era III: We describe in section Hxl the further expansion of the e+e~ plasma, after the engulfmcnt of the 
baryonic remnant of the progenitor star. By direct comparison of the results of integration obtained with the Rome 
and the Livermore codes it is shown how the pair-electromagnetic-baryon (PEMB) plasma further expands and self 
organizes in a sharp pulse of constant length in the laboratory frame (see Fig.EJ. We have examined the formation 
of this PEMB pulse in a wide range of values 10~* < B < 10~^ of the parameter i?, the upper limit corresponding to 
the limit of validity of the theoretical framework developed. 

In section Ixl it is shown how the effect of baryonic matter of the remnant, expressed by the parameter i?, is to 
smear out all the detailed information on the EMBH parameters. The evolution of the PEMB pulse is shown to 
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Figure 6: Left) At the transparent point, the energy radiated in the P-GRB (the sohd hne) and (the dashed hne) the final 
kinetic energy of baryonic matter, EBaryons, in units of the total energy of the dyadosphere (Edya), are plotted as functions of 
the B parameter. Right) The energy corresponding to the peak of the photon number spectrum in the P-GRB as measured 
in the laboratory frame is plotted as function of the B parameter. 

depend only on Edya and B: the PEMB pulse is degenerate in the mass and charge parameters of the EMBH and 
rather independent of the exact location of the baryonic matter of the remnant. 

In section IxH the relevant thermodynamical quantities of the PEMB pulse, the temperature in the different frames 
and the e'^e~ pair densities, are given and the approach to the transparency condition is examined. Particular 
attention is given to the gradual transfer of the energy of the dyadosphere Edya to the kinetic energy of the baryons 
Esaryons during the optically thick part of the PEMB pulse. 

In section FXIII as the condition of transparency is reached, the injector phase is concluded with the emission of a 
sharp burst of electromagnetic radiation and an accelerated beam of highly relativistic baryons. We recall that we have 
respectively defined the radiation burst (the proper GRB or for short P-GRB) and the accelerated-baryonic-matter 
(ABM) pulse. By computing for a fixed value of the EMBH different PEMB pulses corresponding to selected values of 
B in the range [lO~*-10~^] , it has been possible to obtain a crucial universal diagram which is reproduced in FigEl 
In the limit of B — > 10~® or smaller almost all Edya is emitted in the P-GRB and a negligible fraction is emitted in 
the kinetic energy Esaryons of the baryonic matter and therefore in the afterglow. On the other hand in the limit 
B 10^^ which is also the limit of validity of our theoretical framework, almost all Edya is transferred to Esaryons 
and gives origin to the afterglow and the intensity of the P-GRB correspondingly decreases. We have identified the 
limiting case of negligible values of B with the process of emission of the so called "short bursts" . A complementary 
result reinforcing such an identification comes from the thermodynamical properties of the P-GRB: the hardness of 
the spectrum decreases for increasing values of B, see Fig.lHl 

The injector phase is concluded by the emission of the P-GRB and the ABM pulse, as the condition of transparency 
is reached. 

The beam-target phase, in which the accelerated baryonic matter (ABM) generated in the injector phase collides 
with the ISM, gives origin to the afterglow. Again for simplicity we have adopted a minimum set of assumptions: 

1. The ABM pulse is assumed to collide with a constant homogeneous interstellar medium of number density 
^ism ~ Icm^'^. The energy emitted in the collision is assumed to be instantaneously radiated away (fully 
radiative condition). The description of the collision and emission process is done using spherical symmetry, 
taking only the radial approximation neglecting all the delayed emission due to off-axis scattered radiation. 

2. Special attention is given to numerically compute the power of the afterglow as a function of the arrival time 
using the correct governing equations for the space-time transformations in line with the RSTT paradigm. 

3. Finally some approximate solutions are adopted in order to obtain the determination of the power law exponents 
of the afterglow flux and compare and contrast them with the observational results as well as with the alternative 
results in the literature. 

We first consider the above mentioned radial approximation and a spherically symmetric distribution in order to 
concentrate on the role of the correct space-time transformations in the RSTT paradigm and illustrate their impact 
on the determination of the power law index of the afterglow. This topic has been seriously neglected in the literature. 




Figure 7: a) Afterglow luminosity computed for an EMBH of E^ya = 5.29 x lO'^^ erg, Edya = 4.83 x lO"'^ erg, Edya = 4.49 x 10*'' 
erg and B = 3 x 10~^. b) for the Edya = 4.83 x 10*^, we give the afterglow luminosities corresponding respectively to 
B = 9 X 10~^ 6 x 10~^ 3 X 10"^ 1 X 10^^ ? x 10""*, 4 X 10"*. 



We then turn to the f ully r elativistic analysis of the off-axis emission and of the temporal structure in the lo ng; b ursts 
(see also Ruffini et al. (l48j and sections |XXIHXXII|) and of their spectral distribution (see also Rufhni et al. |l5Cl| | and 
section IXXIV|l . Details of the role of beaming are going to be discussed elsewhere (Ruffini et al. 151]). 
We can now turn to the two eras of the beam-target phase: 

The Era IV: the ultrarelativistic and relativistic regimes in the afterglow. In section IXIIII the hydrodynamic 
relativistic equations governing the collision of the ABM pulse with the interstellar matter are given in the form of 
a set of finite difference equations to be numerically integrated. Expressions for the internal energy developed in the 
collision as well as for the gamma factor are given as a function of the mass of the swept up interstellar material 
and of the initial conditions. In section IXVIIII the infinitesimal limit of these equations is given as well as analytic 
power-law expansions in selected regimes. 

The Era V: the approach to the nonrelativistic regimes in the afterglow. In section IXIVI it is stressed that this last 
era often discussed in the current literature can be described by the same equations used for era IV. 

Having established all the governing equations for all the eras of the EMBH theory, we can proceed to compare and 
contrast the predictions of this theory with the observational data. 



F. The Best fit of the EMBH theory to the GRB 991216: the global features of the solution 

As expressed in section lXVl we have proceeded to the identification of the only two free parameters of the EMBH the- 
ory, Edya and B, by fitting the observational data from R-XTE and Chandra on the decaying part of the GRB 991216 
afterglow. The afterglow appears to have three different parts: in the first part the luminosity increases as a function 
of the arrival time, it then reaches a maximum and finally monotonically decreases. In Fig.d we show how such a fit 
is actually made and how changing the two free parameters affects the intensity and the location in time of the peak 
of the afterglow. The best fit is obtained for Edya = 4.83 x 10^^ erg and B = 3 x 10~^. 

Having determined the two free parameters of the theory, we have integrated the governing equations corresponding 
to these values and then obtained for the first time the complete history of the gamma factor from the moment of 
gravitational collapse to the latest phases of the afterglow observations (see Fig. |H1). This diagram clearly shows the 
inadequacy of considering a simple power-law relation 7 oc r~^/^ for the relation between the radius of the source 
and its L orent z gamm a fac tor as assumed in the large majority of cu rrent papers on GRBs (see e.g. Panaitescu & 
Meszaros jll4l |. Piran |ll(ij . Sari |l63i Il6^ . Sari et al. |l6fil |. Waxman |l8fi| | and references therein). Actually, such a 
power-law behaviour is never found to exist. 

We have also determined the different regimes encountered in the relation between the laboratory time and the 
detector arrival time within the RSTT paradigm compared and contrasted with the ones in the current literature 
(see Fig. EJ . The solid curve is computed using the exact formula prescribed by the RSTT paradigm (see Eq. (pITjl in 
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Figure 8: The theoretically computed gamma factor for the parameter values Edya ~ 4.83 x 10^^ erg, B — 3 x W~'^ is given 
as a function of the radial coordinate in the laboratory frame. The corresponding values in the comoving time, laboratory 
time and arrival time are given in Tab. HI The different eras indicated by roman numerals are illustrated in the text (see 
sections IVIIIVIIIIIXIXIIIIXIVfl . while the points 1,2,3,4,5 mark the be ginning and end of each of these eras. The points Pl 
and Pa mark the maximum of the afterglow flux, respectively in emission time and in arrival time (see RufBni et al. Il45| and 
sections |XIIIIX VIII|l . The point 6 is the beginning of Phase D in Era V (see sections |XI VIX VIII|l . At point 4 the transparency 
condition is reached and the P-GRB is emitted. This diagram clearly shows the inadequacy of considering a simple power-law 
relation 7 oc r~^^^ for the relation between the radius of the sourc e an d its L orent z ga mma facto r as assum ed in the large 
majo rity of current papers on GRBs (see e.g. Panaitescu & Meszaros |114(| . Piran jllfill . Sari |l63lll64| - Sari et al. |l66j| . Waxman 
|186| and references therein). Actually, such a power-law behaviour is never found to exist. 
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The dashed-dotted curve is computed using the approximate formula (see Ea. (|41|l 'l 

^^ = (1 + ^)2^' 

often used in the current hterature (see e.g. Fenimore et al. |4S|, Piran jllfi| |. Sari |l6,li Il64l |. Waxman |l8fi| | and 
references therein). The difference between the solid line and the dashed-dotted line clearly shows the inadequacy of 
using such an approximate relation. We like to stress that the difference between the above two curves is especially 
marked in the afterglow region. Note that this difference as been estimated assuming in both curves the correct relation 
between the Lorentz gamma factor and the radial coordinated of the source given in Fig. |H1 In the cas e tha t the 
wrong relation 7 cx r~^/^ is adopted as done in the literature (see e.g. Panaitescu & Meszaros \ll4\, Piran ban 
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Figure 9: Relation between the arrival time (if) measured at the detector and the laboratory time (t) measured at 
the GRB source. The solid curve is computed using the exact formula prescribed by the RSTT paradigm = 

(1 + z)(^t- /,* ^''l^'^p'^ dt' - ^ j (see Eq.lI33 in section The dashed-dotted curve is computed using the approximate 

formula tt = {1 + z) {tl2'y'^ (t)) (see Ea. l|41|l 'l often used in the current literature (see e.g. Fenimore et al. |4^, Piran |llfil |. Sari 
|l63l Il64|l . Waxman (18Q] and references therein). The difference between the solid line and the dashed-dotted line clearly 
shows the inadequacy of using such an approximate relation. We like to stress that the difference between the above two curves 
is especially marked in the afterglow region. Note that this difference as been estimated assuming in both curves the correct 
relation between the Lorentz gamma factor and the radial coordinated of the source give n in Fig.lHl In the case that the wrong 
relati on 7 oc r~^^^ is adop ted as done in the literature (see e.g. Panaitescu & Meszaros |ll4j . Piran |ll(il |. Sari ,163, 164], Sari 
et al. |166(I . Waxman 'l8(f| and references therein) the discrepancy between the two curves will be much larger. It is anyway 
clear that, even knowing quantitatively the exact Lorentz gamma factor curve reported in Fig. |H1 the use of the approximate 
relation given in Eg. 1411 is enough to miss the correct clock synchronization and to obtain a wrong value for the power- law 
index n in the decaying phases of the afterglow (see sections IXVIIIllXIXI and Tab. ITT)! . We distinguish four different phases. 
Phase A: There is a linear relation between t and ff, given by Ea. ll37j in the text (dashed line). Phase B: There is an 
"effective" power-law relation between t and , given by Eg. 114211 (dotted line). Phase C: No analytic formula holds and the 
relation between t and has to be directly computed by the integration of the complete equations of energy and momentum 
conservation fEgs. 110711081 '). Phase D: As the gamma factor approaches 7 = 1, the relation between t and asymptotically 
goes to t = (light gray line). See also Ruffini et al. |144l| . 

|l63l Il64l | , Sari et al. |l66j , Waxman |l86l | and references therein) the discrepancy between the two curves will be 
much larger. It is anyway clear that, even knowing quantitatively the exact Lorentz gamma factor curve reported in 
Fig. |S1 the use of the approximate relation given in Eq. 1)41(1 is enough to miss the correct clock synchronization and 
to obtain a wrong value for the power-law index n in the decaying phases of the afterglow (see sections IXVIIIl lXIXI 
and Tab.ini). 

To be more explicit, from the result given in Figs. |H1 El follows that all existing GRB models, with the exception of 




Figure 10; Best fit of tlie afterglow data of Chandra, RXTE as well as of the range of variability of the BATSE data on the 
major burst, by a unique afterglow curve leading to the parameter values Edya = 4.83 x lO^^erg, B = 3 x 10""^. The horizontal 
dotted line indicates the BATSE noise threshold. On the left axis the luminosity is given in units of the energy emitted at the 
source, while the right axis gives the flux as received by the detectors. 



ours, have the wrong spacetime coordinatization of the GRB phenomenon and they therefore lack the fundamental 
toola to compare the theoretical prediction in the laboratory time to the observations carried out in the asymptotic 
photon arrival time. This extreme situation affects all considerations on GRBs: as an example, all the considerations 
on the afterglow slopes, which drastically depend on the functional dependence between the laboratory time and the 
photon arrival time, are drastically affected (see subsection III HI below and Tab. ITl)l . In turn, all the considerations 
about the possible existence of beaming in GRBs inferred from the afterglow slopes are in this circumstance deprived 
of any meaning. 

We have thus determined the entire space-time grid of the GRB 991216 by giving (see Tab.P) the radial coordinate 
of the GRB phenomenon as a function of the four coordinate time variables. A quick glance to Tab. ^ shows how 
the extreme relativistic regimes at work lead to enormous superluminal behaviour (up to lO^c!) if the classical 
astrophysical concepts are adopted using the arrival time as the independent variable. In turn this implies that 
any causal relation based on classical astrophysics and the arrival time data, as at times found in the current GRB 
literature, is incorrect. 



G. The explanation of the "long bursts" and the identification of the proper gamma ray burst(P-GRB) 

In section IX VII having determined the two free parameters of the EMBH theory, we analyze the theoretical predic- 
tions of this theory for the general structure of GRBs. The first striking result, illustrated in Fig.^| shows that the 
peak of the afterglow emission coincides both in intensity and in arrival time (19.87 s) with the average emission of 
the long burst observed by BATSE. For this we have introduced the new concept of extended afterglow peak emission 
(E-APE). Once the proper space-time grid is given (see Tab-IJ) it is immediately clear that the E-APE is generated 
at distances of 5 x 10^^ cm from the EMBH. The long bursts are then identified with the E-APEs and are not 
bursts at all: they have been interpreted as bursts only because of the high threshold of the B ATSE detectors (see 
Fig. HOB . Thus the lon g sta nding unsolved problem of explaining the long GRBs (see e.g. Piran |117| . Salmonson et 
al. l|32j, Wilson et al. I l87l |) is radically resolved. 

Still in section IXVII the search for the identification of the P-GRB in the BATSE data is described. This 
identification is made using the two fundamental diagrams shown in Fig. 1111 Having established the value of 
Edya = 4.83 X 10^^ erg and of i? = 3 x 10^'^, it is possible from the dashed line and the solid line in Fig. 1111 to 
evaluate the ratio of the energy Ep-qjib emitted in the P-GRB to the energy Esaryons emitted in the afterglow cor- 
responding to the determined value of B, see the vertical line in Fig.^] We obtain Ep-crb / Esaryons = 1-58 x 10"^, 
which gives Ep-grb = 7.54 x 10^^ erg. Having so determined the theoretically expected intensity of the P-GRB, a 
second fundamental observable parameter, which is also a function of Edya and B, is the arrival time delay between 
the P-GRB and the peak E-APE, determined in Fig.[TTl From Tab.|l| we have that the detector arrival time of the 



Table I: Gamma factors for selected events and their space-time coordinates. The points marked 1,2,3,4,5,6, Pl , Pa are the 
same reported in Fig. |H1 while the point F is the endpoint of the simulation. It is particularly important to read the last 
column, where the apparent motion in the radial coordinate, evaluated in the arrival time at the detector, leads to an enormous 
"superluminal" behaviour, up to 9.55 x 10* c. This illustrates well the impossibility of using such a classical estimate in regimes 
with gamma factors up to 310.1. 
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Figure 11: Left) Relative intensities of the E-APE (dashed line) and the P-GRB (solid line), as predicted by the EMBH theory 
corresponding to the values of the parameters determined in Fig. 1101 as a function of B. Details are given in section DC VII The 
vertical line corresponds to the value B = 3 x 10~^. Right) The arrival time delay between the P-GRB and the peak of the 
E-APE is plotted as a function of the B parameter for three selected values of Edya- 

P-GRB occurs at 8.41 x 10"^ s, corresponding to a radial coordinate of 1.94 x 10^^ cm, a comoving time of 21.57 s, a 
laboratory time of 6.48 x 10^ s and an arrival time of 4.21 x 10^^ s. At this point, the gamma factor is 310.1. The 
peak of the E-APE occurs at a detector arrival time of 19.87 s, corresponding to a radial coordinate of 5.18 x 10^^ cm, 
a comoving time of 5.85 x 10^ s, a laboratory time of 1.73 x 10^ s and an arrival time of 9.93 s (see Tab.CJ. The delay 
between the P-GRB and the peak of the E-APE is therefore 19.78 s, see Fig. 1111 The theoretical prediction on the 
intensity and the arrival time uniquely identifies the P-GRB with the "precursor" in the GRB 991216 (see Fig. ISJ. 
Moreover, the hardness of the P-GRB spectra is also evaluated in this section. As pointed out in the conclusions, 
the fact that both the absolute and relative intensities of the P-GRB and E-APE have been predicted within a few 
percent accuracy as well as the fact that their arrival time has been computed with the precision of a few tenths of 
milliseconds, see Tab. and Fig. El can be considered one of the major successes of the EMBH theory. 

H. On the power-laws and beaming in the afterglow of GRB 991216. 

In section IXVIIII a pi ecew ise description of the afterglow by the expansion of the fundamental hydrodynamical 
equations given by Taub jl77| and Landau & Lifshitz ^^87] have allowed the determination of a power-law index for the 
dependence of the afterglow luminosity on the photon arrival time at the detector. It is evident that the determination 
of the power-law index is very sensitive to the basic assumptions made for the description of the afterglow, as well 
as to the relat ions between the different temporal coordinates which have been clarified by the RSTT paradigm (see 
RufRni et al. |144| ). The different power-law indexes obtained are compared and contrasted with the ones in the 
current literature (see Tab. ITll and section. [XIX|I . As a byproduct of this analysis, see also the conclusions, there is 
a perfect agreement between the observational data and the theoretical predictions, implying that the assumptions 
we have adopted for the description of the afterglow (see section |XIII|) must be necessarily all valid and therefore, in 
particular, there is no evidence for a beamed emission in GRB 991216. 

We then summarize in Fig. El the results for the average bolometric luminosity of GRB 991216 with particular 
attention to the striking agreement, both in arrival time and in intensity, for the theoretically predicted structure of 
the P-GRB and the E-APE with the observational data . To show the generality of application of the EMBH theory, 
we have applied it also to GRB 980425 (see RufRni |l4nj ) and the excellent results are also shown, for comparsion, in 
Fig.ini 

I. Substructures in the E-APE due to inhomogeneities in the Interstellar medium 

In section IXXI the role of the inhomogeneities in the interstellar matter has been analyzed in order to explain the 
observed temporal substructures in the BATSE data on GRB 991216. Having satisfactoril y ide ntified the average 
intensity distribution of the afterglow and the relative position of the P-GRB, in Ruffini et al. |147| we have addressed 
the issue whether the fast temporal variation observed in the so-called long bursts, on time scales as short as fraction of 



Table II: We compare and contrast the results on the power-law index n of the afterglow in the EMBH theory with other 
treatments in the current literature, in the limit of high energy and fully radiative conditions. The differences between the 
values of —10/7 ~ —1.43 (Dermer) and the results —1.6 in the EMBH theory can be retraced to the use of the two different 
approximation in the arrival time versus the laboratory time given in Fig. El See details in section IX Villi 
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a second (see e.g. Fishnian & Meegan l5J|), can indeed be explained as an effect of inhomogeneities in the interstellar 
medium. Such a possibility was pioneered in the work by Dermer & Mitman 42], purporting that such a time 
variability corresponds to a tomographic analysis of the ISM. In order to probe the validity of such an explanation, 
we have first considered the simplified case of the radial approximation (Ruffini et al. |147| ). The aim has been to 
explore the possibility of explaining the observed fluctuation in intensity on a fraction of a second as originated from 
inhomogeneities in ISM, typically of the order of 10^® due to apparent superluminal behaviour of roughly lO^c. We 
have shown there that this approach is indeed viable: both the intensity variation and the time scale of the variability 
in the E-APE region can be explained by the interaction of the ABM pulse with inhomogeneities in the ISM, taking into 
due account the apparent superluminal effects. These effects, in turn, can be derived and computed self consistently 
from the dynamics of the source. We have then described the inhomogeneities of the ISM by an appropriate density 
profile (mask) of an ISM cloud. Of course at this stage, for simplicity, only the case of spherically symmetric "spikes" 
with over-density separated by low-energy regions, has been considered. Each spike has been assumed to have the 
spatial extension of lO^^cm. The cloud average density is < ntsm >= 1 particle/cm^. In conclusion, from the data 
of Tab. m and the highly "superluminal" behaviour of the source in the region of the E-APE, it is concluded that 
the observed time variability in the intensity of the emission (AJ//) ~ 5 can be traced to inhomogeneities in the 
interstellar matter: (Ariism/nism) ^ 5. The typical size of the scattering region is estimated to be 5 x 10^^ cm, 
and these are the typical sizes and density contrasts found in interstellar clouds. Since the emission of the E-APE 
occurs at typical dimensions of the order of 5 x 10^^ cm, the observed inhomogeneities are probing the structure of 
the interstellar medium, and have nothing to do with the "inner engine" of the source. 

The big issue was then open if all these results, obtained in the radial approximation, would still be valid in the 
more general case when off-axis emission in the description of the afterglow is taken into a ccoun t. This is the reason 
why we have proceeded to the topic summarized in the next subsections (see Ruffini et al. |l48j l. 




Figure 12: Left) The overall description of the EMBH theory applied to GRB 991216. The BATSE noise threshold is 
represented and the observations both of the P-GRB and of the E-APE are clearly shown in the subpanels. The continuos 
line in the picture represents the theoretical prediction of the EMBH model. Right) The same diagrams are represented for 
GRB 980425. Two aspects are especially important to be mentioned: a) in this source the theoretical prediction of the P-GRB 
intensity is lower than the BATSE noise treshold and is therefore unobservable and unobserved; b) the E -AP E is especially 
smooth as a consequence of the low value of the gamma Lorentz factor (see also section IXXllII and RufRni |14(1 |'). 



J. The definition of the equitemporal surfaces (EQTS) and the afterglow delayed intensity as a function of 

the viewing angle 

While the analysis of the average bolometric intensity of GRB was going on in the radial approximation, we 
have proceeded to develop the full non-radial approximation, taking into account all the relativistic correction s for 
the off-axis emission from the spherically symmetric expansion of the ABM pulse (see RufRni et al. |l48l Il53j | and 
sections IXXIHXXIljl . Photons emitted at the same time but at different angles of displacement from the line of sight 
reach the detector at very different arrival times. Correspondingly, photons detected at the same arrival time are 
emitted at very different times and angles. We have so defined the temporal evolution of the ABM pulse visible area 
as well as the equitemporal surfaces (EQTS), i.e. the locus of points on the ABM pulse emitting surface corresponding 
to a constant value of the photon arrival time at the detector. 

The very same difficulties found in the current literature, relating the laboratory time to the photon arrival time at 
the detector (see Figs.|Hlini), still exists in the present context and are even magnified in the definition of the EQTS. In 
a classical article, Rees jl27i | expressed the relation between the laboratory time and the arrival time at the detector in 
order to explain observations in radio sources with a constant expansion velocity v and Lorentz gamma factor 7 ~ 5. 
He pointed out the EQTS are ellipsoids of constant eccentricity v/c. In the c urren t liter ature the Rees appr oach 
has been adapted to the analysis of GRBs (see e.g. Fenimore et al. [43|, Piran Sari [Ulliel, Waxman 

and references therein). In addition to the very crucial relation between the laboratory time and the photon arrival 
time, which has not been properly treated, there have been a variety of other approximation and averaging processes 
on which we do not agree. Instead of specifically criticizing each assumption which we consider not correct, such 
comparison will be made in a forthcoming paper (RufSni et al. [l55]), we just report here in the following the results 
of the EQTS surfaces (see Fig. I1H|I obtained in conformity with the RSTT paradigm. In the present case of GRBs, 
the gamma factor is not only much larger than the one observed in radio sources, but is also strongly time varying 
(see Fig. IS)). The Rees treatment has to be significantly improved to take into account the huge time variations in 
the Lorentz gamma factor: this is not just a technical point of modifying a formula by the introduction of a new 




Distance from the EMBH (cm) Distance [mm the EM8H (on) 

Figure 13: Left) This figure shows the temporal evolution of visible area of the ABM pulse. The dashed half-circles are the 
expanding ABM pulse at radii corresponding to diflerent laboratory times. The black curve marks the boundary of the visible 
region. The EMBH is located at position (0,0) in this plot. Again, in the earliest GRB phases the visible region is squeezed 
along the line of sight, while in the final part of the afterglow phase almost all the emitted photons reach the observer. This 
time evolution of the visible area is crucial to the explanation of the GRB temporal structure. Right) Due to the extremely 
high and extremely varying Lorentz gamma factor, photons reaching the detector on the Earth at the same arrival time are 
actually emitted at very different times and positions. We represent here the surfaces of photon emission corresponding to 
selected values of the photon arrival time at the detector: the eqmtemporal surfaces (EQTS). Such surfaces differ from the 
ellipsoids described by Rees in the context of the expanding radio sources with typical Lorentz factor 7 ~ 4 and constant. In 
fact, in GRB 991216 the Lorentz gamma factor ranges from 310 to 1. The EQTSes represented here (solid lines) correspond 
respectively to values of the arrival time ranging from 5 s (the smallest surface on the left of the plot) to 60 s (the largest one 
on the right). Each surface differs from the previous one by 5 s. To each EQTS contributes emission processes occurring at 
different values of the Lorentz gamma factor. The dashed lines are the boundaries of the visible area of the ABM pulse and 
the EMBH is located at position (0, 0) in this plot. Note the different scale on the two axes, indicating the very high EQTS 
"effective eccentricity" . The time interval from 5 s to 60 s has been chosen to encompass the E-APE emission, ranging from 
7 = 308.8 to 7 = 56.84. 



integral. There is in the present context the crucial point expressed in the RSTT paradigm that the relation between 
the laboratory time and the arrival time at the detector is a function of all the the previous Lorentz gamma factors in 
the history of the source since 7=1 (see Fig.lSJ. In the definition of each EQTS, therefore, the entire previous past 
history of the source does concur and the EQTS surfaces become therefore a very refined and sensitive test of the 
correct description of the entire spacetime evolution of the source. In this case, we no longer have ellipsoids of constant 
eccentricity ^. Since the velocity is strongly varying from point to point, we have more complicated surfaces like the 
profiles reported in Fig. El where at every point there will be a tangent ellipsoid of a given eccentricity, but such an 
ellipsoid varies in eccentricity from point to point (see Fig. El and section IXXI|I . Any departure from the correct 
equation of motion strongly alters the EQTS surfaces and accordingly modifies all th e res ults of the integrations based 
on the EQTS surfaces, e.g. the spectral distribution or the afterglow (Ruffini et al. |l54j ). 

Having determined the EQTS surfaces we have computed the observed GRB fiux at selected values of the photon 
arrival time at the detector, taking into due account the delayed contributions at different angles and we have presented 
the results in section IXXIII and Fig.lTH 

We have then recomputed the afterglow emission of GRB 991216 taking into account all the effects due to this 
temporal spreading in the arrival time as well as the ones due to the dependency of the photon Doppler shift on the 
angle of displacement from the line of sight of the emission location (see section IXXIIjl . The result is reported in 
Fig.m 

From now on all the afterglow intensities are estimated using this very complex and extensive numerical program 
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Figure 14: Left) The predicted afterglow curve for GRB 991216 assuming a constant ISM density equal to 1 particle/cm'^ and 
taking into account all the effects due to off-axis emission (solid line). For comparison we plot also the corresponding curve 
obtained in the simple radial approximation (dashed line). We see that this last curve falls sharply to zero when the ABM 
pulse reaches 7 = 1, while the first one has a much smoother behavior due to the time delay in the arrival of the photons 
emitted at large i9. Recall that when 7 tends to 1, the maximum allowed values of i9 tend to 90°. Right) This figure shows 
how the radiation emitted from different angles contributes to the afterglow luminosity. The solid line on the top of the picture 
is the total luminosity as in the previous plots. The other dashed and dotted curves represent the radiation components 
corresponding to selected values of n in Ea. lfTSTji . From the upper to the lower one they corresponds respectively to n = 0, 
n = 0.05A'', n = 0.25A'', n = G.bN , n — N, where in this plot = 200. We can easily see that the radiation emitted at large 
angles (n = A'') is time shifted with respect to that emitted near the line of sight (n = 0). 



which is rooted in all previous history of the source: the general considerations on simple analytic expansion expressed 
in section FX Villi are kept only as an heuristic procedure as a guideline to comprehend these more complex results. 



K. The E-APE temporal substructures taking into account the off-axis emission 

Having determined the EQTS surfaces, we have reconsidered the E-APE temporal substructure taking into due 
account the off-axis emission contribution (see Fig. 1151 and section IXXIIIII . 

We can distinguish two different regimes corresponding respectively to 7 > 150 and to 7 < 150. In the E-APE 
region (7 > 150) the GRB substructure intensities indeed correlate with the ISM inhomogeneities. In this limited 
region (see peaks A, B, C) the Lorentz gamma factor of the ABM pulse ranges from 7 ~ 304 to 7 ^ 200. The 
boundary of t he v i sible region is smaller than the thickness AR of the inhomogeneities (see Figs. ll5tH^ Tab. II VI and 
Ruffini et al. jl48i Il5,'j |). Under these conditions the adopted spherical symmetry for the density spikes is not only 
mathematically simpler but also fully justified. The angular spreading is not strong enough to wipe out the signal 
from the inhomogeneity spike. 

As we descend in the afterglow (7 < 150), a border-line case occurs at peak D where 7 ^ 140. There the visible 
region is comparable to the thickness AR: to fit the observed data a three dimensional description would be necessary, 
breaking the spherical symmetry and making the computation more difficult, but we do not foresee any conceptual 
difficulty. For the peaks E and F we have 7 ^ 50: under these circumstances the boundary of the visible region 
becomes much larger than the thickness AR. The spherically symmetric description of the inhomogeneities is already 
enough to prove the overwhelming effect of the angular spreading and no three dimensional description is needed 
(Ruffini et al. 148, 153.])- 

From our analysis we can conclude that Dermer's expectations do indeed hold for 7 > 150. However, as the gamma 
factor drops from 7 '--^ 150 to 7 ~ 1 the intensity due to the inhomogeneities markedly decreases due to the angular 
spreading (events E and F). The initial Lorentz factor of the ABM pulse 7 ~ 310 decreases very ra pidly to 7 ~ 150 
as soon as a fraction of a typical ISM cloud is engulfed (see Figs. 115151 Tab. II VI and Ruffini et al. |l48l Il53| ). We 
conclude that the "tomography" is indeed effective, but uniquely in the first ISM region close to the source and for 
GRBs with 7 > 150. 

It is then clear that no information on the nature of the GRB source can be inferred by the analysis of the T90, 
nor by the intensity variability structure of the so-called "long hurts" : the only indirect information can be obtained 




I 

f 



\.o> Ifl" 



J 

F 

a s.[b:itr" 
3 

l 

E 




lIL 



r (-rm) 



i W 



"^ypcfkmitiij" 



A 




Lis > li}* 


1.3D !< Ity* 


T.OO 


IS.S 


a.iiio 




B.a ^< m^e 




S JO 10" 


5,74 s 10' 


1.73 X Ity* 


fl-.5tl 


10.0 








C 




10" 


LDD X 10^ 


ILd 


22.B 


1.13 




[1.3 K lO-'c 


D 


^,3t)x It)'* 


7,K4 lif 




15,n 


HO. 1 


S.Hi 




9.B y ItHr 


E 


e,w >: 10'* 


^,32 K IC' 


S.17 St 1(^» 


J7.g 


&5.1? 


10.2 




3.9 M lO^f 


F 


fi.SO J< 10'* 


1.10 y lO* 


1.27 « 10^ 


43,7 


S7.^ 


lD.fi 




2 11 X inh- 



Figure 15: In this figure we summarize the main results of the fit obtained by the EMBH model for the E-APE intensity in 
the case of GRB 991216 taking into account all off-axis contributions. The upper two diagrams represent respectively the 
observational data and the corresponding theoretically computed results. On the lower left the "mask" of the spherically 
symmetric density inhomogeneities with average < Uism >= 1 particle/cm^ is represented. The table summarizes all the 
parameters corresponding to the inhomogeneities including the vary large apparent superluminal effect up to ~ 10^ 
in section 23. 



from the value of Lorentz gamma factor, which has to be 7 > 150 in presence of significant observed substructure. 
In this sense compare and contrast the t wo cases of GRB 991216 and GRB 980425 where the 7 value in the E-APE 
is found to be 7 ^ 120 (see Ruffini |l40j '). The intensity substructures in the E-APE only carry information on the 
structure of the ISM clouds. 



L. The observation of the iron hnes in GRB 991216: on a possible GRB-supernova time sequence 

In section IXXVI the program of using GRBs to further explore the region surrounding the newl y formed EMBH 
is carried one step further by using the observations of the emitted iron hnes (Piro et al., |l2f]l |). This gives us 
the opportunity to introduce the GRB-supernova time sequence (GSTS) paradigm and to introduce as well the novel 
concept of an induced supernova explosion. The GSTS paradigm reads: A massive GRB-progenitor star Pi of mass Mi 
undergoes gravitational collapse to an EMBH. During this process a dyadosphere is formed and subsequently the P- GRB 
and the E-APE are generated in sequence. They propagate and impact, with their photon and neutrino components, 
on a second supernova-progenitor .star P2 of mass M2. Assuming that both stars were generated approximately at 
the same time, we expect to have M2 < Mi. Under some special conditions of the thermonuclear evolution of the 
supernova-progenitor star P2, the collision of the P-GRB and the E-APE with the star P2 can induce its supernova 
explosion. 

Using the result presented in Tab. U and in all preceding sections, the GSTS paradigm is illustrated in the case of 



GRB 991216. Some general considerations on the nature of the supernova progenitor star are also advanced. 

Some general considerations on the EMBH formation are presented in section IXXVII The general conclusions are 
presented in section IXXIXI 

We now proceed to a more detailed presentation of the results and we refer to the already published material for 
the complete details. 



III. THE ZEROTH ERA: THE PROCESS OF GRAVITATIONAL COLLAPSE AND THE FORMATION 

OF THE DYADOSPHERE 

We first recall the three theoretical results which lie at the basis of the EMB H theory. 

In 1971 in the article "Introducing the Black Hole" (Ruffini & Wheeler |l4lj 'l. the theorem was advanced that the 
most general black hole is characterized uniquely by three independent parameters: the mass-energy M , the angular 
momentum L and the charge Q making it an EMBH. Such an ansatz, which came to be known as the "uniqueness 
theorem" has turned out to be one of the most difficult theorems to be proven in all of physics and mathematics. 
The progress in the proof has been authoritatively summarized by Carter The situation can be considered 

satisfactory from the point of view of the physical and astrophysical considerations. Nevertheless some fundamental 
mathematical and physical issues concer ning the most general perturbation analysis of an EMBH are still the topic 
of active scientific discussion (Bini et al. [l^)- 

In 1971 it was shown that the energy extractable from an EMBH is governed by the mass-energy formula 
(Christodoulou & Ruffini [13), 



= = ( M,,c^ + ^ j (1) 



2 \ 2 f2„2 

1 



;4(7r)(Q' + 4i'^')^i' (2) 



with 



where 

S = 4npl = An{rl + = 16^ (^^) M^,, (3) 

is the horizon surface area, Mjr is the irreducible mass, r+ is the horizon radius and p+ is the quasi-spheroidal 
cylindrical coordinate of the horizon evaluated at the equatorial plane. Extreme EMBHs satisfy the equality in Eq.(|21l. 
Up to 50% of the mass-energy of an extreme EMBH can in p rinciple be extracted by a special set of transformations: 
the reversible transformations (Christodoulou & Ruffini 1301'). 

In 1975, generalizing some previous results of Zaumen '191], and Gibbons ^67], Damour & Ruffini "3^ showed that 
the vacuum polarization process a la Heisenberg-Euler-Schwinger (Heisenberg & Euler |78| |. Schwinger _172j ) created 
by an electric field of strength larger than 

can indeed occur in the field of a Kerr-Newmann EMBH. Here and e are respectively the mass and charge of 
the electron. There Damour and Ruffini considered an axially symmetric EMBH, due to the presence of rotation, 
and limited themselves to EMBH masses larger then the upper limit of a neutron star for astrophysical applications. 
They purposely avoided all complications of black holes with mass smaller then the dual electron mass of the electron 



ch 



which may lead to quantum evaporation processes (Hawking 76]). They pointed out that: 



The vacuum polarization process can occur for an EMBH mass larger than the maximum critical mass for 



neutron stars all the way up to 7.2 x lO^Mg 



2. The process of pair creation occurs on very short time scales, typically , and is an almost perfect reversible 
process, in the sense defined by Christodoulou-Ruffini, leading to a very efiicient mechanism of extracting energy 
from an EMBH. 

3. The energy generated by the energy extraction process of an EMBH was found to be of the order of lO^'' erg, 
released almost instantaneously. They concluded at the time "this work naturally leads to a most simple model 
for the explanation of the recently discovered "f-ray bursts". 



After the discovery of the afterglow of GRBs and the determination of the cosmological distance of their sources 
we noticed the coincidence between the theoretically predicted energetics and the observed ones in Damour & Rufhni 
we retu rned to our theore tical resul ts developing some new basic theoretical concepts (Preparata et al. fl2^ 
[l2^ . RufHni [iS^, Ruffini et al. [mKm^), which have led to the EMBH theory. 

As a first simplifying assumption we have developed our considerations in the absence of rotation with spherically 
symmetric distributions. The space-time is then described by the Reissner-Nordstrom geometry, whose spherically 
symmetric metric is given by 



d^s = gtt{r)dH + grr{r)d^r + r^d^O + sit/ ed^(l) , 



(5) 



where gu{r) = - [l - + tS\ = and grr{r) = a-^[r). 

The first new result we obtained is that the pair creation process does not occur at the horizon of the EMBH: it 
extends over the entire region outside the horizon in which the electric field exceeds the critical value given by Eq. ^ 
Since the electric field in the Reissner-Nordstrom geometry has only a radial component given by (see Ruffini jl36j) 



this region extends from the horizon radius 



r+ = 1.47- 10^Ai(l + Vl -C) cm 

out to an outer radius (Ruffini |l37l| ') 



(6) 



(7) 



(8) 



where we have introduced the dimensionless mass and charge parameters /i = -^g^, ^ = ^^^^^^ < 1, see Fig. 0J 

The second new result has been to realize that the local number density of electron and positron pairs created in 
this region as a function of radius is given by 



ne+e-{r) 



Q 



Anr^ -S- e 

V mc/ 

and consequently the total number of electron and positron pairs in this region is 



Q-Qc 



1 



(9) 



(10) 



where Qc — £ci"\- 

The total number of pairs is larger by an enormous factor r*/ (h/mc) > 10^® than the value Q/e which a naive 
estimate of the discharge of the EMBH would have predicted. Due to this enormous amplification factor in the number 
of pairs created, the region between the horizon and r* is dominated by an essentially high density neutral plasma of 
electron-positron pairs. We have defined this region as the dyadosphere of the EMBH from the Greek duas, d uadso s 
for pairs. Consequently we have called r* the dyadosphere radius r* = r^s (Preparata et al. |l2.lLll2^ . Ruffini [iSTll). 
The vacuum polarization process occurs as if the entire dyadosphere are subdivided into a concentric set of shells of 
capacitors each of thickness h/rrieC and each producing a number of e~^e~ pairs on the order of ~ Q/e (see Fig. 0J|. 
The energy density of the electron-positron pairs is given by 



6(r) 



1 - 



(11) 



(see Figs. 2-3 of Preparata et al. [l23j). The total energy of pairs converted from the static electric energy and 
deposited within the dyadosphere is then 
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Figure 16: The energy extracted by the process of vacuum polarization is plotted (solid lines) as a function of the mass M in 
solar mass units for selected values of the charge parameter ^ = 1,0.1, 0.01 (from top to bottom) for an EMBH, the case ^ = 1 
reachable only as a limiting process. For comparison we have also plotted the maximum energy extractable from an EMBH 
(dotted lines) given by eq. Q. Details in Preparata et al. [125|[. 

As we will see in the following this is one of the two fundamental parameters of the EMBH theory (see Fig. I17|l . 
In the hmit ^ ^ 0, Ea. (|12|) leads to i?dya ^ which coincides with the energy extractable from EMBHs by 

reversible processes (Afir — const.), namely Ebh — Mr ~ \^ (Christodoulou & RufRni [sODj see Fig. El Due to the 
very large pair density given by Eq.JSJ and to the sizes of the cross-sections for the process e"'"e~ <-> 7 + 7, the system 
is expected to thermalize to a plasma configuration for which 

ne+ = ne- ^ ^ K+e- > (13) 

where n°_^^ _ is the total number density of e+e^-pairs created in the dyadosphere (see Preparata et al. |l23lll24 ^. 

The third new result which we have introduced for simplicity is that for a given Edya we have assumed either a 
constant average energy density over the entire dyadosphere volume, or a more compact configuration with energy 
density equal to the peak value. These are the two possible initial conditions for the evolution of the dyadosphere 
(see Fig.inj. 

These three old and three new theoretical results permit a good estimate of the general energetics processes origi- 
nating in the dyadosphere, assuming an already formed EMBH. In reality, if the data become accurate enough, the full 
dynamical description of the dyadosphere formation mentioned above will be needed in order to follow all the general 
relativistic eff ects and characteristic ti me s cales of the approach to the EMBH horizon (Cherubini et al. Ruffini 
& Vitaghano |l56lll57f . RufRni et al. |l59j see also section IXXVUl . 

Below we shall concentrate on the dynamical evolution of the electron-positron plasma created in the dyadosphere. 
We shall first examine in the next three sections the governing equations necessary to approach such a dynamical 
description. 
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Figure 17: Left) Selected lines corresponding to fixed values of the Edya are given as a function of the two parameters jj, ^, 
only the solutions below the continuous heavy line are physically relevant. The configurations above the continuous heavy lines 
correspond to unphysical solutions with Vds < r+. Right) Two different approximations for the energy density profile inside 
the dyadosphere. The first one (dashed line) fixes the energy density equal to its peak value, and computes an "effective" 
dyadosphere radius accordingly. The second one (dotted line) fixes the dyadosphere radius to its correct value, and assumes an 
uniform energy density over the dyadosphere volume. The total energy in the dyadosphere is of course the same in both cases. 
The solid curve represents the real energy density profile. 



IV. THE HYDRODYNAMICS AND THE RATE EQUATIONS FOR THE PLASMA OF e+e -PAIRS 



The evolution of the e+e^-pair plasma generated in the dyadosphere has been treated in two papers (Ruffini et 
al. ^42, 143J). We recall here the basic governing equations in the most general case in which the plasma fluid is 
composed of e+e^-pairs, photons and baryonic matter. The plasma is described by the stress-energy tensor 



T^'' =re^'' + (p + p)C/^[/^ 



(14) 



where p and p are respectively the total proper energy density and pressure in the comoving frame of the plasma fluid 
and is its four-velocity, satisfying 



where and are the radial and temporal contravariant components of the 4- velocity. 

The conservation law for baryon number can be expressed in terms of the proper baryon number density ub 

= {nBU')^t + ^{r''nBU^)^r = . 
The radial component of the energy-momentum conservation law of the plasma fluid reduces to 



(15) 
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(16) 



(17) 



The component of the energy-momentum conservation law of the plasma fluid equation along a flow line is 

1 



p 



= 



(18) 



Defining the total proper internal energy density e and the baryonic mass density ps m the comoving frame of the 
plasma fluid, 



PB, 



Pb = riBrnc 



(19) 



and using the law (|16|l of baryon-number conservation, from Eq. (|18() we have 

{eUn^.+piUn^.^O ■ (20) 

Recalling that ^ = y(C/^);^, where y is the comoving volume and r is the proper time for the plasma fluid, we have 
along each flow line 

d(Ve) dV dE dV ^ 
dr ar dr dr 

where E = Ve is the total proper internal energy of the plasma fluid. We express the equation of state by introducing 
a thermal index T{p, T) 

r = 1 + ^ . (22) 

We now turn to the second set of governing equations describing the evolution of the e~ pairs. Letting Uf.- 
and be the proper number densities of electrons and positrons associated with pairs and n''_ the proper number 
densities of ionized electrons, we clearly have 

rie- = ne+ = npair, n\-^ZnB, (23) 

where ripair is the number of e+e^ pairs and Z the average atomic number i<Z<l(Z = l for hydrogen atom and 
Z = ^ for general baryonic matter). The rate equation for electrons and positrons gives, 

1 

mj[{n,-{T)+nl^{T))n,+ {T) 

[n,- +nl-)n,+ ], (24) 
1 

= CTW K-(T)ne+(r) - ne-?ie+] , (25) 

= av [n\- {T)n,^ (T) - ] , (26) 

where Wv is the mean of the product of the annihilation cross-section and the thermal velocity of the electrons 
and positrons, rte±(r) are the proper number densities of electrons and positrons associated with the pairs, given 
by appropriate Fermi integrals with zero chemical potential, and n''„ (T) is the proper number density of ionized 
electrons, given by appropriate Fermi integrals with non-zero chemical potential /^e at an appropriate equilibrium 
temperature T. These rate equations can be reduced to 

{n,±U^'),^ = (n,±[/*)^t-f -^(r^neit/^),. 

= w[n,-(T)n,+ (r)-ne-n,+ ], (27) 
{n^^U^),^ = {nl_U')^t + ^{r^nl^U^),r^Q, (28) 

Frac ^ ^ = (29) 
Equation (|28|l is just the baryon-number conservation law and (|29|l is a relationship satisfied by rig± , n^± (T) and 

The equilibrium temperature T is determined by the thermalization processes occurring in the expanding plasma 
fluid with a total proper energy density p governed by the hydrodynamical equations H16I17I18|I . We have 

P = Pi + Pe+ + Pe- + Pe- + PB, (30) 

where is the photon energy density, ps — msc^riB is the baryonic mass density which is considered to be non- 
relativistic in the range of temperature T under consideration, and p^± is the proper energy density of electrons and 
positrons pairs given by 

Pe± = , Pe±{T), (31) 



where n^i is obtained by integration of Eq.|57J and Pi^±{T) is the proper energy density of electrons(positrons) 
obtained from zero chemical potential Fermi integrals at the equilibrium temperature T. On the other hand p^^ is 
the energy density of the ionized electrons coming from the ionization of baryonic matter 

n'' 

pl- = ^fUn (32) 

where n^- is obtained by integration of Ea. (|28(l and (T) is the proper energy density of ionized electrons obtained 
from an appropriate Fermi integral of non-zero chemical potential fj,e at the equilibrium temperature T. 

Having intrinsically defined the equilibrium temperature T in Ea. (|30|) . we can also analogously evaluate the total 
pressure 

P=P-f+Pe++Pe-+p''e-+PB, (33) 

where is the photon pressure, Pe± and are given by 

Pe± = "^"^ . Pe±{T), (34) 



n 



Pe-iT), (35) 



the pressures p^±{T) are determined by zero chemical potential Fermi integrals, and p^-{T) is the pressure of the 
ionized electrons, evaluated by an appropriate Fermi integral of non-zero chemical potential pe at the equilibrium 
temperature T. In Ea. (|33|l . the ion pressure pB is negligible by comparison with the pressures p^ e± e-(T'), since 
baryons and ions are expected to be nonrelativistic in the range of temperature T under consideration. Finally using 
Eas. H3()l33|l we compute the thermal factor F of the equation of state (|^ . 

It is clear that the entire set of equations considered above, namely Eas- ljltillTllHI) with equation of state given by 
Ea. (|22|l and the rate equation H27|l . have to be integrated satisfying the total energy conservation for the system. The 
boundary conditions adopted here are simply purely ingoing conditions at the horizon and purely outgoing conditions 
at radial infinity. The calculation is initiated by depositing a proper energy density Hll|l between the Reissner- 
Nordstrom horizon radius r+ and the dyadosphere radius r^s, following the approximation presented in Fig lltil The 
total energy deposited is given by Ea. H12() . 



V. THE EQUATIONS LEADING TO THE RELATIVE SPACE-TIME TRANSFORMATIONS 

In order to relate the above hydrodynamic and pair equations with the observations we need the governing equations 
relating the comoving time to the laboratory time corresponding to an inertial reference frame in which the EMBH 
is at rest and finally to the time measured at the detector, which must also include the effect of the cosmological 
expansion. Th ese t ransformations have been the object of the Relative space-time Transformations (RSTT) Paradigm, 
(Ruffini et al. [Til). 

For signals emitted by a pulse moving with velocity v in the laboratory frame (see also Ruffini et al. |l44j l. we have 
the following relation between the interval of arrival time At^ and the corresponding interval of laboratory time At 
(see Fig. CHI: 

Ma = {to + At + -{t^ + ^^^M-"-. (36) 

For simplicity in what follows we indicate by ta the interval of arrival time measured from the reception of a light 
signal emitted at the onset of the gravitational collapse. Analogously, t indicates the laboratory time interval measured 
from the time of the gravitational collapse. In this case, Eg . H36|l can be written simply as: 

t.^t-!:^t-^°-(^')-^^' + --^t- r^g|pdt--!i^, (37) 

c c Jo 7 (^ ) c 



where, as usual, 7 (t') = — v"^ (f) jc? and the dyadosphere radius r^s is the value of r at t = 0. It is important 

to stress that, although there is the presence of the Lorentz gamma factor, Eq.JJTJ is not a Lorentz transformation. 
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Figure 18: This qualitative diagram illustrates the relation between the laboratory time interval At and the arrival time interval 
Ata for a pulse moving with velocity v in the laboratory time (solid line). We have indicated here the case where the motion 
of the source has a nonzero acceleration. The arrival time is measured using hght signals emitted by the pulse (dotted lines). 
Ro is the distance of the observer from the EMBH, to is the laboratory time corresponding to the o nset of the gravitational 
collapse, and r is the radius of the expanding pulse at a time t = to + At. See also Rufhni et al. |144| |. 

which by its own nature is Unear and refers to a specific value of the Lorentz gamma factor at a given laboratory 
time. The transformation in Ea. (|37|) is nonlinear in the Lorentz gamma factor and do depend on all the values of 
the gamma factor of the source from the time t = to the laboratory time t. This transformation is the price to 
pay to relate the laboratory time t, relativistically correct, to the "highly pathological" time usually considered by 
the astronomers, even in the case of object moving close to the speed of light, against the correct synchronization 
procedures established by Einstein in his classical paper of 1905 (Einstein |43). We consider here only the photons 
emitted along the line of sight from the external surface of the pulse. The arrival time spreading d ue to the a ngular 
dependence and that due to the thickness has also been given (see section IXXII and Ruffini et al. |l48i lisij ). The 
solution of Eq. H37|l has the expansion: 

c c 2 c 3 c 

so the relation between ta and t in the specific case of GRBs is very highly nonlinear: it is sufficient to recall that in 
the early GRB phases we are witnessing the strongest acceleration ever recorded in the universe, since the PEM pulse 
goes from Lorentz factor 7 = 1 to Lorentz factor 7 — 1000 in 10^ seconds in the laboratory time (sec scction lVIll . The 
series in Ea. (|38|l will definitely converge, but the number of terms needed to reach a good approximation will strongly 
depend on the variability of the functions around the initial values 7 = 1. It is clear that the precise knowledge of ta 
as a function of the laboratory time, which is indeed essential for any physical interpretation of GRB data, depends 
on the definite integral given in Eq. H37|) whose limits in the laboratory time extend from the onset of the gravitational 
collapse to the time t relevant for the observations. Such an integral depends on all previous values of the Lorentz 
gamma factor in the history of the source and is not generally expressible by a simple linear relation or even by any 
explicit analytic relation since we are dealing with processes with variable gamma factor unprecedente d in the entire 
realm of physics (see Figs.|Hland Fig.O. This is the crucial point of the RSTT paradigm (Ruffini et al. |l44j ) and this 
is the reason why we have spent a very large amount of work to develop the exact equations of motion of all different 
eras of the GRB phenomenon, starting from the onset of gravitational collapse and the creation of dyadosphere (see 
the following sections). It is clear then that, in order to express the arrival time ta and the radial coordinate of 
the source at the start of the afterglow phase, we need the explicit knowledge of all the previous eras of the GRB 



phenomenon, starting from 7=1 (Ruffini et al. [14J|). 

What has been currently done in the hterature, is an extremely different approach. First they have assumed 7 
constant. Therefore Ea. l37|l has been modified in: 



i^ = t_ VZEi f dt'-'-^^t-^^t^ (39) 
1 Jo c 7 

where in the last approximation the contribution of the initial size of the source has been neglected. Even the validity 
of this last approximation has to be actually carefully verified since it is only valid in the late phases of the GRB 
expansion. They have further assumed 7 3> 1 and obtained: 

272 j * ^ 272 



^•^-^-ll-^ (40) 



At this stage, they emphasize the existence of a linear relation between the arrival time ta and the laboratory time 

t. After this they proceed in two dif feren t directions; One to assume (see e.g. Fenimore et al. Fenimore 

m, Fenimore et al. [13 , Sari & Piran Waxman [l8|) 

ta = t/ (272 (t)) , (41) 

concurrently advancing the belief that the relation between the arrival time and the laboratory time does not depend 
from an integral on all the previous values of the gamma Lorentz factor of the source but from the instantaneous 
value of the gamma Lorentz factor at the time t, much like in a Lorentz transformation. This claim is clearly absurd 
from a physical point of view. 

They further assume (see e.g. Panaitescu & Meszaros |113| . Piran |ll6j |. Sari |l63lll64 and references therein) 

5ta = 5t/ (27^ it)) or, alternatively, dta = dtj (27^ (i)) , (42) 

and they proceed to develop all the observable quantities of the GRB phenomenon by integrating using the "differ- 
ential" given in Ea. l42|l . reaching clearly meaningless results. As we show later, this also leads to the unfortunate 
attempt to obtain the gamma Lorentz factor and its time variability from the astrophysical data of the afterglow, 
neglecting all previous GRB source history what is clearly physically and astrophysically impossible. 

Having established the correct relations between the laboratory time t and the arrival time ta in Eq.|j23, we now 
proceed to relate the time in the laboratory frame t to the time in the detector frame if. We have to do one additional 
step: the two frames are related by a transformation which is a function of the cosmological expansion. We recall 
that the geometry of the space-time of the universe is described by the Robertson- Walker metric: 

/ dr^ \ 

ds^ = dt"^ - 7^2(^) — + r^di?^ + r^sind^difi^ , (43) 

\ 1 — fcr^ / 

where TZ (t) is the cosmic scale factor and fc is a constant related to the curvature of the three-dimensional space 
{k — 0, +1,— 1 corresponds to flat, close and open space respectively). The wavelength of an electromagnetic wave 
traveling from the point Pi {ti , ri , i3i , ) to the point Po {to ,ro,'do,^o) where the observer is located is related to the 
red-shift parameter z by 

z = , (44) 

where Ao is the wavelength of the radiation for the observer and Ai for the emitter. We have the following general 
relation: 

l + z = (l + z„)(l + Zo)(l + ^.), (45) 

where z is the total redshift due to the motion of the source Zg, the motion of the observer Zo and the cosmological 
redshift Zu- In the following we will assume Zq << I and Zg << 1 so z = Zu- In terms of the scale factor TZ{t) the 
relation H44|l gives 

Al iZ{ti) ujQ 



where lui and loq are the frequencies associated to Ai and Aq respectively. This frequency ratio then relates the time 
elapsing at the source with the time elapsing at the detector due to the cosmological expansion. 



We can now define the corrected arrival time measured at the detector, which is related to ta, clearly defined by 
Eq.ljSH), by 

tt^ta{l + z), (47) 

where z is the cosmological redshift of the GRB source. In the case of GRB 991216 we have z ~ 1.00. 

The observed flux is the flux which crosses the surface 4:Tt{TZ (to) r)^ but this flux is lower by a factor 1 + z due to 
the redshift energy of the photons and by another factor 1 + z due to the fact that the number of photons at reception 
is less than the number at emission. Thus we can define a luminosity distance by: 

dl=nyil + zy. (48) 

Then the observed flux is related to the absolute luminosity of the GRB by the following relation: 

where the luminosity distance is simply related to the proper distance dp = TZo^ by d^ — dp{l + z). The observed 
total fluence / is related to the total energy E of the GRB by the following relation: 

Then the cosmological effect is taken into account by the definition of the proper distance TZof which depends 
on the cosmological parameters: the Hubble constant Ho ~ TZ{to) /TZ{to) at time to and the matter density po or 

Po/pcrit, where pcnt = 
The computation of the proper distance is then simply given by the relation : 



where F{z) = yj9.M{'^ + zY- 

In the case of the Friedman fiat universe, — ^ and we have: 



dpiz) = -FT 
rlr, 



1 

1 - 



(52) 



So the measurement of the redshift gives us the luminosity distance via a cosmological scenario. With the measure- 
ment of the fiux we can deduce the proper luminosity of the burst and from the measurement of the total fluence the 
total energy so we are then able to find the E^ya- 



VI. THE NUMERICAL INTEGRATION OF THE HYDRODYNAMICS AND THE RATE EQUATIONS 

A. The Livermore code 



A computer code (Wilson et al. j l88l Il89l |'l has been used to evolve the spherically symmetric general relativistic 
hydrodynamic equations starting from the dyadosphere (Ruffini et al. jl42l |). 

We define the generalized gamma factor 7 and the radial 3-velocity in the laboratory frame 



From Eas.l(5l [T5|l. we then have 



Following Ea. p9|l . we also define 



7^v/TTT7^, v"-^^. (53) 



{Uy = - — {l + grr{Un')^^^'. (54) 
gtt a-^ 



i? = £7, D = pbJ, and p= pj 



(55) 



so that the conservation law of baryon number IjlGf) can then be written as 



dD 

'dt 



or " 



a 



Ea. (|18|l then takes the form, 



dE a d r^ 
ot r-^ or a 



ot or a 



Defining the radial momentum density in the laboratory frame 

Sr = a{j) + p)U*Ur = {D + TE)Ur, 
we can express the radial component of the energy- momentum conservation law given in Eq. H17|l by 



(56) 



(57) 



(58) 



dSr a d r^ dp 

at r-^ Or a or 
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dgr 



dr 



cx d r^ 
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dp 
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D + TE\ 
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r\2 



(59) 



In order to determine the number-density of e+e~ pairs, we turn to Ea. (|27|l . Defining the e+e~-pair density in the 
laboratory frame Ng± = 7ne± and Nf.±{T) = 7ng±(T), where the equilibrium temperature T has been obtained from 
Eas. ipnji and ||2U, and using Ea. (|^ . we rewrite the rate equation given by Ea. l(?7jl in the form 



dt 



(60) 



These equations are integrated starting from the dyadosphere distributions given in Fig. 1171 and assuming as usual 
ingoing boundary conditions on the horizon of the EMBH. 



B. The Rome code 



In the follo wing we recall a zeroth order approximation of the fully relativistic equations of the previous section 
(Ruffini et al. [m^): 

(i) Since we are mainly interested in the expansion of the e+e^ plasma away from the EMBH, we neglect the 
gravitational interaction. 

(ii) We describe the expanding plasma by a special relativistic set of equations. 

(iii) In contrast with the previous treatment where the evolution of the density profiles given in Fig. El are followed 
in their temporal evolution leading to a pulse-like structure, selected geometries of the pulse are a priori adopted and 
the correct one validated by the complete integration of the equations given by the Livermore codes. 

Analogously to Ea. (|21|l . from Ea. l|16|) we have along each fiow line in the general case in which baryonic matter is 
present 

%^ = 0. (61) 
dr 

For the expansion of a shell from its initial volume AVo to the volume AT^, we obtain 

_ Ay _ AV7(r) 
" AK ~ AVo7o(r)' 



(62) 



where AV is the volume of the shell in the laboratory frame, related to the proper volume Ay in the comoving frame 
by Ay = 7(r)AV, where 7(r) defined in Ea. (|53|l is the gamma factor of the shell at the radius r. 



Similarly from Ea. (|21|l . using the equation of state H22|l . along the flow lines we obtain 

d\ne + Td\nV = 0. (63) 
Correspondingly we obtain for the internal energy density e along the flow lines 



AVo) \AvJ ^7o(r) 



(64) 



where the thermal index T given by (|22(l is a slowly- varying function with values around 4/3. It can be computed for 
each value of e,p as a function of AV. 

The overall energy conservation requires that the change of the internal proper energy of a shell is compensated by 
a change in its bulk kinetic energy. We then have (Ruflini et al. jl42| ) 

dK = [7(r) - l]{dE + pBdV). (65) 

In order to model the relativistic expansion of the plasma fluid, we assume that E and D as defined by Ea. H55|l are 
constant in space over the volume AV . As a consequence the total energy conservation for the shell implies (Rufhni 
et al. ill) 

(eo + p^)7o'('^)AVo = (e + PB)l\r)AV, (66) 

which leads the solution 



Corresponding to Ea. H6U|l we obtain the equation for the evolution of the number-density as seen by an observer 
in the laboratory frame 

I (iV.± ) . - AT^, A. ^ + av-L-^ [Nl^ (T) - Nl^ ) . (68) 

Eas. H62|l . H64(l . H67|l and Hfi8|l are a complete set of equations describing the relativistic expansion of the shell. If 
we now turn from a single shell to a finite distribution of shells, we can introduce the average values of the proper 
internal-energy, baryon-mass, baryon-number and pair-number densities (e, ps, n^, ngi) and E = je, D = ^ps, 
Nf,± = ^{r)nf,± for the PEM-pulse, where the average 7-factor is defined by 

l=^J^l{r)dV, (69) 

an d V i s the total volume of the shell in the laboratory frame. The corresponding equations are given in Ruffini et 
al. '142]. Having defined all its governing equations we can now return to the description of the different eras of the 
GRB phenomena. 

VII. THE ERA I: THE PEM PULSE 

We have assumed that, following the gravitational collapse process, a region of very low baryonic contamination 
exists in the dyadosphere all the way to the remnant of the progenitor star. 

Recalling Eq.@ the limit on such baryonic contamination, where pB^ is the mass-energy density of baryons, is 
given by 



1- ( — 

rds 



2 

r ' 



(g/cm^). (70) 



PB^ « mpne+e-(r) = 3.2 • 10« (^)^ 
Near the horizon r ~ r_|_, this gives 

PB, < mpTie+e- (r) = 1.86 • lO" ( ^ ) (g/cm^) , (71) 



and near the radius of the dyadosphere rds '■ 



Pb^ ^ mpHe+e- {r) = 3.2 • 10^ 



(72) 



Such conditions can be easily satisfied in the collapse to an EMBH, but not necessarily in a collapse to a neutron star. 

Consequently we have solved the equations governing a plasma composed solely of e"'"e~-pairs and electromagnetic 
radiation, starting at time zero fro m th e dyadosphere configurations corresponding to constant density in Fig. 1171 
The Livermore code (Ruffini et al. jl42l |) has shown very clearly the self organization of the expanding plasma in 
a very sharp pulse which we have defined as the pair-electromagnetic pulse (PEM pulse), in analogy with the EM 
pulse observed in nuclear explosions. In order to further examine the structure of the PEM pulse with the simpler 
procedures of the Rome codes we have assumed (Ruffini et al. |M2j) three alternative patterns of expansion of the 
PEM pulse on which to try the simplified special relativistic treatment and then compared the results with the fully 
general relativistic hydrodynamical results: 



• Spherical model: we assume the radial component of the four-velocity Ur{r) = U^, where U is the radial 
component of the four-velocity at the moving outer surface r = TZ{t) of the PEM pulse and the 7-factor and the 
velocity Vr are 



3 

8LP 



2[/(l + C/2)2 -Uil + U^' 



f r — _ ; 

7 



lii{U + \/l + lP) 



(73) 



this distribution expands keeping an uniform density profile which decreases with time similar to a portion of a 
Friedmann Universe. 

• Slab 1: we assume U{r) — Ur = const., the constant width of the expanding slab T) = Ro '\n the laboratory 
frame of the PEM pulse, while 7 and Vr are 



Ur 



1 



(74) 



this distribution does not need any averaging process. 



• Slab 2: we assume a constant width R2 — Ri = Ro of the expanding slab in the comoving frame of the PEM 
pulse, while 7 and Vr are 



7 = v/TTc^, 



v-^ 

7 



(75) 



This distribution needs an averaging procedure and i?i < f < i?2, i.e. f is an intermediate radius in the slab. 

These different assumptions lead to three different distinct slopes for the monotonically increasing 7-factor as a 
function of the radius (or time) in the laboratory frame, having assumed for the energy of dyadosphere E^ya — 3.1 x 10'^'* 
erg (see Fig. I19II . In principle, we could have an infinite number of models by defining arbitrarily the geometry of 
the expanding fluid in the special relativistic treatment given above. To find out which expanding pattern of PEM 
pulses is the physically realistic one, we need to compare and contrast the results of our simplified models (performed 
in Rome) with the numerical results based on the hydrodynamic Eas. (|56l57l59() (obtained at Livermore) (Ruffini et 
al. Details of the iterative method used to solve the special relativistic equation can be found in Ruffini et al. 

It is manifest from the results (see Fig. I19|l that the slab 1 approximation (constant thickness in the laboratory 
frame) is in excellent agreement with the Livermore results (open squares). 

The remarkable validation of the special relativistic treatment of the PEM pulse (Ruffini et al. |l42l |'l. allows us to 
easily estimate the related quantities of physical and astrophysical interest in the model, like the e^e~-pair densities 
as a function of the laboratory time, the temperature of the plasma in the comoving and laboratory frames, the 
reheating ratio as a function of the e^e~-pair annihilation for a variety of initial conditions (Ruffini et al. (142!|). 
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Figure 19: gamma factor as a function of radius. Three models for the expansion pattern of the PEM-pulse are compared with 
the results of the one dimensional hydrodynamic code for an energy of dyadosphere Edya = 3.1 x 10^* erg. The 1-D code has 
an expansion pattern that strongly resembles that of a shell with constant thickness in the laboratory frame. 



VIII. THE ERA II: THE INTERACTION OF THE PEM PULSE WITH THE REMNANT OF THE 

PROGENITOR STAR 



The PEM pulse expands initially in a region of very low baryonic contamination created by the process of gravita- 
tional collapse. As it moves further out the baryonic remnant (see Fig. ^ of the progenitor star is encountered. As 
discussed in section FXX VII below . the existence of such a remnant is necessary in order to guarantee the overall charge 
neutrality of the system: the collapsing core has the opposite charge of the remnant and the system as a whole is 
clearly neutral. The numb er of extra charges i n th e baryonic remnant negligibly affects the overall charge neutrality 
of the PEM pulse (Ruffini [Hi, Ruffini et al. [HI). 

The baryonic matter remnant is assumed to be distributed well outside the dyadosphere in a shell of thickness A 
between an inner radius rin and an outer radius Tout = + A at a distance from the EMBH at which the original 
PEM pulse expanding in vacuum has not yet reached transparency. For the sake of an example we choose 

Tin = lOOrds, A = lOrd.. (76) 

The total baryonic mass Mb = Nsmp is assumed to be a fraction of the dyadosphere initial total energy (E'dya)- The 
total baryon-number Nb is then expressed as a function of the dimensionless parameter B given by 



NBrUpC^ 

Edya 



B = il^^ , (77) 



where S is a parameter in the range 10^* — 10^^ and is the proton mass. We shall see below the paramount 
importance of B in the determination of the features of the GRBs. We will see in section |3 the sense in which B and 
Edya can be considered to be the only two free parameters of the EMBH theory for the entire GRB family, the so 
called "long bursts". We shall see in section IXIII that for the so called "short bursts" the EMBH theory depends on 
the two other parameters fj,, ^, since in that case B ^ 0. The baryon number density is assumed to be a constant 



-o J^-B -o -o 2 

Ub = — , Pb = mpHgC . 



(78) 



As the PEM pulse reaches the region Tin < r < Tout, it interacts with the baryonic matter which is assumed to be 
at rest. In our simphfied quasi-analytic model we make the following assumptions to describe this interaction: 

• the PEM pulse docs not change its geometry during the interaction; 

• the collision between the PEM pulse and the baryonic matter is assumed to be inelastic, 

• the baryonic matter reaches thermal equilibrium with the photons and pairs of the PEM pulse. 

These assumptions are valid if: (i) the total energy of the PEM pulse is much larger than the total mass-energy of 
baryonic matter Mb, 10^^ < B < 10~^, (ii) the ratio of the comoving number density of pairs and baryons at the 
moment of collision n^+^- /n°g is very high (e.g., 10^ < ng+g-/n^ < 10^^) and (iii) the PEM pulse has a large value 
of the gamma factor (100 < 7). 

In the collision between the PEM pulse and the baryonic matter at rout > f > , we impose total conservation 
of energy and momentum. We consider the collision process between two radii r2,ri satisfying rout > r2 > ri > rju 
and r2 — ri <C A. The amount of baryonic mass acquired by the PEM pulse is 

AM = ^f(ri-r?), (79) 

where Mb/Vb is the mean-density of baryonic matter at rest. The conservation of total energy leads to the estimate 
of the corresponding quantities before (with "o") and after such a collision 

(Peo + P°b)iIVo + AM - (Pe + + ^ + rAe)7'V, (80) 

where Ae is the corresponding increase of internal energy due to the collision. Similarly the momentum-conservation 
gives 

(Peo -f P°b)1oU:Vo = {Te + PB + ^+ V^e)^UrV, (81) 



where the radial component of the four- velocity of the PEM pulse is U° — a/t^— T and P is the thermal index. We 
then find 



, -o, 7o^°Vo ^ AM, 



(82) 



a ^lo ^ AM 

7 = , a= ^- ^ . (83) 

' U° (Peo + p?j)7oC/°Vo ^ ' 

These equations determine the gamma factor 7 and the internal energy density e — lo + Ae in the capture process of 
baryonic matter by the PEM pulse. 

The effect of the collision of the PEM pulse with the remnant leads to the following results (RufRni et al. |l43j l as 
a function of the B parameter defined in Ea. fTTjl : 

1) an abrupt decrease of the gamma factor given by 

where 70 is the gamma factor of the PEM pulse prior to the collision and B is given by Ea. fTTjl . 

2) an increase of the internal energy in the comoving frame Ecoii developed in the collision given by 



Edya 7o V 7o 

3) a corresponding reheating of the plasma in the comoving frame but not in the laboratory frame, an increase of the 
number of e+e~ pairs and correspondingly an overall increase of the opacity of the pulse. See details in section IXTI 



IX. THE ERA III: THE PEMB PULSE 



After the engulfment of the baryonic matter of the remnant the plasma formed of e+e^-pairs, electromagnetic 
radiation and baryonic matter expands again as a sharp pulse, namely the PEMB pulse. The calculation is continued 
as the plasma fluid expands, cools and the e^e~ pairs recombinc until it becomes optically thin: 

dr{n^± + ZnB)<TT ^ 0(1), (86) 



R 



where ax = 0.665 • 10~^*cm^ is the Thomson cross-section and the integration is over the radial interval of the PEMB 
pulse in the comoving frame. We have first explored the general problem of the PEMB pulse evolution by integrating 
the general relativistic hydrodynamical equations with the Livermore codes, for a total energy in the dyadosphere of 
3.1 X 10^^ erg and a baryonic shell of thickness A = lOrds at rest at a radius of lOOrds and B ~ 1.3 • 10~^. 

In t otal analogy with the special relativistic treatment for the PEM pulse, presented in section Ivnl fsec also Rufhni 
et al. obtain for the adiabatic expansion of the PEMB pulse in the constant-slab approximation described 

by the Rome codes the following hydrodynamical equations with ps ^ 
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(87) 
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(89) 



^(iVe±) - -N,i^^ + av^{N^^{T)-N^^). (90) 

In these equations (r > Tout) the comoving baryonic mass- and number densities are pB = Mb/V and ub = Nb/V, 
where V is the comoving volume of the PEMB pulse. 

We compare and contrast (see Fig. [SJ the bulk gamma factor as computed from the Rome and Livermore codes, 
where excellent agreement has been found. This validates the constant-thickness approximation in the case of the 
PEMB pulse as well. On this basis we easily estimate a variety of physical quantities for an entire range of values of 
B. 

For the same EMBH we have considered five different cases: a shell of baryonic mass with (1) B ~ 1.3 • 10^^; (2) 
B ~ 3.8 • 10"''; (3) B ^.3 • lO'^ (4) B ~ 3.8 • lO'^ (5) B ~ 1.3 ■ IQ-^). The results of the integration given in 
detail in RufFini et al. jl4,'l | show that for the first parameter range the PEMB pulse propagates as a sharp pulse of 
constant thickness in the laboratory frame, but already for B ~ 1.3 • 10~^ the expansion of the PE MB p ulse becomes 
much more complex and the constant-thickness approximation ceases to be valid; see Ruffini et al. |l43f for details. 

It is particularly interesting to evaluate the final value of the gamma factor of the PEMB pulse when the transparency 
condition given by Ea. lj86|l is reached as a function of B, see Fig. 1201 For a given EMBH, there is a maximum value of 
the gamma factor at transparency. By further increasing the value of B the e ntire Edya is transferred into the kinetic 
energy of the baryons; see also section IXIII Details are given in Ruffini et al. |l43l |. 

In Fig.EOIwe plot the gamma factor of the PEMB pulse versus the radius for different amounts of baryonic matter. 
The diagram extends to values of the radial coordinate at which the transparency condition given by Ea. H86(l is 
reached. The "asymptotic" gamma factor 

7asym = (91) 



is also shown for each curve. The closer the gamma value approaches the "asymptotic" value (I91|l at transparency, 
the smaller the intensity of the radiation emitted in the burst and the larger the amount of kinetic energy left in the 
baryonic matter. 



X. THE IDENTIFICATION OF THE FREE PARAMETERS OF THE EMBH THEORY 



Within the approximation presented in section ITTll the EMBH is characterized by two parameters: p and ^. The 
energy of the dyadosphere is expressed in terms of these two parameters by Eq. (|12|l . 
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Figure 20: Left) The gamma factor (the soUd Une) at the transparent point is plotted as a function of the B parameter. The 
asymptotic value (the dashed line) Edy?,/ {Mb(?) is also plotted. Right) The gamma factors are given as functions of the radius 
in units of the dyadosphere radius for selected values of B for the typical case Edya = 3.1 x 10^* erg. The asymptotic values 
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r /rds = 100 where the jump occurs and the PEMB pulse starts. 



Edya./ {Mbc^) = 10*, 10^, 10^ are also plotted. The collision of the PEM pulse with the baryonic remnant occurs at 



There is an entire family of EMBH solutions with different values of /i and ^ corresponding to the same value of Edya 
(see Fig. ll7|) . These solutions are physically different with respect to the density of electron-positron pair distributions 
given by Eq.(|5J, as well as to their energy density given by Ea. (|ll() . A clear example of such a degeneracy is given 
in Fig. 1211 where the two limiting energy density profiles approximating the dyadosphere as introduced in Fig. 1171 are 
given for three different EMBH configurations corresponding to the same value of Edya = 3.1 x 10^* erg. The three 
configurations correspond respectively to the three different pairs (/.t,^): (10,0.76), (10^,0.27), (10*^,0.10). 

The corresponding dynamical evolution of the PEM pulse introduced in section IVIII and Ruffini et al. |l42l | is 
clearly different in the three cases. It is remarkable that when the collision with the remnant of the progenitor star is 



considered all these differences disappear. As usual (see section lVIIip we describe the baryonic content of the remnant 
by the parameter B. The PEMB pulse generated after the collision with the baryonic matter depends uniquely on 
the two parameters Edya and B. In Fig.|221the temperature in the laboratory frame is given for the PEM pulse and 
the PEMB pulse corresponding to the three configurations of Figl21l and B = A x 10~^. It is clear that while for 
the PEM pulse era the three configurations are markedly different, they do converge to a common behaviour in the 
PEMB pulse era. 

If we turn now to the effect of the distance between the EMBH and the baryonic remnant, we see that this 
degeneracy is further extended: while the three PEM pulse eras are quite different, the common PEMB pulse era is 
largely insensitive to the location of the baryonic remnant, see Fig. EHl We have plotted the three gamma factors 
in the PEM pulse era corresponding to the different configurations of Fig. 1211 and B = 10~^, in the two cases the 
baryonic remnant is positioned at different distances from the EMBH. 

If the PEM pulse has reached extreme relativistic regimes, the common value jcoii to which the three gamma factors 
drop in the collision with the baryonic matter of the remnant can be simply expressed by the large gamma limit of 
Eq.(|H3 

B + 1 

while the internal energy EcoU developed in that collision is simply given by the corresponding limit of Eq. (|85(l 

= -B + + 2B . (93) 

Edya 



This approximation applies when the final gamma factor at the end of the PEM pulse era is larger than jcoU, upper 
panel in Fig. |2S1 
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Figure 21: Three different dyadospheres corresponding to the same value of Edya = 3.1 x 10 erg and with different values of 
the two parameters fj, and ^ are given. The three different configurations are markedly different in their spatial extent as well 
as in their energy-density distribution. 
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Figure 22: The temperature of the plasma during the PEM pulse and PEMB pulse eras, measured in the laboratory frame, 
corresponding to the three configurations presented in Fig. I2il is given as a function of the laboratory time. The three different 
curves converge to a common one in the PEMB pulse era, which is therefore only a function of the Edya and B. The difference 
among the three curves in the early part of the PEMB pulse follows from having located the baryonic matter at a distance of 
50(r-ds — r+), which is different in the three cases. Such difference become negligible at large distances in the later phases of 
the evolution. 
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Figure 23: The gamma factors for the three configurations considered in Fig. l21l are given as a function of the radial coordinate 
in the laboratory frame. The two figures correspond to a baryonic remnant positioned respectively at Tin = 50(rds — r+) (left) 
and at ri„ = 5(rds — r+) (right). Again the convergence to a common behaviour, uniquely a function of E^ya and B for the 
late stages of the PEMB pulse, is manifest. 



Turning from these general considerations to the GRB data, this degeneracy in the PEMB pulse eras and their 
dependence on only two parameters E^ya and B has far reaching astrophysical implications for the identification of the 
source of GRBs. As we will see in the conclusions all the information obtainable from GRBs with a large value of the 
parameter B will lead to the determination of the above two parameters. An entire family of degenerate astrophysical 
solutions in the range of charges and masses given in Fig. 1171 are possible. The direct knowledge of the mass and 
charge of the EMBH can only be gained from the PEM pulse or from GRBs with very small values of i? — the so 
called "short bursts", see section IXIII and the conclusions. 



XI. THE APPROACH TO TRANSPARENCY: THE THERMODYNAMICAL QUANTITIES 



As above we shall consider as a typical case an EMBH of Edya = 3.1 x 10^^ erg and _B = 10 ^. The considerations 



As the condition of transparency expressed by Eq. (|86(l is reached the injector phase terminates. The electromagnetic 
energy of the PEMB pulse is released in the form of free-streaming photons — the proper GRB. The remaining energy 
of the PEMB pulse is released as an accelerated-baryonic-matter (ABM) pulse. 

We now proceed to the analysis of the approach to the transparency condition. It is then necessary to turn from 
the pure dynamical description of the PEMB pulse described in the previous sections to the relevant thermody- 
namic parameters. Also such a description at the time of transparency needs the knowledge of the thermodynamical 
parameters in all previous eras of the GRB 

ise an EMBH of Ldy 

will refer to a dyadosphere configuration described by the two limiting approximations shown in Fig. 1171 

One of the key thermodynamical parameters is represented by the temperature of the PEM and PEMB pulses. It is 
given as a function of the radius both in the comoving and in the laboratory frames in Fig. 1241 Before the collision the 
PEM pulse expands keeping its te mperature in the laboratory frame constant while its temperature in the comoving 
frame falls (see Ruffini et al. |142| 1 In fact Eas. (j66lfi7|l are equivalent to 

where the baryon mass-density is ps — and the thermal energy-density of photons and e+e -pairs is e = (TbT'*(1 + 
/e+e-)i is the Boltzmann constant and fe+e- is the Fermi-integral for e+ and e~ . This leads to 

e^^v = Eiy^, T^-f^V = const. (95) 

Since and in the PEM pulse are extremely relativistic, we have the equation of state p ~ e/3 and the thermal 
index F ~ 4/3 in the evolution of PEM pulse. Eq.(|n3) is thus equivalent to 

r^7V ~ const. (96) 

These two equations (I94|l and H96|l result in the constancy of the laboratory temperature T7 in the evolution of the 
PEM pulse. 

It is interesting to note that Eas. H95() and (|96|l hold as well in the cross-over region where T ~ rrieC^ and e+e~ 
annihilation takes place. In fact from the conservation of entropy it follows that asymptotically we have 

exactly for the same reasons and physics scenario discussed in the cosmological framework by Weinberg, see e.g. 
Eq. (15.6.37) of Weinberg (1972). The same considerations when repeated for the conservation of the total energy 
e'yV = e7^V following from Eq. I|94|l then lead to 

The ratio of these last two quantities gives asymptotically 

To = {T-f)T>m^c^ = {T'^)T<m^c^, (99) 

where To is the initial average temperature of the dyadosphere at rest. 

During the collision of the PEM pulse with the remnant we have an increase in the number density of e'^e~ pairs 
(see Fig. I24|) . This transition corresponds to an increase of the temperature in the comoving frame and a decrease of 
the temperature in the laboratory frame as a direct effect of the dropping of the gamma factor (see Fig. I2()|l . 

After the collision we have the further acceleration of the PEMB pulse (see Fig. 120(1 . The temperature now decreases 
both in the laboratory and the comoving frame (see Fig. I24|l . Before the collision the total energy of the e+e~ pairs 
and the photons is constant and equal to -Edya. After the collision 

Edya = E^Baryons + E^+^- + E^photons, (100) 

which includes both the total energy -Eg+e- -l~i?photons of the nonbaryonic components and the kinetic energy i^Baryons 
of the baryonic matter 

E^Baryons - PbV{^ - 1). (101) 

In Fig. 1251 we plot both the total energy i?e+e- -I- iJphotons of the nonbaryonic components and the kinetic energy 
-EBaryons of the baryonic matter as functi ons o f the radius for the typical case i?dya = 3.1 x 10''^'* erg and B = 10~^. 
Further details are given in Ruffini et al. |l43j . 
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Figure 24: Left) The temperature of the plasma in the comoving frame T'(MeV) (the sohd line) and in the laboratory frame 
-yT' (the dashed line) are plotted as functions of the radius in the unit of the dyadosphere radius rds- Right) The number 
densities n^+^- (T) (the solid line) computed by the Fermi integral and n^+^- (the dashed line) computed by the rate equation 
(see section lrv| are plotted as functions of the radius. T' <^ nieC? , two curves strongly divergent due to e'''e~-pairs frozen out 
of the thermal equilibrium. The peak at r ~ lOOrds is due to the internal energy developed in the collision. 

XII. THE P-GRBS AND THE "SHORT BURSTS". THE END OF THE INJECTOR PHASE. 

We now analyze the approach to the transparency condition given by Ea. (|86|l . For selected values of B we give the 
energy Ep-ghb of the P-GRB, and E'earyons of the ABM pulse. We clearly have 

Edya = Ep 

-GRB -^Baryons • 

(102) 

Taking into account the results shown in Figs. I24H25I we can repeat all the considerations for selected values of 
B. We shall examine values of B ranging from B = 10^^ only up to _B = 10^^: for larger values of B our constant 
slab approximation breaks down. We will sec in the following that this range does indeed cover the most relevant 
observational features of the GRBs. 

As clearly shown in Fig. 1201 both the final value of the gamma factor and the radial coordinate at which the 
transparency condition is reached depend very strongly on B. Therefore a strong dependence on B is also found in 
the relative values of Ep-grb and -EBaryons- 

We are now finally ready to give in Fig. El the crucial diagram representing the values of Ep-grb and i?Baryons 
in units of the E^ya as functions of B. This diagram, a universal one, is very important and is essential for the 
understanding of the GRB structure. 

We find that for small values of B (around 10^^) almost all the Edya is emitted in the P-GRB (see also our previous 
paper Ruffini et al. [l42]) and very little energy is left in the baryons. While for B ~ 10^^ roughly only 10^^ of the 
total initial energy of the dyadosphere is radiated away in the P-GRB and almost all energy is transferred to the 
baryons. 
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Figure 25: The energy of the non baryonic components of the PEMB pulse (the sohd hne) and the kinetic energy of the baryonic 
matter (the dashed line) in unit of the total energy are plotted as functions of the radius in the unit of the dyadosphere radius 



This behaviour is at the heart of the fundamental difference between the so called short bursts and long bursts. We 
have proposed in Ruffini et al. 145] that the short bursts must be identified with the P-GRBs in the case of very 
small B. There are a variety of reasons supporting this identification: 

1. For small values of B, i?Baryons is negligible, see Fig. |H| and consequently the intensity of the afterglow is also 
negligible and the entire energy E^ya is released into the P-GRB. This is clearly consistent with the absence of 
observed afterglows in the short bursts. 

2. The temperature of the P-GRB in the laboratory frame at the transparency point is a strongly decreasing 
function of i?, see Fig.|Sl 7T is rela ted to the energy corresponding to the peak of the photon- number spectrum, 
as described in Ruffini et al. 142]. This is also in very good agreement with the observed decrease of the 
hardness ratio between the short bursts and the long bursts (Kouveliotou et al. j84l]). 

3. The time rem , the duration of 90% of the energy emission as used in the current literature and discussed in 
Ruffini et al. jl4.ll | is plotted in Fig. 1211 for selected values of Edya and for different values of B. 

Before concluding a word of caution is needed about how to use the above results: all these considerations are based 
on the drastic approximations in the description of the dyadosphere presented in section IIIII see also Fig. |^ This 
treatment is very appropriate in estimating the general dependence of the energy of the P-GRB, the kinetic energy 
of the ABM pulse and consequently the intensity of the afterglow. Especially powerful is the establishment of the 
dependence of Ep-grb and i?Baryons on B (see Fig. EJ. As we will see in the next sections, this approximation is 
similarly powerful in determining the overall time structure of the GRB and especially the time of the release of the 
P-GRB with respect to the moment of gravitational collapse and the afterglow. 
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Figure 26: The duration computed with the Tgo criterion is represented as a function of the B parameter for three selected 
EMBH respectively with E^ya = 4.4 x 10^^ erg, Edya = 3.1 x 10^'* erg, Edya = 4.1 x 10^* erg going from the lower curve to the 
upper one. 

If, however, we turn to the detailed temporal structure of the P-GRB and its detailed spectral distribution, it is 
clear that the approximations given in section ITTTl is no longer valid. The detailed description of the formation of the 
dyadosphere as qualitatively expressed in Fig.0H|is now needed in all mathematical rigour with the full development 
of all its governing equations. Progre ss in this direction is being made at this moment (Cherubini et al. '28'|, RufRni 
& Vitagliano [l56, 157], Ruffini et al. |l59j 'l. This situation, however, provides a unique opportunity to follow in real 
time the general relativistic effects of the approach to the EMBH horizon as it occurs. In other words all direct general 
relativistic effects of the GRBs are encoded in the fine structure of the P-GRB. For the reasons given in section Ixl the 
information on the EMBH mass and charge can only come from the short bursts. 

This terminates the injector phase. We now turn to the Beam-Target phase in which the ABM pulse collides with 
the interstellar medium target and the afterglow is generated. We shall in the following sections review the basic 
theoretical treatment necessary for the description of these remaining eras and proceed then to the confrontation of 
the EMBH theory with the data. 

XIII. THE ERA IV: THE ULTRARELATIVISTIC AND RELATIVISTIC REGIMES IN THE 

AFTERGLOW 

In the introduction we have already expressed the basic assumptions which we have adopted for the description 
of the collision of the ABM pulse with the ISM. In analogy and by extension of the results obtained for the PEM 
and PEMB pulse cases, we also assume that the expansion of the ABM pulse through the ISM occurs keeping its 
width constant in the laboratory frame, although the results are quite insensitive to this assumption. We assume 
then that this interaction can be represented by a sequence of inelastic collisions of the expanding ABM pulse with a 
large number of thin and cold ISM spherical shells at rest with respect to the central EMBH. Each of these swept up 
shells of thickness Ar has a mass AMjsm and is assumed to be located between two radial distances ri and r2 (where 
r2 — ri = Ar <^ ri) in the laboratory frame. These collisions create an internal energy Ai^int. 

We indicate by Ae the increase in the proper internal energy density due to the collision with a single shell and 
by Pb the proper energy density of the swept up baryonic matter. This includes the baryonic matter composing the 
remnant around the central EMBH, already swept up in the PEMB pulse formation, and the baryonic matter from 
the ISM swept up by the ABM pulse: 

{Mb + MsnO 

PB = y • (103) 

Here V is the ABM pulse volume in the comoving frame, Mb is the mass of the baryonic remnant and Mism is the 



ISM mass swept up from the transparency point through the r in the laboratory frame: 

4:71 

Mism = TOpHism Y {r^ " To^) , (104) 

where the proton mass and riism the number density of the ISM in the laboratory frame. 

The energy conservation law in the laboratory frame at a generic step of the collision process is given by 

PB,7i'Vi + AM,,„c2 = (pB, ^ + ^MlEI^ + Ae) ^^V^, (105) 

where the quantities with the index "1" are calculated before the collision of the ABM pulse with an elementary shell 
of thickness Ar and the quantities with "2" after the collision, 7 is the gamma factor and V the volume of the ABM 
pulse in the laboratory frame so that V — 7V. 

The momentum conservation law in the laboratory frame is given by 

PB,liUr,Vi = (pB,y + ^^^^ + Ae") 72f/r.V2, (106) 



where Ur = •\/7^ — 1 is the radial covariant component of the four- velocity vector (see Ruffini et al. |l42l Il43l | and 
EqISSI). 

We thus obtain 

We can use for Ae the following expression 

^~~2 ~ V2 V2 ^ ' 

because we have assumed a "fully radiative regime" and so Sinti = 0. Substituting Ea. (|l()8|l in Ea. H107|l and applying 
Ea. (|109(l . we obtain: 



AEint = PBiVix 1 + 2^1 1- -— -psiVi lH 7^ ' (110) 



PBiVi \ PBiVi J ' V PBiV^] 



71 



72 = , ^ ■ (111) 



These relativistic hydrodynamic (RH) equations have to be numerically integrated. 

These are the actual set of equations we have integrated in the EMBH theory. In order to compare and contrast 
our results with the ones in the current literature, in section rXVIIII we have introduced the continuous limit of our 
equations and proceeded to have piecewise approximate power law solutions. We examine as well in section FXXI still 
under the above assumptions, the effects of a possible departure from homogeneity in the interstellar medium, still 
keeping the average density riism = const. Although these inhomogeneities are not relevant for the overall behaviour 
of the afterglow which we address here, they are indeed important for the actual observed flux and its temporal 
struc tures (see Rufhni et al. |l47j ). Also these considerations are affected by the angular spreading (Ruffini et al. 

[Hi). 



XIV. THE ERA V: THE APPROACH TO THE NONRELATIVISTIC REGIMES IN THE AFTERGLOW 



The only reason for addressing this last era is that the issue of the approach to nonrclativistic behaviour has been 
extensively discussed in the literature. In our treatment these results do not show any particular problems and the 



relativistic equations of the previous section continue to hold. In the specific example of GRB 991216 we will present 
in section FX Villi some analytic asymptotic expansions of these equations. 

This concludes the exposition of the different eras of the EMBH theory. It goes without saying that for the 
description of each era, all the preceding eras must necessarily be known in order to determine the space-time grid 
in the laboratory frame and its relation to the arrival times as seen by a distant observer. This is the basic message 
expressed in the RSTT paradigm. 

We can now turn to the comparison of the EMBH theory with the observational data. 

XV. THE BEST FIT OF THE EMBH THEORY TO THE GRB 991216: THE GLOBAL FEATURES OF 

THE SOLUTION 

For reasons already explained in the introduction, we use the GRB 991216 as a prototype. We will then later apply 
the EMBH theory to other GRBs. The relevant data of GRB 991216 are reproduced in Fig. 13 the data on the burst 
as recorded by BATSE Rapid Burst Response and t he data on the afterglow from the RXTE satellite (Corbet & 
Smith 31]) and the Chandra satellite (Piro et al., [l2dj ) , see also Halpern et al. f?^ . 

The data fitting procedure relies on three basic assumption: 

1. In the E-APE region, the source luminosity is mainly in the energy band 50-300 KeV, so we consider the flux 
observed by BATSE a good approximation of the total flux. 

2. In the decaying part of the afterglow, we assume that during the R-XTE and Chandra observations the source 
luminosity is mainly in the energy band 2-10 KeV, so we can again assume that the flux observed by these 
satellites is a good approximation of the total one. 

3. We have neglected in this paper the optical and radio emissions, since they are always negligible with respect 
to the X and 7 ray fluxes. In fact, even in the latest afterglow phases up to where the X-ray data are available, 
they are one order of magnitude smaller then the X-ray flux. 

These assumptions were initially adopted for the sake of simplicity, but have now also been justifled on the basis of 
the spectral description of the afterglow (Ruffini et al. 150]). 

As already emphasized in the previous sections, in the EMBH theory there are only two free parameters character- 
ising the afterglow: the energy of the dyadosphere, Edya, and the baryonic matter in the remnant of the progenitor 
star, parametrized by the dime nsionless parameter B. The location of the remnant has been assumed ~ 10^" cm. As 
discussed in Ruffini et al. |l44j and section Q the results are rather insensitive to t he ac tual density and location of 
the baryonic component but they are very sensitive to the value of B (Ruffini et al. jl43l ] ) . 

In Fig. Owe present the actual first results of fitting our EMBH theory to the data from the R-XTE and Chandra 
satellites, corresponding to selected values of Edya and B. There are three distinct features which are clearly evident 
as a function of the arrival time at the detector: an initial rising part in the afterglow luminosity which reaches a 
peak followed by a monotonically decreasing part. 

We have then proceeded to fine tune the two parameters in Fig. 1271 The main conclusions from our model are the 
following: 

1) The slope of the afterglow in the region where the experimental data are present is n = —1.6 and is in perfect 
agreement with the observational data. The index n in this region is rather insensitive to the values of the parameters 
Edya and B. The physical reason for this universality of the slope is rather remarkable since it depends on a variety 
of factors including the ultrarelativistic energy of the baryons in the ABM pulse, the assumption of constant average 
density in the ISM, the "fully radiative" conditions leading predominantly to X-ray emission, as well as all the different 
relativistic effects described in the RSTT paradigm (see also section IXVIIIjl . 

2) The afterglow fit does not depend directly on the parameters /z,^ but only through their combination Edya- 
Thus there is a 1-parameter family of values of the pair (/i, ^) allowed by a given viable value of Edya (see Figll7land 
section |X)| . 

3) By fine tuning the parameters of the best fit of the luminosity profile and time evolution of the afterglow the 
following parameters have been found: 

Edya = 4.83 X lO^^erg, B = 3x 10"^ . (112) 

After fixing in Eq. (|112(l the two free parameters of the EMBH theory, modulo the mass-charge relationship which 
fixes Edya, we can derive all the space-time parameters of the GRB 991216 (see Tab. IJ) as well as the explicit 
dependence of the gamma factor as a function of the radial coordinate (see Fig.UI). 




Detector arrival time (t^^) (s) 

Figure 27: Fine tuning of the best fit of tlie afterglow data of Chandra, RXTE as well as of the range of variability of the BATSE 
data on the major burst by a unique afterglow curve leading to the parameter values Edya = 4.83 x B = 3 x 10~^. 
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Figure 28: The distribution of the burst dura tions clearly shows two different classes of events: the "short bursts" and the 
"long bursts" (reproduced from Paciesas et al. |l07l |'l. 

Of special interest is the fundamental diagram of Fig. Its role is essential in interpreting all quantities measured 
in arrival time (the time of an observer in an incrtial frame at the detector) and their relations to the ones measured 
in the laboratory time by an observer in an inertial frame at the GRB source. The two times are clearly related by 
light signals (see Fig. I18|l and expressed by the integral Eq.(j33) and are also affected by the cosmological expansion 
(see section Ivji . 

XVI. THE EXPLANATION OF THE "LONG BURSTS" AND THE IDENTIFICATION OF THE 

PROPER GAMMA RAY BURST (P-GRB) 

Having determined the two free parameters of the EMBH theory, any other feature is a new prediction. An 
unexpected result soon became apparent, namely that the average luminosity of the main burst observed by BATSE 
can be fit by the afterglow curve (see Fig. ll()|l . This led us to the identification of the long bursts observed by BATSE 
with the extended afterglow peak emission (E-APE). The peak of this E-APE occurs at ~ 19.87 s and its intensity 
and time scale are in excellent agreement with the BATSE observations (see also Rufhni et al. |l47j '). It is clear 
that this E-APE is not a burst, but is seen as such by BATSE due to its high noise threshold (see also Ruffini et al. 
|l47l| 'l. Thus the ou tstan ding unsolved problem of explaining the long GRBs (see e.g. Piran 117], Salmonson et al. 
[162], Wilson et al. \l8l\ ) is radically resolved: the so called "long bursts" do not exist, they are just E-APEs (see 
Fig.EHl). 

We now turn to the most cogent question to be asked: where does one find the burst which is emitted when the 




Figure 29: A qualitative diagram showing the full picture of the model, with both P-GRB and E-APE. 



condition of transparency against Thomson scattering is reached? We have referred to this as the proper gamma ray 
burst (P-GRB) in order to distinguish it from the global GRB phenomena (see Bianco et al. [i^J, Ruffini et al. (144 1). 
We are guided in this search by two fundamental diagrams (see Fig. and Fig. 111(1 : 



1. In Ruffini et al. |14,'^ it is shown that for a fixed value of Edya the value of B uniquely determines the energy 
Ep-cRB of the P-GRB and the kinetic energy EBaryons of the ABM pulse which gives origin to the afterglow 
(see Fig. Illfl. For the particular values of the parameters given in Eq. H112|) . we find 

Ep-GRB = 7.54 X lO^^erg , Esaryons = 9.43 X lO^^erg (113) 

and then: 

Ep-GRB ^158^10-2^ (114) 

EBaryons 



2. One important additional piece of information comes from the differences in arrival time between the P-GRB 
and the peak of the E-APE, see Fig. ^2 Using the results of this figure and the numerical values given in Tab.lH 
we can retrace the P-GRB by reading off the time parameters of point 4 in Fig. |H1 Transparency is reached at 
21.57 s in comoving time at a radial coordinate r = 1.94 x 10^'' cm in the laboratory frame and at 8.41 x 10'^ s 
in arrival time at the detector. 

All this, namely the energy predicted in Eg. ((113(1 for the intensity of the burst and its time of arrival, leads to 
the unequivocal identification of the P-GRB with the apparently inconspicuous initial burst in the BATSE data. We 
have estimated from the BATSE data the ratio of the P-GRB to the E-APE over the noise threshold to be ^ 10~^, 
in excellent agreement with the result in Eq. 1(114(1 . see Fig. 1291 

It is important to emphasize that the diagrams in Fig. Eland Fig. II H are not universal, but depend on the dyadosphere 
energy. The corresponding diagrams for three selected Edya values {Edya — 5.29 x 10^"'^ erg, Edya = 4.83 x 10^^ erg 
and Edya — 4.49 x 10^^ erg) are given in Fig. l3Ub where we have plotted the energy of the P-GRB and of the E-APE 
as a function of B. The crossing of the intensity of P-GRB and E-APE occurs respectively at Bi — 6.0 x 10~^, 
B2 = 2.5 X 10^^ and B^ = 1.2 x 10^^ where Bi > B2 > B^. In Fig. 130b the same quantities are plotted as a function 
of the baryon mass Mb in units of solar masses and the opposite dependence occurs: Mi < M2 < M3. 

The physical reasons beyond these results is the following. We recall that the kinetic energy EBaryons and mass 
Mb of PEMB pulse are 

EBaryons 1)Mb Mb = BEdya (115) 

at the crossing point defined by 

EBaryons = Ep-GRB ~ ^rEdya- (116) 





Ej,s,.5.29*10=' erg 
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Figure 30: a) The same diagram of Fig. |^ is plotted for three different Edya values: Edya = 5.29 x 10^^ erg (dashed lines), 
Edya = 4.83 X 10^'^ erg (solid lines) and Edya ~ 4.49 x 10^^ erg (dotted lines), b) Same as in a) but plotted as a function of 
the baryonic mass Mb in units of solar masses instead of B. 



From these two equations, we obtain 



B=^^ z-—, (117) 

2(7o-l) 27o' ^ ^ 



7o is the Lorentz gamma factor of the PEMB pulse at the transparency point, where (see section IXlj) 

{npa'i.r + nB)(TT ^ UbCTt = 1, ns = - — 5-r — , (118) 

47rrjA7o 

At is the PEMB pulse thickness and ro the radial position at the transparency point. In addition, from the total 
energy conservation, we have 

(e + nB)7^47rr^A = const., (119) 
where e is the thermal energy of the PEMB pulse. In the regime ns 3> e, we have 

7o (120) 

and in the regime ^ e, we have 

70-^^0. (121) 

Considering the crossing point to occur in the second regime, we obtain at the crossing point 

B ^ {Edyayi, Mb - {Edyaf^. (122) 

These results are plotted in Figs. K-ilh -b. The agreement with the computed results is quite satisfactory. The differences 
can be attributed to the approximation adopted in Ea. l|121l) which is modified for high B values. 

The conclusion is that for increasing E^ya also the baryonic mass corresponding to the cross increases, but in 
percentage it increases less than Edya- 



XVII. CONSIDERATIONS ON THE P-GRB SPECTRUM AND THE HARDNESS OF THE SHORT 

BURSTS 



Regarding the P-GRB spectrum, the initial energy of the electron-positron pairs and photons in the dyadosphere 
for given values of the parameters can be easily computed following the work of Preparata et al. |l24j . We obtain 
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Figure 31: a) The B values corresponding to the crossings in Fig. l3Ub are plotted versus Edya (solid line). The function 
B oc -B^j^Y* obtained from a qualitative theoretical estimate (see Eg. 112211 ') is also plotted (dashed line), b) The AIb values 
corresponding to the crossings in Fig. I8()b are plotted versus Edya (solid line) . The function Ais oc -B^^^ obtained from a 
qualitative theoretical estimate (see Ea. p22^ 'l is also plotted (dashed line). 
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Figure 32: The temperature of the pulse in the laboratory frame for the first three eras of Fig. 1 of Ruffini et al. |l44ll is given 
as a function of the laboratory time. The numbers 1, 2, 3, 4 represent the beginning and end of each era. The two curv es refer 
to two extreme approximations adopted in the description of the dyadosphere. Details are given in Ruffini et al. |143|| and in 
section |3 



respectively T = 1.95 MeV and T — 29 A MeV in the two approximations we have used for the average energy density 
of the dyadosphere (see section EJ- It is then possible to follow in the laboratory frame the time evolution of the 
temperature of the electron-positron pairs and photons through the different eras, see Fig. The condition of 
transparency is reached at temperatures in the range of ~ 15 — 55 KeV at the detector, in agreement with the BATSE 
results. We emphasize that in the limit of B going to 10^^ in which the P-GRB coincides with the "short bursts" 
the spectrum of the P-GRB becomes harder in agree ment with the observational data (see Fig. and Band et al. 
0, Dermer et al. [13, Frontera et al. Norris et al. [T05l |). 

All the above are average values derived from the two approximations used in Fig. 1171 If one wishes to compare the 
EMBH theoretical results with the fine temporal details of the observational data on the P-GRB, a departure from 
this average approach will be needed and the fully time varying relativistic analysis outlined in Fig. 021 applies as will 
be further discussed in section IXXVII 



XVIII. APPROXIMATIONS AND POWER LAWS IN THE DESCRIPTION OF THE AFTERGLOW 



In addition to the BATSE data, there is also clearly perfect agreement with the decaying part of the afterglow data 
from the RXTE and Chandra satellites. 

We can also establish at this point a first set of conclusions on the luminosity power law index "n" which is a function 
depending strongly on the transformation i — > ta ^ if (see Fig. O. In the current literature such transformations 
and the corresponding n values are incorrect. Our theoretical value ntheo ~ —1.6 obtained for spherical symmetry for 
fully radiative conditions and constant density of the ISM is in agreement with observed riobs = — 1.616 ± 0.067. No 
evidence of beaming is found in GRB 991216. We shall return to this point in the conclusions. 

An extremely large number of papers in the literature deal with the power law index in the afterglow era. This 
issue has been particularly debated in connection with the aim of decreasing the energy requirements of GRBs by the 
effect of beaming (see e.g. Davies et al. 34], Mao & Yi It is currently very popular to infer the existence of 

beaming from the direct observations of breakings in the power-law index of the afterglow (see e.g. Derm er & Chiang 
[sgf , Gou et al. FH , Halper n et al. l73ll , Meszaros fc Re es |95l | , Meszaros et al. , Panaitescu et al. |ll2| , Panaitescu 
& Meszaros |11,'t| . R,hoads jl29t ll8l[ 7Sari. at al. |l67j ). Our aim here is to underline an often neglected point that 
the power law index of the afterglow is the result of a variety of factors including the very different regimes in the 
relation between the laboratory time t and the detector arrival time if presented in Fig. O No meaningful statements 
on the values of the power-law index of the afterglow can be made neglecting these necessary considerations expressed 
in the RSTT paradigm. This becomes particularly transparent from the power law expansion in the semianalytic 
treatments we present below. It is therefore not so surprising, as we will show in the next session, that the results 
obtained in the EMBH theory differ from the ones in the current literature. 



A. The approximate expression of the hydrodynamic equations 



We proceed to a first approximation and expand Eas. 11101 llll|l to second order in the quantity 



We obtain the following expressions: 



<1. (123) 



AiJint = (71 - 1) AMi,„,c2 - i ^ (AMism)' , (124) 



Mb + Mism 2 (Ms Msm) 

where we set A7 = 72 — 71 and have used the fact that pbiVi = {Mb + Mism)c'^- In the limit AEint — > dEint, 
A7 d'-f, and AMism dMism, neglecting also second order terms, where 

dMism = inr'^mpnismdr = inr'^mpnismvdt, v ^ (126) 

and where the ISM number density nism is assumed for simplicity to be nism = 1 cm~'^, we obtain: 

dEint = (7 - 1) dMismC^ , (127) 



Mb + JWism 

Eqs. 1)1271 | 128ll are limiting cases of Taub's hydrodynamical equations (Boccaletti et al. '25|, Landau & Lifshitz 
[83, Taub (177j '). They have been at times referred into the GRB literature as the Blandford-McKee equations (see 
Blandford & McKee [13 )■ It is clear that the application of these equations holds if Ea. H123|l applies. The behaviour 
of ^^'^^'"f as a function of the radius when Mism <ti Mb is: 

^ '-^ (129) 
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Figure 33: The factor ^p^^'y'^ is represented as a function of the radial coordinate. It is manifestly an increasing function. 

The condition Mism ^ Mb holds for GRB 991216 during the entire evolution of the system and so Eg. 1)123(1 is valid 
Eas. (|127ll28|) can be simply solved analytically (see e.g. Blandford & McKee [13). We then have: 



^ (Mb + M,^)^ + C 



(130) 



where 



! 7° - 1 
7o + l' 



(131) 



where we recall that To and 70 are the radial coordinate and the gamma factor at the transparency point and Mb is 
the initial baryonic mass of the ABM pulse. 

Eq. H13Q|1 is a differential equation for r {t) , namely 



dr 
cdt 



(Mb + Afism)' + C 



(Afs + A/ism)2 - Cj 

which can be integrated analytically with solution (see e.g. Abramowitz & Stegun Q) 



(132) 
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where m° = ^TTmpnismfo 1 B ~ ^ ) ^"^^ we recall that to is the laboratory time at the transparency point. 

Clearly the fulfilment of Eq. p23() has to be checked to ensure the validity of this solution. 
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Figure 34: a) The GRB flux emitted in laboratory time, b) the flux emitted in the arrival time, measured by an observer at 
rest with respect to the detector (see section IXVIIIfl . 

B. The approximate expression of the emitted flux 
From Eas. (|127ll28(l . it follows that the emitted flux in the laboratory frame is given by (see Fig. \'A4h) 



dE 
'dt 



Anr'^nismmpVj (7 - 1) c^, 



(134) 



and the corresponding flux in detector arrival time (see Fig. 134b ) by 



dE 



dt dE 
dtf,~dt 



t=t(ti) 



AimisairripC 



(135) 



t=t(t^) 



For the solution of these equations we distinguish four different phases (A-D). The first two correspond to era V. 



Phase A 



Just after the transparency condition is reached, the ISM matter involved is so small that we can approximately 
neglect the Mism term in Ea. (|130() and we have: 



7 ^ lo. 



(136) 



In the specific case of GRB 991216 we have 70 = 310.1, To = 1.94 x lO" cm, to = 6.48 x 10^ s, ta, ^ 4.21 x 10"^ s and 
tf^ ~ 8.41 X 10~^ s, where the index "o" refers to the quantities at the transparency point. We can then establish 
the following equation describing the ABM pulse motion in this phase: r [t) = vt with w ~ c. We can than use the 
following relation between laboratory time and arrival time: 



t = 2^oh, = p^tt (137) 



which is in perfect agreement with the full numerical computation (see Fig.lHl). 
We can substitute these equations into Eas. (|134ll35(l . obtaining: 



(138) 



in laboratory time and 



(X (139) 



dti (1 + z) 

in arrival time, assuming 7(7—1) ~ 7^. The results of the numerical integration of Eas. pU7ll08|l are in perfect 
agreement with these approximations (see Fig. I34() . 

Points P — the two maxima of the energy flux 

Since the contribution of the ISM mass in Eas. H13Qlll31|l can no longer be neglected, the value of 7 starts to 
significantly decrease (see Fig.lHJ and the flux reaches a maximum value. We integrate Ea. H134|l and Ea. H135|l using 
Ea. (|130|l for 7, assuming r (t) = vt with u ~ c and Ea. H137|l for the relation between the laboratory time and the 
arrival time (see Figs. I^5H ^. We can now obtain the point where the emitted flux reaches its maximum. In general, 
the location of the maximum of the flux, point P in Ruffini et al. |144| . will occur at different events, if considered in 
the arrival time (Pa) or in the laboratory time (Pl)- In this second case, the point P^ is determined by equating to 
zero the first derivative of Ea. H134|l . and we have: 



2 Mb 

«> -"'-'ism 

which in the case of GRB 991216 gives jp^ = 206.7 and 



270, (140) 

Pl 



~ 620.2. The maximum of the observed flux is 

Pl 

determined by equating to zero the first derivative of Ea. H135|l . We obtain: 



5 Mb 
IPa - 7t7o, 



Mi 

Mb 



isin 



^ 570, (141) 

Pa 



1550.5. 

Pa 



which in the case of GRB 991216 gives 7p^ ~ 258.4 and 

The results of the numerical integration of Eas. p07ll08|l are in perfect agreement with these approximations (see 
Fig. El- 



Phase B — the "golden value" n — —1.6 

In this phase 7 can no longer be considered constant and strongly decreases (see Fig. (S)). Afism is increasing, but 
V is still almost constant, equal to c. As a consequence, we can still say that r (t) = vt with w = c, but the relation 
between laboratory time and arrival time given in Ea. p37|l is no longer valid, and also Eg. 141(1 is no longer applicable 
in this phase (see Fig. inj. We can instead write the following "effective" relation: 

t^{tir\ (142) 

which is a result of a best fit of the numerical data in this region. Expanding the squares in Ea. H130|l . neglecting Mj^^^ 
with respect to M^ but retaining the terms in Afism and assuming 70 3> 1 we obtain: 

- 7pA=7p.%, (143) 



where rp^ and tp^ are the values of r and t at point Pl- Substituting this result into Eas. ()134|l . we obtain the emitted 
flux in the laboratory frame, given by 

dE 

— cx -flt^m^^t-^ , (144) 
dt 




Laboratory time (t) (s) 



Figure 35: The exact numerical solution for r (t) (solid line), together with the line r = ct (dotted line) and the fitting function 
given m 

Eq.dHHl (dashed line). 



and this is in good agreement with the full numerical computation (see Fig. I34|l . 

To obtain an analytic formula for the observed flux on the detector, we can still try to use the approximate relation 
between t and given by Eq.(@TJ: 



where 7 (t) is given by Ea. H143|l . We obtain: 



t = 2j{t)\^^^tt (145) 



■ (146) 



Using this formula in Ea. p85(l . we finally obtain: 



dE Iptp n 
(X — — — 

dti (1 + zY 



{tty (147) 



where we again assumed 7 (7 — 1) ~ 7^^. This results are not in agreement with the observational data, because the 
power-law index for the observed flux is —10/7 ~ —1.43, instead of the observed value —1.6. 

This is a confirmation that Ea. H145|l cannot be applied in this phase, as instead has been done by many authors in 
the current literature. We instead have to use Ea. H142|l . In fact, doing so we obtain the correct value: 



§^n,.n.{ti)-'\ (148) 



The results of the numerical integration of Eas. (|107llU8|l are in perfect agreement with these approximations (see 
Fig. l34|l . which implies that the approximate Ea. H127ll28() can still be used in this regime, but not Ea. (|41|l . which has 
to be replaced by an "effective" local power-law behaviour (see Ea. H142|l ). 



Phase C 



This new phase begins when 7 has decreased so much that the approximation r = is no longer valid (see Fig. l35|l . 
In the case of GRB 991216 this happens when 7 ~ 3.0, t ~ 1.5 x lO'^ s, ~ 2.9 x 10^ s and r ~ 4.4 x 10^"^ cm. In 
this entire phase, r {t) manifests the following behaviour typical of damped motion: 



•(i) = f (1 - e 



(149) 



where f, t* and r are constants that can be determined by the best fit of the numerical solution. In the present case 
of GRB 991216 we obtain: 



f ~ 1.101 X lO^^CTO, T ~ 2.072 X lO^s, ~ 4.52 x lOS. (150) 

It is important to note that this interesting behaviour, typical of a damped motion, does not lead to any power-law 
relationship for the emitted flux as a function of the laboratory time (see Fig. I34II . However, if we look at the observed 
flux as a function of the detector arrival time, we see that a power-law relationship still can be established, fltting the 
numerical solution. The result is: 

dE , J. -1.36 

^ « {<) • (151) 

This quite unexpected result can be explained because the relation between t and depends on r (t) in a nonpower- 
law behaviour. This fact balances the complex behaviour of the emitted flux as a function of the laboratory time, 
leading flnally again to a power-law behaviour arrival time. 

In this last phase, however, the flux decreases markedly, and from the point of view of the GRB observations, the 
most relevant regions are phases A and B described above, as well as the peak separating them. 



Phase D 



This last phase starts when the system approaches a Newtonian regime. In the case of GRB 991216 this occurs 
when 7 ~ 1.05, t ~ 5.0 x lO'^ s, if ~ 3.1 x 10^ s and r ~ 1.0 x 10^^ cm. In this phase r (t) is again approaching a 
linear behaviour, due to the velocity decreasing less steeply than in Phase C. The emitted flux as a function of the 
laboratory time still does not show a power-law behaviour, while the observed flux as a function of detector arrival 
time does, with an index n — —1.45 (see Fig. I34|l . 



XIX. THE POWER-LAW INDEX OF THE AFTERGLOW AND INFERENCES ON BEAMING IN GRBS 

The results obtained in the previous sections have emphasized the relevance of the proper application of the RSTT 
paradigm to the determination of the power-law index of the afterglow. Particularly interesting is the subtle interplay 
between the different regimes in the relation between the laboratory time and the arrival time at the detector clearly 
expressed by Fig. Inland the corres pond ing different regimes encountered in the flrst order expansion of the relativistic 
hydrodynamic equations of Taub |l77l | (see section IXVIII|I . It is interesting to compare and contrast our treatment 
with selected results of the current literature, in order to illustrate some relevant points (see Tab. lIII|l . We will consider 
the results in the literature only with reference to the limiting case which we address in our work: the condition of 
fully radiative emission. 

The first line of Tab. IIIII describes the ultrarelativistic regime, corresponding to an increasing energy flux of the 
afterglow as a function of the arrival time (phase A in previous section) . Our treatment and the results in the literature 
by Dermer et al. (see e.g. Bottcher & Dermer C hiang fc D ermer |2^ . De rmer et al. ^J) coincide. They agree 
as well with the results by Piran et al. (see e.g. Piran |ll6llll7| . Sari & Piran '168]). 

The second line corresponds to the relativistic regime, in which the energy flux of the afterglow, after having reached 
the maximum (point P in previous section), monotonically decreases (phase B in previous section). The dependence 
we have found of the gamma factor on the radial coordinate of the expanding ABM pulse does coincide with the 
one given by Dermer et al. and Piran et al. Our power law index n in this regime, which perfectly flts the data, 
however, is markedly different from the others. Particularly interesting is the difference between our results and those 
of Dermer et al: the two treatments coincide up to the last relation between the laboratory time and the arrival time 
at the detector. As explained in Eas. H147ll48|) . the two treatments differ in the approximation adopted in relating 
the laboratory time to the arrival time at the detector, illustrated in Fig. |2| Dermer et al. incorrectly adopted 
the approximation represented by the lower curve in Fig. and consequently they do not flnd agreement with the 
observational data. We have not been able to retrace in the treatment by Piran et al. the steps which have led to 
their different results. Special mention must be made of a result stated by Halpern et al. ^7^, the last entry in line 2, 



Table III: We compare and contrast the results on the power-law index n of the afterglow in the EMBH theory with other 
treatments in the current literature, in the limit of high energy and fully radiative conditions. The differences between the 
values of —10/7 ~ —1.43 (Dermer) and the results —1.6 in the EMBH theory can be retraced to the use of the two different 
approximation in the arrival time versus the laboratory time given in Fig. El See details in section IX Villi 
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Dermer et al. 
Bottcher & Dermer [20] 


Piran fll6j 
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Vietri [IgQj 


Halpern et al. 


Ultra-relativistic 
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Relativistic 


7 ~ r 
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n = -1.6 
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n = = -1.375 




n > -1.47 


Non-relativistic 


n = -1.36 
1.05 < 7 < 3.0 






n = -1.7 




Newtonian 


n = -1.45 
1< 7 < 1.05 











that an absolute lower limit for the power-law index n — 1.47 can be established on theoretical grounds. Such a result, 
clearly not correct also on the basis of our analysis, has been erroneously used ti support the existence of beaming in 
GRBs, as we will see below. 

The third line in Tab. IIIll is also in teres ting, treating the nonrelativistic limit (Phase C in previous section). This 
regime has been analysed by Vietri |l8(Tl |. avoiding the exact integration of the equations and relying on simple 
qualitative arguments. These results are not confirmed by the integration of the equations we have performed. This 
is an interesting case to be examined for its pedagogical consequences. Having totally neglected the relation between 
the laboratory time and the time of arrival at the detector, which we have illustrated in Fig.|51 and identifying = t, 
Vietri reaches a very different power law from our. Moreover, his solution brings to an underestimation of the radial 
coordinate: he estimated a radial coordinate of 1.1 x 10^^ cm at — 3.5 x 10'* s, while the exact computation shows 
a result greater than 3.0 x 10*^ cm (see Tab.lJ). On the other hand if one assumes, from the above mentioned identity 
t'^ = t, t = 3.5 X 10^ s, one obtains a gamma factor of ^ 300 (see Tab.lJ) in total disagreement with the nonrelativistic 
approximation adopted by Vietri. Quite apart from this pedagogical value, this nonrelativistic phase is of little interest 
from the observational point of view, due to the smallness of the flux emitted. 

For completeness, we have also shown our estimates of the index n as the Newtonian phase approaches in the last 
line of Tab. UTIl 

The perfect agreement between our theoretically predicted value for the power-law index, ntheo, and the observed 
one. Hobs, 

ntheo - -1.6, nabs = -1.616 ± 0.067, (152) 
confirms the validity of our major assumptions: 

1. The fully radiative regime. 

2. The constant average density of the ISM {riism = 1 proton /cm^). 

3. The spherical symmetry of the emission and the absence of beaming in GRB 991216. 



After the work of Mao & Yi [91| pointing to the possibility of introducing beaming to reduce the energetics of GRBs 
and after the discovery of the afterglow, many articles have appeared trying to obtain theoretical and observational 
evi dence for beamed emission in GRBs. The observations have ranged from radio (see e.g. Frail etal. R-ol et 
al. 1134^ 1 to optical (see e.g. Garnavich, et al 0|, Halpern et al. [tJI, Sagar et al. |l6l| . Schaefer |l69f ') all the way 
to X-rays. Particular attention has been devoted to relating the existence of beaming to possible breaks in the light 
curve slope, generally expected at a value of the gamma factor 

7=^, (153) 

where i^o is the beam opening angle. There are many arti cles on this s ubject: to men tion only th e mo st popular 
ones, we recall Meszaros et al. |97J, Panaitescu & Meszaros 'll5], Rhoads '129'. ' 13(1 Il3lj . Sari, at al. |l67j . Far from 
having reached a standard formulation, these approaches differ f rom each other in the expected time at which the 
break should take place up to a factor of 20 (see e.g. S ari, at al. |167| ). They differ as well for the opening angle of 
the beam, up to a factor of 3 (see e.g. Sari , at al. |l67j ). Disagreement still ex ists o n the number of breaking points: 
two i n the cas e of Panaitescu & Meszaros |ll5j . one in the case of Sari, at al. |l67j . one again in the case of Rhoads 
but differing in position from the one of Sari, at al. ,1673 . It has also been noticed that other authors 
have shown through numerical simulations that such a transition, if visible at all, is not very sharp (see e.g. Halpern 
et al. 73] ). 

Ample observational data have been obtained for the GRB 991216, in addition to the X-ray band, also in the optical 
and radio. For the reason mentioned at the beginning of section [XVI we only address in this article the problem of 
the 7- and the X-ray emission. In that respect, the main article addressing the issue of beaming in the X-rays for 
GRB 991216 is the one of Halpern et al. j23|- The key argument is based on the theoretical inequality claimed to 
exist for the power-law index n > —1.47 (see above). The fact that the observed X-ray decay rate is found to be 
n-obs = 1.616 ± 0.067 is interpreted by the authors as evidence for beaming. Moreover, the fact that the decay rate 
n = —1.6 has been observed before a steepening in the optical decay occurred at approximately 1 day of arrival time 
authorized an even more extreme proposal of a narrower beam in the X-rays within the optical beam. 

It is clear from the entire treatment which we have presented and the results of the EMBH theory given by 
nth en = —1.6 that there is no evidence for such a beaming, as already stated above. The motivation by Halpern et 
al. |2d| stems from the incorrect theoretical assumption of the existence of a lower limit in the afterglow power-law 
index n > —1.47. From our theoretical analysis the existence of n = —1.6 is clear proof of isotropic emission in the 
GRB 991216 and a clear test of the complete relativistic treatment of the source. The fact that the break in the index 
should be "achromatic" and the absence of beaming in the X-rays imply an absence of beaming also in the optical and 
radio bands. The observed steepening in the optical decay has to find an alternative explanation. Although this is 
not the subject of our pres ent work for the above mentioned reasons, we have found interesting the considerations by 
Panaitescu & Kumar |lllj . which find that "there are some major difficulties to apply a jet model to GRB 991216". 
They also state, still for GRB 991216, that "the steepening of the optical decay of a few days is not due to a jet effect, 
as suggested by Halpern et al. [l^ j but to the passage of a spectral break" . 

Concerning our own position on the possibility of beaming in GRBs, we would like just to remark that, from 
a preliminary analysis of beamed emission within the EMBH model, we have found some new features which are 
not encompassed by the results in the current literature, and they cou ld b ecome a distinctive signature for the 



discrimination of the existence or nonexistence of beaming (RufRni et al. |151|). T he study of the steepening in the 
optical and radio decay is addressed within the EMBH theory in RufRni et al. |150| |. 



XX. SUBSTRUCTURES IN THE E-APE DUE TO INHOMOGENEITIES IN THE INTERSTELLAR 

MEDIUM 

The afterglow is emitted as the ABM pulse plows through the interstellar matter engulfing new baryonic material. 
In our previous articles we were interested in explaining the overall energetics of the GRB phenomena and in this 
sense, we have adopted the very simplified assumption that the interstellar medium is a constant density medium 
with riisrn = l/crn^. Consequently, the afterglow emission obtained is very smooth in time. We are now interested in 
seeing if in this framework we can also explain most of the time variability observed by BATSE (see e.g. Fishman & 
Meegan H^), all of which except for the P-GRB should correspond to the beam-target phase in the IBS paradigm. 

We pursue this treatment neglecting the angular spreading due to off-axis scattering in the radiation of the afterglow, 
which will be presented in sections IXXIIIXXIIII 

Our goal is to focus in this simplified model on the basic energetic parameters as well as on the drastic consequences 
of the space-time variables expressed in the RSTT paradigm. 
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Figure 36: a) Flux of GRB 991216 observed by BATSE. The enlargement clearly shows the P-GRB (see Ruffini et al. |l45l|'l. 
b) Flux computed in the collision of the ABM pulse with an ISM cloud with the density profile given in Fig. m\ The dashed 
line indicates the emission from an uniform ISM with n = lcni~^. The dotted line indicates the BATSE noise level. 
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Figure 37: The density contrast of the ISM cloud profile introduced in order to fit the observation of the burst of GRB991216. 
The dashed line indicates the average uniform density n = lcm~^. 



Having obtained the two results presented in Fig. |H1 and Fig. |^ we can proceed to attack the specific problem of 
the time variability observed by BATSE. 

The fundamental point is that in both regimes the flux observed in the arrival time is proportional to the interstellar 
matter density: any inhomogeneity in the interstellar medium Auism/nism will lead correspondingly to a proportional 
variation in the intensity A/// of the afterglow. This result has been erroneously interpreted in the current literature 
as a burst originating in an unspecified "inner engine" . 

In particular, for the main burst observed by BATSE (see Fig. VSQh) we have 

(A///) = (An^sm/rMsm) - 5. (154) 

There are still a variety of physical circumstances which may lead to such density inhomogeneities. 

The additional crucial parameter in understanding the physical nature of such inhomogeneities is the time scale of 
the burst observed by BATSE. Such a burst lasts Ata — 20s and shows substructures on a time scale of ~ Is (see 
Fig. 136b ). In order to infer the nature of the structure emitting such a burst we must express these times scales in the 
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Figure 38: a) Same as Fig. 186b with the ISM cloud located at a distance of 3.17 x lO^^cm from the EMBH, the time scale of 
the burst now extends to ^ 1.58 x lO^s. b) Same as a) with the ISM cloud at a distance of 4.71 x lO^^cm from the EMBH, 
the time scale of the burst now extends to ~ 1.79 x lO^s. 



laboratory time (see RufRni et al. |l44j V Since we are at the peak of the GRB we have 7p^ ~ 258.5 (see Ea. H141|l ^ 
and Ata corresponds in the laboratory time to an interval 

Ai-l.OxlO^s, (155) 

which determines the characteristic size of the inhomogeneity creating the burst AL ^ 5.0 x lO^^cm (see Tab. 13 and 
Fig.EJ. 

It is immediately clear from Ea. (|154|l and Ea. H155|l that these are the typical dimensions and density contrasts 
corresponding to a small interstellar cloud. As an explicit example we have shown in Fig. 1371 the density contrasts 
and dimensions of an interstellar cloud with an average density < n >= \/cm?. Such a cloud is located at a distance 
of ~ 8.7 x lO^^cm from the EMBH, gives rise to a signal similar to the one observed by BATSE (see Fig. 136b'). 

It is now interesting to see the burst that would be emitted, if our present approximation would still apply, by 
the interaction of the ABM pulse with the same ISM cloud encountered at later times during the evolution of the 
afterglow. Fig. 138b shows the expected structure of the burst at a distance 4.1 x lO^^cm, corresponding to an arrival 
time delay of ~ 2 days, where the gamma factor is now 7* ~ 3.6. It is interesting that the overall intensity would 
be smaller, the intensity ratio of the burst relative to the average emission would remains consistent with Ea. H154|l . 

but the time scales of the burst would be longer by a factor 

{j^^ - 5 X 10^. Fig. EHb shows the corresponding 

quantities for the same ISM cloud located at a distance 6.4 x lO^^cm from the EMBH, corresponding to an arrival 
time delay of 1 month, where the gamma factor is ^ 1.5. 

We are going to analyze in the coming sections the modifications of this basic theory by the effect of the angular 
spreading: it will increase the accuracy of the fit obtained in Fig.|2niand will wash away all the features at late arrival 
time in the afterglow (see Fig.|2HI)- 



XXI. CONSIDERATIONS ON THE RELATIVISTIC BEAMING ANGLES AND ON THE ARRIVAL 

TIME 

We now generalize the results obtained in section [V] to consider also the effects due to the size of the emitting 
surface and of its curvature. The frequency lu and wave- vector k of photons emitted from the ABM pulse (see Fig. l39() 
expressed in the laboratory frame are: 

k = - (- sin i9u + cos i9v) , |k| = - , (156) 
c c 

where d is the angle (in the laboratory frame) between the radial expansion velocity and the line of sight, v is a unit 
vector along the radial expansion velocity of the ABM pulse, and u is a unit vector orthogonal to v oriented toward 



Figure 39: Qualitative description of the kinematics of the system. The big sphere is the expanding ABM pulse interacting 
with the ISM (not shown in the picture). The radius of the ABM pulse at time t is r (t). The generic point P on the ABM 
pulse, from which the photon is emitted, corresponds to a displacement angle from the line of sight. L is the distance of 
P from the observer. _Rt is the distance of the EMBH from the observer, r^s is the dyadosphere radius. Ro is defined by 
i?o = Rt — Tda- V is a unit vector along the radial expansion velocity, u is a unit vector orthogonal to v oriented toward rising 
i9. k is the momentum of the photons emitted toward the observer. Note that we have assumed i9 = A, i.e. k || _Rt (see text). 



rising We are assuming here that k and Rt are parallel, also for photons emitted with d ^ 0, so that X = 
This is clearly a good approximation, because the distance Rt corresponds to a redshift z ^ 1, while the radius of 
the emitting region is less than a light year in order of magnitude. Then the Lorentz boost along v to the comoving 
frame of the ABM pulse yields the corresponding comoving quantities: 

t^o = ( 1 - - cosz9) , Wo = |ko|c, (157) 



ko = - |k| sindu + 7 |k| [cosi^ - - j v, (158) 

In the comoving frame photons radiating out of the ABM pulse must have (see Ea. H158|l ): 

cos7?>-, (159) 
c 

because the component of the photon momentum in the comoving frame along the radial expansion velocity direction 
must be positive in order to escape. There will then be a maximum allowed 'd value 'dmax defined by cos "dmax = {v/c) 
(see Figs.lMEIll- 

Due to the high value of the Lorentz gamma factor (~ 300) for the bulk motion of the expanding ABM pulse, the 
spherical waves emitted from its external surface appear extremely distorted to a distant observer. Let us indicate by 
ta the arrival time at a detector of a photon emitted at a laboratory time t by the spherical surface of the relativistically 
expanding shell (see also sectional. Photons arriving at the same time ta will be emitted at different t as a function 
of the angle d (see Fig. I39II . The r elation between t and ta in the case of a constant 7 ~ 5 for expanding radio sources 
was found by Rees (see Rees |l27j |): 

ta = t (1 - ^cos??) . (160) 

For a constant expansion speed, the radius r (t) of the source is given by: 

r{t)=vt. (161) 
From Eas. H160lll61|l wc find the equation describing the "surface" emitting the photons detected at arrival time t^,: 
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Figure 40: Left) Not all values of i9 are allowed. Only photons emitted at an angle such that cosi9 > iv/c) can be viewed by 
the observer. Thus the maximum allowed i9 value Tl^max corresponds to cos'drnax ~ (v/c). In this figure we represent Tl^max (i.e. 
the angular amplitude of the visible area of the ABM pulse) in degrees as a function of the arrival time at the detector for the 
photons emitted along the line of sight (see text). In the earliest GRB phases w ~ c and so "dmax ~ 0. On the contrary, in 
the latest phases of the afterglow the ABM pulse velocity decreases and 'dmax tends to the maximum possible value, i.e. 90°. 
Right) The diameter of the visible area is represented as a function of the ABM pulse radius. In the earliest expansion phases 
(7 ~ 310) "dmax is very small (see left pane and Fig. I4H . so the visible area is just a small fraction of the total ABM pulse 
surface. On the other hand, in the final expansion phases "dmax —> 90° and almost all the ABM pulse surface becomes visible. 



which describes an ellipsoid of eccentricity ^ (see Rees [121 

In our case the ABM pulse Lorentz gamma factor is not constant (see Fig.|Hll, and so we must generalize Eas. (|160ll62|l 
to nonconstant expansion velocity. This can be done using the geometry of Fig. 1211 We set t = when the plasma 
starts to expand, so that r (0) = r^s, i-e. the dyadosphere radius. Let a photon be emitted at time t from the point 
P. Its distance from the observer is L. The time it takes to arrive at the detector is of course — . Thus its arrival 

c 

time, measured from the arrival of the first photon a time after its emission at i — 0, is: 

t, = t+--^, (163) 
c c 

where we have defined ia = when a photon emitted at t — and — reaches the observer. L is clearly given by: 



^Rl + r{ty -2RTr{t)cos^, (164) 

where at any given value of emission time t, cosi? can assume any value between ^-^^^ and 1 as noted above, where 

V (t) is the expansion speed of the ABM pulse at time t (see Eg. HISyH ). Now r (t) is less than one light year in order 
of magnitude while i?T corresponds to a redshift z ~ 1. Thus we can expand the right hand side of equation H164fl in 
powers of to first order: 

i ~ i?T (^1 - ^ cosz?^ , (165) 

which corresponds to assuming L to be equal to its projection on the line of sight (see Fig. 1^. Substituting H165fl 
into p63|l yields: 

ta = i+ — -— cosi^, (166) 
c c 

where we have used the fact that Rt = Rq + r^s (see Fig. I39|) . For r (t) we can use the following expression: 



r{t)= f V {t') dt' + n,, 
Jo 



(167) 
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Figure 41: Left) This figure shows the temporal evolution of visible area of the ABM pulse. The dashed half-circles are the 
expanding ABM pulse at radii corresponding to diflerent laboratory times. The black curve marks the boundary of the visible 
region. The EMBH is located at position (0,0) in this plot. Again, in the earliest GRB phases the visible region is squeezed 
along the line of sight, while in the final part of the afterglow phase almost all the emitted photons reach the observer. This 
time evolution of the visible area is crucial to the explanation of the GRB temporal structure. Right) Due to the extremely 
high and extremely varying Lorentz gamma factor, photons reaching the detector on the Earth at the same arrival time are 
actually emitted at very different times and positions. We represent here the surfaces of photon emission corresponding to 
selected values of the photon arrival time at the detector: the eqmtemporal surfaces (EQTS). Such surfaces differ from the 
ellipsoids described by Rees in the context of the expanding radio sources with typical Lorentz factor 7 ~ 4 and constant. In 
fact, in GRB 991216 the Lorentz gamma factor ranges from 310 to 1. The EQTSes represented here (solid lines) correspond 
respectively to values of the arrival time ranging from 5 s (the smallest surface on the left of the plot) to 60 s (the largest one 
on the right). Each surface differs from the previous one by 5 s. To each EQTS contributes emission processes occurring at 
different values of the Lorentz gamma factor. The dashed lines are the boundaries of the visible area of the ABM pulse and 
the EMBH is located at position (0, 0) in this plot. Note the different scale on the two axes, indicating the very high EQTS 
"effective eccentricity" . The time interval from 5 s to 60 s has been chosen to encompass the E-APE emission, ranging from 
7 = 308.8 to 7 = 56.84. 



so that equation IjlGGfl can be written in the form: 

Jo V (f) dt' + y-ds T-ds .,„„s 

ts. = t COS 17 H , (168) 

c c 

which reduces to Eq. IjlfiOII only if v is constant and r^s is negligible with respect to r (t) . 

Also from Eq. (|168|l we can obtain the equation describing the surface that emits the photons detected at an arrival 
time ta- In this case, we no longer have ellipsoids of constant eccentricity ^. Since the velocity is strongly varying 
from point to point, we have more complicated surfaces like the profiles reported in Fig. 1411 where at every point there 
will be a tangent ellipsoid of a given eccentricity, but such an ellipsoid varies in eccentricity from point to point. 

For a fixed time t of emission in Ea. l|168|l . the allowed angular interval | < cosi? < 1 leads to a corresponding 
smearing of the arrival time ta over the interval 

72c (1 + -) 

We need now to correct Ea. Hlf)8|l for the cosmological expansion effects to get the wanted relation between t and 
if. We recall that (see section IVTl 



if = (l + z)t„. 



(170) 



where z is the cosmological redshift. Our final relation is therefore: 



*a = (1 + 2) i - " ^ ^ cosi9 + . (171) 



XXII. THE EMISSION PROCESS TAKING OFF-AXIS CONTRIBUTIONS INTO ACCOUNT 

We now take into consideration the contributions of the off-axis emission to the afterglow to see if the previous 
positive results still hold and if some of the problems just stated can be overcome by a more detailed and relativistic 
treatment. The corresponding computation for the P-GRB structure will be presented elsewhere, where the time 
evolutions of the dyadosphere formation and its c onsequenc es on the P-G RB s tructures are presented following the 
work of Cherubini et al. 28], Ruffini & Vitaghano |l5filll57j . Ruffini et al. |l59f . The effects on the P-GRB structure 
of the dyadosphere formation dominate those due to the angular spreading. 

Following Eas. HllUHlll|l . we recall that in the comoving frame of the expanding ABM pulse we suppose that the 
internal energy due to kinetic collision is instantly radiated away and that the corresponding emission is isotropic. As 
in section^ let Ae be the internal energy density developed in the collision. In the comoving frame the energy per 
unit of volume and per solid angle is simply 

^) (172) 

due to the fact that the emission is isotropic in this frame. The total number of photons emitted is an invariant 
quantity independent of the frame used. Thus we can compute this quantity as seen by an observer in the comoving 
frame (which we denote with the subscript "o") and by an observer in the laboratory frame (which we denote with 
no subscripts). Doing this we find 

dN^ f / dN^ \ , . „ 

- I ' ^1 A^^cosi^, (173) 



dtdfldT, J shell \dtdfldY] 

where cosi^ comes from the projection of the elementary surface of the shell on the direction of propagation and 
A = 7(1 — /3cosi?) is the Doppler factor introduced in the two following differential transformation 

dflo ^dnx A-^ (174) 

for the solid angle transformation and 

dto ^dtx A^i (175) 

for the time transformation. The integration in dT, is performed over the visible area of the ABM pulse at laboratory 
time t, namely with < "0 < "dmax and 'dmax defined in section lXXIl fsee Ec ljlSyf) and Figs. l4(JH41|l . An extra A factor 
comes from the energy transformation: 

Eo^ExK. (176) 
See also Chiang & Dermer l^^. Thus finally we obtain: 



dE _ r / dE \ 
dtdndY, ~ J^^^ji \dtdndY. ) ^ ' 



A-'^cosi^. (177) 



Doing this we clearly identify ( dtdUdT. ) o '"^ energy density in comoving frame up to a factor ^ (see Ea. ()172|) '). 
Then we have: 

[ — i; cost^A-^ds, (178) 

dtdn J shell 4^ ' ^ ' 

where the integration in is performed over the ABM pulse visible area at laboratory time i, namely with Q <'d < 
'&max and dmax defined in section IXXII 

Ea. H178|l gives us the energy emitted toward the observer per unit solid angle and per unit laboratory time t in the 
laboratory frame. But what we really need is the energy emitted per unit solid angle and per unit detector arrival 



time t^, so we must use the complete relation between and t given in Ea. (|171() . First we have to multiply the 
integrand in Ea. ll78|l by the factor [dt/df^ to transform the energy density generated per unit of laboratory time 
t into the energy density generated per unit arrival time f^. Then we have to integrate with respect to over the 
equitemporal surface (EQTS, see section IXXI|l of constant arrival time instead of the ABM pulse visible area at 
laboratory time t. The analog of Ea. H178|) for the source luminosity in detector arrival time is then: 

dE^ f Ae „ , A dt ,^ , „„, 

-rzrh= — w cosi? A"-* — dS. (179 

dtidn J EQTS 47r dti 

It is important to note that, in the present case of GRB 991216, the Doppler factor A"** in Ea. ll79|l enhances the 
apparent luminosity of the burst, as compared to the intrinsic luminosity, by a factor which at the E-APE is in the 
range between 10^° and lO^^i 

To perform the numerical integration of Eq. H179|l we have implemented the following procedure for each fixed value 
of the laboratory time t: 

1. We fix the laboratory time t. 

2. We divide the interval of the allowed values {y (t) /c) < cosd < 1 into N small steps, each one of amplitude 

A,icos^)^'-MlM, (180) 

3. We select n directions defined by: 

cosi9„ = l- nA (cosi9) , (181) 

where n is an integer, < rt < and so i9o = and — "dmax- 

4. For each we compute with Ea. H179| ) the contribution to the afterglow luminosity arising from an angular 
aperture corresponding to A^r (cost?) around such a direction. 

5. We compute for each value of n the corresponding values of the arrival time using Ea. (|171|l . 

To obtain the total luminosity at arrival time we sum together all the above contributions corresponding to the 
same f^^. 

We first apply this treatment to the analysis of the afterglow using assumptions 1 and 2 of section. ^1 namely that 
the ISM density is constant riism =< nism >— 1 particle/cm"^ and that the ABM is spherically symmetric. 

Fig. 021 compares the new result for the afterglow luminosity as a function of the detector arrival time with the 
previous one obtained in section FXVI bv neglecting off-axis emission. The main conclusions are: 

1. The total energy emitted both in the radial approximation and in the full computation with the off-axis emission 
is conserved. This is a necessary condition for checking the consistency of the model. 

2. The slope of the decreasing part of the afterglow is unchanged. We emphasize once more the great advantage 
of the radial approximation which has allowed to obtain an analytic expression for this slope. 

3. The final phase of the afterglow (7 < 2) is largely affected by the late arrival of the radiation emitted at large 
angles. In fact in the radial approximation the luminosity goes abruptly to zero when 7 reaches 1 while in the 
new complete treatment the behavior is much smoother due to the delayed arrival of the radiation emitted at 
large angles. Consequently, enforcing the energy conservation, in the rising part of the afterglow the luminosity 
in the new treatment is shown to be slightly smaller than in the radial case. 

In order to acquire a better understanding of the effects of angular spreading, we have found it helpful to analyze the 
radiation emitted from selected angles between and "dmax ■ This is in addition to the integration results presented 
in Fig. 221 In Fig- 021 we show the results of such an analysis plotting the contributions to the total luminosity 
corresponding to selected values of n in Ea. H181() . We easily see that radiation emitted at large angles is time shifted 
with respect to that emitted near the line of sight. In fact the afterglow peak occurs later going to higher n values 
(see Fig. 1221 ■ 
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Figure 42: Left) The predicted afterglow curve for GRB 991216 assuming a constant ISM density equal to 1 particle/cm'^ and 
taking into account all the effects due to off-axis emission (solid line). For comparison we plot also the corresponding curve 
obtained in the simple radial approximation (dashed line). We see that this last curve falls sharply to zero when the ABM 
pulse reaches 7 = 1, while the first one has a much smoother behavior due to the time delay in the arrival of the photons 
emitted at large i9. Recall that when 7 tends to 1, the maximum allowed values of i9 tend to 90°. Right) This figure shows 
how the radiation emitted from different angles contributes to the afterglow luminosity. The solid line on the top of the picture 
is the total luminosity as in the previous plots. The other dashed and dotted curves represent the radiation components 
corresponding to selected values of n in Ea. lfTSTji . From the upper to the lower one they corresponds respectively to n = 0, 
n = 0.05A'', n = 0.25-/V, n = 0.5-/V, n — N, where in this plot = 200. We can easily see that the radiation emitted at large 
angles (n = A'') is time shifted with respect to that emitted near the line of sight (n = 0). 



Table IV: For each ISM density peak represented in Fig. 1431 we give the initial radius r, the corresponding comoving time r, 
laboratory time t, arrival time at the detector t^, diameter of the ABM pulse visible area d„, Lorentz factor 7 and observed 
duration Atf of the afterglow luminosity peaks generated by each density peak. In the last column, the apparent motion in the 
radial coordinate, evaluated in the arrival time at the detector, leads to an enormous "superluminal" behavior, up to 9.5 x 10^ c. 



Peak 


r(cm) 


t{s) 


tis) 


ti{s) 


dy{cm) 


Atiis) 


7 


"Superluminal" 


A 


4.50 X 10'" 


4.88 X lO-" 


1.50 X 10" 


15.8 


2.95 X 10'* 


0.400 


303.8 


9.5 X lO^c 


B 


5.20 X 10'" 


5.74 X 10^ 


1.73 X lO*' 


19.0 


3.89 X 10"' 


0.622 


265.4 


9.1 X lO^c 


C 


5.70 X 10'*^ 


6.54 X 10^ 


1.90 X 10*^ 


22.9 


5.83 X 10"' 


1.13 


200.5 


8.3 X lO^c 


D 


6.20 X 10'*^ 


7.64 X 10^ 


2.07 X 10" 


30.1 


9.03 X lO'"* 


5.16 


139.9 


6.9 X lO^c 


E 


6.50 X 10'" 


9.22 X 10^ 


2.17 X lO'^ 


55.9 


2.27 X 10'^ 


10.2 


57.23 


3.9 X lO^c 


F 


6.80 X 10'" 


1.10 X 10* 


2.27 X 10" 


87.4 


2.42 X 10'^ 


10.6 


56.24 


2.6 X lO^c 



XXIII. THE E-APE TEMPORAL SUBSTRUCTURES TAKING INTO ACCOUNT THE OFF-AXIS 

EMISSION 

We are now ready to reconsider the problem of the ISM inhomogeneity generating the temporal substructures in 
the E-APE by integrating on the EQTS surfaces and improving on the considerations based on the purely radial 
approximation. We have created (see details in Ruffini et al. 153]) an ISM inhomogeneity "mask" (see Fig. 1431 and 
Tab. ^ with the main criteria that the density inhomogeneities and their spatial distribution still fulfill < riism >= 
1 particle/cm'^. 

The results are given in Fig. 1441 We obtain, in perfect agreement with the observations: 

1. the theoretically computed intensity of the A, B, C peaks as a function of the ISM inhomogneities; 

2. the fast rise and exponential decay shape for each peak; 



3. a continuous and smooth emission between the peaks. 
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Figure 43: The density profile ("mask") of an ISM cloud used to reproduce the GRB 991216 temporal structure. As before, 
the radial coordinate is measured from the black hole. In this cloud we have six "spikes" with overdensity separated by low 
density regions. Each spike has the same spatial extension of lO'^^ cm. The cloud average density is < riism >~ 1 particle/cm^. 
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Figure 44: Left) The BATSE data on the E-APE of GRB 991216 (source: BATSE GRB hght curves Isl) together with an 
enlargement of the P-GRB data (source: BATSE Rapid Burst Response 0). For convenience each E-APE peak has been labeled 
by a different uppercase Latin letter. Right) The source luminosity connected to the mask in Fig. 1431 is given as a function 
of the detector arrival time (solid "spiky" line) with the corresponding curve for the case of constant niem = 1 particle/cm"^ 
(dashed smooth line) and the BATSE noise level (dotted horizontal line) . The "noise" observed in the theoretical curves is 
due to the discretization process adopted, described in Rufiini et al. |l53j| . for the description of the angular spreading of the 
scattered radiation. For each fixed value of the laboratory time we have summed 500 different contributions from different 
angles. The integration of the equation of motion of this system is performed in 22, 314, 500 contributions to be considered. 
An increase in the number of steps and in the precision of the numerical computation would lead to a smoother curve. 



Interestingly, the signals from shells E and F, which have a density inhomogeneity comparable to A, are undetectable. 
The reason is due to a variety of relativistic effects and partly to the spreading in the arrival time, which for A, 
corresponding to 7 = 303.8 is 0.4s while for E (F) corresponding to 7 = 57.23 (56.24) is of 10.2 s (10.6 s) (see Tab. II VI 



In the case of D, the agreement with the arrival time is reached, but we do not obtain the double peaked structure. 
The ABM pulse visible area diameter at the moment of int eraction w ith the D shell is ~ 1.0 x 10^^ cm, equal to the 
extension of the ISM shell (see Tab. II VI and Ruffini et al. |l48i Il53j ). Under these conditions, the concentric shell 
approximation does not hold anymore: the disagreement with the observations simply makes manifest the need for a 
more detailed description of the three dimensional nature of the ISM cloud. 

The physical reasons for these results can be simply summarized: we can distinguish two different regimes corre- 
sponding in the afterglow of GRB 991216 respectively to 7 > 150 and to 7 < 150. For different sources this value may 
be slightly different. In the E-APE region (7 > 150) the GRB substructure intensities indeed correlate with the ISM 
inhomogeneities. In this limited region (see peaks A, B, C) the Lorentz gamma factor of the ABM pulse ranges from 
7 ~ 304 to 7 ~ 200. The boundary of the visible region is smaller than the thickness AR of the inhomogeneities (see 
Fig. 1411 and Tab. llVjl . Under this condition the adopted spherical approximation is not only mathematically simpler 
but also fully justified. The angular spreading is not strong enough to wipe out the signal from the inhomogeneity 
spike. 

As we descend in the afterglow (7 < 150), the Lorentz gamma factor decreases markedly and in the border line 
case of peak D 7 ~ 140. For the peaks E and F we have 7 50 and, under these circumstances, the boundary of 
the visible region becomes much larger than the thickness AR of the inhomogeneities (see Fig. 1411 and Tab. II V|) . A 
three dimensional description would be necessary, breaking the spherical symmetry and making the computation more 
difficult. However we do not need to perform this more complex analysis for peaks E and F: any three dimensional 
description would a fortiori augment the smoothing of the observed flux. The spherically symmetric descrip tion of 
the inhomogeneities is already enough to prove the overwhelming effect of the angular spreading (Ruffini et al. |l53j |). 

On this general issue of the possible explanation of the observed substructures with the ISM inhomogeneities, there 
exists in the literature two extreme points of view: the one by Fenimore and collab orators (s ee e.g. Fenin iore et al. 
[4^. Feni more [49l| . Fenimore et al. [H^l) and Piran and collaborators (see e.g. Piran |ll6llll7l |. Piro et al., |l20l |. Sari 

6 Piran Il65ll ') on one side and the one by Dermer and collaborators (Dermer 0|, Dermer et al. |43|, Dermer & 
Mitman [43|) on the other. 

Fenimore and collaborators have emphasized the relevance of a specific signature to be expected in the collision of 
a relativistic expanding shell with the ISM, what they call a fast rise and exponential decay (FRED) shape. This 
feature is confirmed by our analysis (see peaks A, B, C in Fig. I44() . However they also conclude, sharing the opinion 
by Piran and collaborators, that the variability observed in GRBs is inconsistent with causally connected variations 
in a single, symmetric, relativistic shell interacting with the ambient material ("external shocks") (Fenimore et al. 
[SO'I). In their opinion the solution of the short time variability has to be envisioned within the protracted activity of 
an unspecifi ed "i nner engine" (Sari & Piran 165j); see as well Meszaros & Rees ^98, 99J, Meszaros (lOOij , Panaitescu 
fc Meszaros [Tl^. Rees fc Meszaros 128l]. 

On the other hand, Dermer and collaborators, by considering an idealized process occurring at a fixed 7 — 300, 
have reached the opposite conclusions and they purport that GRB light curves are tomographic images of the density 
distributions of the medium surrounding the sources of GRBs (Dermer & Mitman 0|). 

From our analysis we can conclude that Dermer's conclusions are correct for 7 ^ 300 and do indeed hold for 7 > 150. 
However, as the gamma factor drops from 7 ^ 150 to 7 ~ 1 (see Fig (SI, the intensity due to the inhomogeneities 
markedly decreases also due to the angular spreading (events E and F). The initial Lorentz factor of the ABM pulse 

7 ~ 310 decreases very rapidly to 7 ~ 150 as soon as a fraction of a typical ISM cloud is engulfed (see Tab. II Vp . We 
conclude that the "tomography" is indeed effective, but uniquely in the first ISM region close to the source and for 
GRBs with 7 > 150. 

One of the most striking feature in our analysis is clearly represented by the fact that the inhomogeneities of a mask 
of radial dimension of the order of 10^^ cm give rise to arrival time signals of the order of 20 s. This outstanding result 
implies an apparent "superluminal velocity" of ~ lO^c (see Tab. lIVIl . The "superluminal velocity" here considered, 
first introduced in Ruffini ct al. 144], refers to the motion along the line of sight. This effect is proportional to 7^. It 
is much larger than the one usu ally c onsidered in the literature, within the context of radio sources and microquasars 
(see e.g. Mirabel & Rodriguez [lOll|'l. referring to the component of the velocity at right ang les t o the line of sight 
(see details in Ruffini et al. |153| V This second effect is in fact proportional to 7 (see Rees Il27|). We reca ll tha t 
this "superluminal velocty" was the starting point for the enunciation of the RSTT paradigm (Ruffini et al. |l44j ). 
emphasizing the need of the knowledge of the entire past worldlines of the source. This need has been further clarified 
here in the determination of the EQTS surfaces (see Fig. which indeed depend on an integral of the Lorentz gamma 
factor extended over the entire past worldlines of the source. In turn, therefore, the agreement between the observed 
structures and the theoretical predicted ones (see Figs. IHlfl^ is also an extremely stringent additional test on the 




values of the Lorentz gamma factor determined as a function of the radial coordinate within the EMBH theory (see 
Fig.lHl. 



XXIV. ON THE ISTANTANEOUS SPECTRUM OF GRBS 

Variability on the shortest time scale ever observed in nature is the main message we have acquired from the 
theoretical understanding of GRB astrophysical phenomena (see sections IVIVIII - IXljl . This situation is made even 
more extreme by the fact that astronomical and astrophysical observations are carried out in the "pathological" time 
coordinate of the photon arrival time at the detector (see section 0), whereby the first 10^ seconds of the GRB 
phenomena are further compressed in ~ 0.1 seconds (see Tab. ^) and further enhanced. The understanding that in 
these first 10^ seconds four different physical eras of the GRB phenomena occur has led us to a sentiment of natural 
skepticism toward any global or average description of the GRB phenomenon. We start to realize that such average 
descriptions mediate on totally different physical processes and lead to very questionable results. Such skepticism was 
even strengthened as soon as we realized that the characteristic quantities usually adopted for the description of the 
bursts, e.g. T50 and Tgo, which so many tried for years t o exp lain wit hin the context of the internal shock model (see 
e.g. Fenimore j4£|, Fenimore et al. 50J, Paczyiiski & Xu |l09j . Piran |116| . Rees & Meszaros [128j . Sari & Piran .IGSj l 
and references therein) were actually referring not at all to the bursts but to the extended emission from the peak of 
the afterglow: the E-APE! In this sense they were quite irrelevant for understanding the nature of the GRB source and 
were at most of interest for inquiring the structure of the ISM a few light months away from the source! It has been 
then with this sentiment of marked skepticism toward a global approach that we have started to consider the problem 
of the spectrum of GRBs and the validity of the band relation (Band et al. 0). To attempt an integral description 
of the spectra of the GRBs extending over 10^ seconds in arrival time is clearly meaningless. It mediates on tw o 
conceptually physically different phases of GRBs: the injector phase and the beam-target phase (RufFini et al. |145| ). 
In addition, in each of these phases many specific eras are present and each one of these eras needs due attention and 
can lead in principle to a different instantaneous spectrum. The fact that the spectral distribution observed by Band 
was a non-thermal one has been a very strong objection to consider any thermal spectrum. The situation became so 
extreme in the recent years that the sole appearance of a thermal spectrum in any part of a theoretical paper was 
considered a good reason for rejecting the paper by a refereed journal and to discard the validity of that work. 

Ha ving developed the very powerful theoretical tool of the EQTS surfaces (see section IXXII and Ruffini et al. 
jEHUssJ) and having been successful in having established the substructure of the E-APE, in addition t o the features 
of the afterglow, we have decided to approach the instantaneous spectra of the GRBs in Ruffini et al. |l5d| . In the 
abstract of that paper, we summarize as follows the results: "A theoretical attempt to identify the physical process at the 
basis of the afterglow emission of GRBs is presented, assuming GRB 991216 as a prototype. Such a physical process is 
identified in a mechanism leading to a thermal emission occurring in the comoving frame of the shock wave originating 
the GRBs. For the determination of the actually observed GRB luminosities and spectra at a given arrival time, the 
concept of equitemporal surfaces (EQTS, see Ruffini et al. 'IAS]) has to be implemented: the final results comprehend 
an integration over an infinite number of planckian spectra, weighted by appropriate relativistic transformations, each 
one corresponding to a different viewing angle in the past light cone of the observer. The relativistic transformations 
have been computed on the gr ound of the kno wledge of the already determined equations of motion of GRBs within the 
EMBH theory (Ruffini et al. \1AA . fiZl . \li §]). The only free parameter of the present theory is then the dimension of 
the "effective cavity" where the thermalization process occurs. A precise fit (x^ — 1.08) of the observed luminosity in 
the 2-10 keV band of GRB 991216 is presented as well as a detailed estimate of the observed luminosity in the 50-300 
keV band and of the expected one in the 10-50 keV band. The long awaited explanation of the observed hard-to-soft 
transition in GRBs is also presented" (Ruffini et al. ,150]). It is interesting that this theoretical result, which up 
to few years ago were hardly testable due to the paucity of photons collected by the detectors, have now become a 
necessity in order to interpret the splendid observational results o f the new families of space observatory like Chandra 
and XMM (see e.g. Borozdin & Trudolyubov Watson et al. [isl Il85l | 1 . 

Prior to our work, the possibility that the non-thermal looking spectrum of GRBs can be found as a superposition of 
a set of thermal blackbody spectrum was forcefully expressed in a simple paper by Blinnikov, Kozyreva & Panchenko 
[Tsl. These three authors have expressed in an analytic treatment that indeed the time integration of the black 
body planckian spectrum with a temperature varying with time following a simple power-law and expanding with 
another power-law can lead to a non-thermal spectrum in agreement with the observed Band relations. To obtain 
this result, they use two indexes for their qualitative analysis to be fitted by the observational data. Toward the 
end o f the ir paper they finally quoted "In reality, not o nly t ime, but also space integration takes place. As shown by 
Rees \l2'i J. (see also Drozdova & Panchenko Sari \l6A l) in the case of an expanding emitting shell an observer 
simultaneously detects radiation produced in different moments of time (thus, with different temperatures) on the 
ellipsoidal or egg-like surface. The integration over this surface can give the same effect as the integration over time 



done in this paper, hut we do not perform this here because the result strongly depends on the unknown geometry of 
the emitting surface" (Blinnikov, Kozyreva & Panchenko 18]). This treatment which they outhne but they discard 
due to the difficulty of defining the geometry of the EQTS is exactly what we have done. Our treatment has only 
one free parameter and can fit the data of GRBs in a range between a few seconds all the way up to 10^ seconds. 
There is a basic observational feature between our treatment and the one by Blinnikov, Kozyreva & Panchenko TS]: 
their instantaneous spectral distribution has necessarily to be a blackbody one, while in our case is represented by an 
integration over an infinite number of planckian spectra, weighted by appropriate relativistic transformations, each one 
corresponding to a different viewing angle in the past light cone of the observer. The difference between such unique 
spectra should be simply discernible using the observations of XMM, Chandra, and of future space observatories. 

XXV. THE OBSERVATION OF THE IRON LINES IN GRB 991216: ON A POSSIBLE 

GRB-SUPERNOVA TIME SEQUENCE 

We have seen in the previous sections how the time structure of the E-APE gives information on the composition of 
the interstellar matter at distances of the order of 5 x 10^^ cm from the source. We w ould like now to point out that 
the data on the iron lines from the Chandra sat ellite on the GRB 9912 16 (Piro et al., jl20l |') and similar observations 
from other sources (Amati et al. Piro et al. |119| |. Piro et al., |l2fl| |) make it possible to extend this analysis to a 
larger distance scale, possibly all the way out to a few light years, and consequently probe the distribution of stars in 
the surroundings of the newly formed EMBH. 

Most importantly, the se co nsiderations lead to a new paradigm for the interpretation of the supernova-GRB corre- 
lation (see Ruffini et al. |l4(ij ). Indeed a correlation between the occurrence of GRBs and supernova events exists and 
has b een e stablished by t he works of Bloom e t al. [ l9l| . Galama et al. [5l,|5^|63|, Kulkarni et al. [s^, Piran |116| |. Piro 
et al. |ll8j . Rhoads 'l3l], van Paradijs et al. |179| |. 

Such an association has been assumed to indicate that GRBs are generated by supernova explosions (see e.g. 
Kulkarni et al. 85]). In turn, such a point of view has implied further consequences: the optical and radio data of the 
supernova have been attributed to the GRB afterglow, and many theorists have tried to encompass these data and 
explain them as a genuine component of the GRB scenario. 

We propose instead an alternative point of view implying a very clear distinction between the GRB phenomenon 
and the supernova: if relativistic effects presented in the RSTT paradigm are properly taken into account, then a 
kinematically viable explanation can be given of the supernova-GRB association. We still use GRB 991216 as a 
prototypical case. 

The G RB- Supernova Time Sequence paradigm, which we have indicated for short as GSTS paradigm (see Ruffini et 
al. |l46l ]). states that: A massive GRB-progenitor star Pi of mass Mi undergoes gravitational collapse to an EMBH. 
During this process a dyadosphere is formed and subsequently the P-GRB and the E-APE are generated in sequence. 
They propagate and impact, with their photon and neutrino components, on a second supernova-progenitor star P2 
of mass M2. Assuming that both stars were generated approximately at the same time, we expect to have M2 < Mi. 
Under some .special conditions of the thermonuclear evolution of the supernova-progenitor star P2, the collision of the 
P-GRB and the E-APE with the .star P2 can induce its supernova explosion. 

Especially relevant to our paradigm are the following data from the Chandra satellite (see Piro et al., |l2Cl| ): 

1. At the arrival time of 37 hr after the initial burst there is evidence of iron emission lines for GRB 991216. 

2. The emission lines are present during the entire observation period of 10^ s. The iron lines could also have been 
produced earlier, before Chandra was observing. Thus the times used in these calculations are not unique: they 
do serve to provide an example of the scenario. 

3. The emission lines appear to have a peak at an energy of 3.49 ± 0.06 keV which, at a redshift z — 1.00 ± 0.02 
corresponds to an hydrogen-like iron line at 6.97 keV at rest. This source does not appear to have any significant 
motion departing from the cosmological flow. The iron lines have a width of 0.23 keV consistent with a radial 
velocity field of 0.1c. The iron lines are only a small fraction of the observed flux. 

On the basis of the explicit computations of the different eras presented in the above sections, we make three key 
points: 

1. An arrival time of 37 hr in the detector frame corresponds to a radial distance from the EMBH travelled by the 
ABM pulse of 3.94 x 10^^ cm in the laboratory frame (see Tab.QJ. 

2. It is likely that a few stars are present within that radius as members of a cluster. It has become evident from 
observations of dense clusters of star-forming regions that a stellar average density of typically lO^pc"'^ (Beck 




Figure 45: A qualitative simplified space-time diagram (in arbitrary units) illustrating the GSTS paradigm. The EMBH, 
originating from the gravitational collapse of a massive GRB-progenitor star Pi, and the massive supernova-progenitor star 
P2-neutron star (P2-NS) system, separated by a radial distance Dp^, axe assumed to be at rest in in the laboratory frame. 
Their worldlines are represented by two parallel vertical lines. The supernova shell moving at 0.1c generated by the P2-NS 
transition is represented by the dotted line cone. The solid line represents the motion of the pulse, as if it would move with 
an "effective" constant gamma factor 71 during the eras reaching the condition of transparency. Similarly, another "effective" 
constant gamma factor 72 < 71 applies during era IV up to the collision with the P2-NS system. A third "effective" constant 
gamma factor 73 < 72 occurs during era V after the collision as the nonrelativistic regime of expansion is reached. The dashed 
lines at 45 degrees represent signals propagating at speed of light. 



et al. 9"|) should be expected. There is also the distinct possibility for this case and other systems that the stars 
Pi and P2 are members of a binary system. 

3. The possible observations at different wavelengths of the supernova crucially depend on the relative intensities 
between the GRB and the supernova as well as on the value of the distance and the redshift of the source. In 
the present case of GRB 991216, the expected optical and radio emission from the supernova are many orders of 
magnitude smaller than the GRB intensity. The opposite situation will be encountered in GRB 980425 (Ruffini 
et al. [HI). 

In order to reach an intuitive understanding of these complex computations we present a schematic very simplified 
diagram (not to scale) in Fig. 021 

We now describe the sequence of events and the specific data corresponding to the GSTS paradigm: 

1. The two stars Pi and P2 are separated by a distance Dp^ = 3.94 x 10^^ cm in the laboratory frame, see 
Fig. 1451 Both stars are at rest in the inertial laboratory frame. At laboratory time t — and at comoving time 



r = 0, the gravitational collapse of the GRB-progenitor star Pi occurs, and the initial emission of gravitational 
radiation or a neutrino burst from the event then synchronizes this event with the arrival times = at the 
supernova-progenitor star P2 and if = for the distant observer at rest with the detector. The electromagnetic 
radiation emitted by the gravitational collapse process is instead practically zero, due to the optical thickness 
of the material at this stage (Bianco et al. see Tab.Hl). 

2. From Tab.Hl at laboratory time ti = 6.48 x 10^ s and at a distance from the EMBH of Di = 1.94 x 10^'* cm, the 
condition of transparency for the PEMB pulse is reached and the P-GRB is emitted (see section llXjl . This time 
is recorded in arrival time at the detector tf^ = 8.41 x 10^^ s, and, at P2, at iai = 4.20 x 10^^ s. The fact that 
the PEMB pulse in an arrival time of 8.41 x 10~^ s covers a distance of 1.94 x 10^"* cm gives rise to an apparent 
"super luminal" effect. This app arent paradox can be straightforwardly explained by introducing an "effective" 
gamma factor, see Ruffini et al. |146| . 

3. At laboratory time t = 1.73 x 10^ s and at a distance from the EMBH of 5.18 x 10^^ cm in the laboratory frame, 
the peak of the E-APE is reached which is recorded at the arrival time ta ~ 9.93 s at P2 and = 19.87 s at the 
detector. This also gives rise to an apparent "superluminal" effect. 

4. At a distance Dp,^ = 3.94 x 10^^ cm, the two bursts described in the above points 2) and 3) collide with the 
supernova-progenitor star P2 at arrival times = 4.20 x 10^^ s and ta ~ 9.93 s respectively. They can then 
induce the supernova explosion of the massive star P2. 

5. The associated supernova shell expands with velocity 0.1c. 

6. The expanding supernova shell is reached by the ABM pulse generating the afterglow with a delay of = 18-5 hr 
in arrival time following the arrival of the P-GRB and the E-APE. This time delay coincides with the interval 
of lab orato ry time separating the two events, since the P2 is at rest in the inertial laboratory frame (see Ruffini 
et al. 'l46]). The ABM pulse has travelled in the laboratory frame a distance Dp^ — Di ~ Dp^ = 3.94 x 10^^ cm 
in a laboratory time ^2 — — ^2 = 1-32 x lO'' s (neglecting the supernova expansion). 

The collision of the pulse with the supernova shell occurs at 7 ~ 4.0. By this time the supernova shell has reached 
a dimension of 1.997 x lO^"' cm, which is consistent with the observations from the Chandra satellite. 

In these considerations on GRB 991216 the supernova remnant has been assumed to be close to but not exactly 
along the line of sight extending from the EMBH to the distant observer. If such an alignment should exist for other 
GRBs, it would lead to an observation of iron absorption lines as well as to an increase in the radiation observed 
in the afterglow corresponding to the crossing of the supernova shell by the ABM pulse. In fact, as the ABM pulse 
engulfs the baryonic matter of the remnant, above and beyond the normal interstellar medium baryonic matter, the 
conservation of energy and momentum implies that a larger amount of internal energy is available and radiated in 
the process (see section IXIIip . This increased energy- momentum loss will generally affect the slope of the afterglow 
decay, approaching more rapidly a nonrelativistic expansion phase (details are given in section rXVIII|l . 

It is quite clear that as soon as the relativistic transformations of the RSTT paradigm are duly taken into account, 
the sequence of events between the supernova and the GRB occur rences are exa ctly the opposite of the one postulated 
in the so-called "supranova" scenario (Vietri & Stella |l8ll II8 j |. Vietri et al. |l83j 'l. This can be considered a very 
appropriate pedagogical example of how classical nonrelativistic applied to ultrarelativistic regimes can indeed subvert 
the very causal relation between events. 

If we now turn to the possibility of dynamically implementing the scenario, there are at least three different 
possibilities: 

1. Particularly attractive is the possibility that a massive star P2 has rapidly evolved during its thermonuclear 
evolution to a white dwarf (see e.g. Chandrasekhar |2^). It it then sufficient that the P-GRB and the E- 
APE implode the star sufficiently as to reach a central density above the critical density for the ignition of 
thermonuclear burning. Consequently, the explosion of the star P2 occurs, and a significant fraction of a solar 
mass of iron is generated. These configurations are currently generally considered precursors of some type I 
supernovae (see e.g. Filippenko [^Jl and references therein). 

2. Alternatively, the massive star P2 can have evolved to the condition of being close to the point of gravitational 
collapse, having developed the formation of an iron-silicon core, type II supernovae. The above transfer of 
energy momentum from the P-GRB and the E-APE may enhance the capture of the electrons on the iron nuclei 
and consequently decrease the Fermi energy of the core, leading to the onset of gravitational instability (see e.g. 
Bethe P- 270 and foUowings). Since the time for the final evolution of a massive star with an iron-silicon 
core is short, this event requires a well tuned coincidence. 



3. The pressure wave may trigger massive and instantaneous nuclear burning process, with corresponding changes 
in the chemical composition of the star, leading to the collapse. 

The GSTS paradigm has been applied to the case of the GRB 980425 - SN1998bw which, with a red shift of 0.0083, 
is one of the closest and weaker GRBs observed. In this case, the radio and the optical emission of the supernova is 
distinctively observed. For this particular case, the EMBH appears to have a significantly lowe r value of the parameter 
^ and the validity of the GSTS paradigm presented here is confirmed (see Ruffini et al. |l52j V 



XXVI. GENERAL CONSIDERATIONS ON THE EMBH FORMATION 



Before concluding let us consider the problem of the EMBH formation. Such a problem has been debated for many 
years since the earliest discussions in 1970 in Princeton and has been finally clarified and addressed in general terms 
to justify the plausibility of the hypothesis in RufRni |l39l |. There has been a basic change of paradigm. All the 
cons iderations on the electric charge of stars were traditionally directed, following the classical work by Shvartsman 
all the way to the fundamental book by Punsly {126 j , to the presence of a net charge on the star surface in a 
steady state condition. The star can be endowed of rotation and magnetic field and surrounded by plasma, like in 
the case of Goldreich & Julian 70], or, in the case of absence of both magnetic field and rotation, the e lectrostatic 
processes can be related to the depth of the gravitational well, like in the treatment of Shvartsman |l74j . However, 
in neither cases it is possible to reach the condition of the overcritical field needed for pair creation nor have the 
condition of no baryonic contamination discussed in sections IHII IVHI and essential for the dyadosphere formation. 
The basic conceptual point is that GRBs are maybe the most violent transient phenomenon occurring in the universe 
and so the condition for the dyadosphere creation have to be searched in a transient phenomenon. The solution is 
related to the most transient phenomenon occurring in the life of a star: the process of gravitational collapse. 

Having acquired such a fundamental understanding, the next step is to estimate the amount of polarization needed 
in order to reach the fully relativistic condition 

^ ^ . (182) 
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Recalling that the charge to mass ratio of a proton is g^/ (^mp\/(fj = 1.1 x 10^^, it is enough to have an excess of 

one quantum of charge every 10^® nucleons in the core of the collapsing star to obtain an extreme EMBH after the 
occurrence of the gravitational collapse. Physically this means that we are dealing with a process of charge segregation 
between the core and the outer part of the star which has the opposite sign of net charge in order to enforce the 
overall charge neutrality condition. We here emphasize the name "charge segregation" instead of the name "charge 
separation" in order to contrast a very mild charge surplus created in different part of the star, keeping the overall 
charge neutrality, from the much more extreme condition of charge separation in which all the charges of the atomic 
component of the star are separated. It is indeed reassuring that such a core, endowed with charge segregation, is 
indeed stable with respect to the Fermi-Chandrasekhar criteria for the stability of self-gravitating stars duly extended 
from the magnetic to the electric case: the electric energy of such a core is consistently smaller than its gravitational 
energy (see Boccaletti et al. 0|). 

Such a condition of charge segregation between the core and the oppositely charged star surface layer can be reached 
under a very large number of physical conditions. We consider, for simplicity, one of the oldest example: the one of a 
star endowed with both a magnetic field and rotation. It is proved that a typical magnetic field expected for the ISM 
is Bo ^ 10^^ G (Ferriere [HJl). We further assume, consistently with the data which we have acquired and verified 
in the present article (see sections IXIIII |XIX|) . that also in the galaxy where GRB 991216 occurred the ISM has an 
average density of Uism = 1 proton/cm^. From this value of density we have that an ISM cloud with mass M ^ IOMq 
occupies a sphere of radius Ro ~ 1.4 x 10^^ cm. If this sphere collapse to a star with radius R = Rq, from the flux 
conservation we obtain that it is enough for this star to rotate with the most reasonable angular speed 



to conclude that the progenitor star core is endowed of a charge to mass ratio equal to ^. In the extreme case 
of Ea. (|182|) we have ^ = 1 and so the angular speed is ~ 1.1 x 10~^ radjs — i.e. one round in 1.5 hr — and 
correspondingly we have smaller Q values for ^ < 1 (see Boccaletti et al. jljl). Clearly the overall neutrality is 
guaranteed by the oppositely charged baryonic matter which is the one measured by the B parameter in the EMBH 
model (see sections IVIIII IIX)l . The smallness of the B value clearly points to the absence of an extended envelope of 
the progenitor star. 



The formation process of such an electromagnetised progenitor star will be clearly affected by the presence of differ- 
ential rotation, the consequent amplification of the magnetic field and a variety of magnetohydrodynamical problems 
which will affect somewhat the simplicity of the heuristic Ea. H183(l . Similarly the process of gravitational collapse 
of such a progenitor star endowed with rotation will lead to complex phenomena of "gravitationally induced electro- 
magnetic radiation" (Johnston et al. jsj) and of "electromagnetically induced gravitational radiation" (Johnston et 
al. which will tend to reduce both the eccentricity and the angular velocity of the collapsing core. The general 
outcome of gravitational collapse will be a Kerr-Newmann spacetime. I t is interesting that such a general case will 
break the degeneracy in (^,^) described in section |X] (see Ruffini et al. |153| ). In this article we have addressed the 
much simpler case of a solution in which (cL) / (GAf^) ^ 1 and the treatment can be well approximated by a collapse 
described by a Reissner-Nordstrom geometry. 

In addition to this scenario, based on the role of magnetic field and rotation, we are as well pursuing the possible 
generation of the charge segregation by quantum effects at the surface of the Fermi semi-degenerate core. In this 
framework, it is particularly interesting to consider the purely electric analog of the Chandrasekhar & Fermi 
paper on the gravitational stability of self-gravitating magnetized stars. The stability condition, based on the virial 
theorem, is simply that the Coulomb energy of the inner core of a charged star should be smaller or equal than the 
gravitational energy of the star (Boccaletti et al. ^^). Previous to the collapse, the gravitational energy can be 
much smaller than the rest energy of the star and be amplified dur ing the process of gravitational collapse reaching 
overcritical intensity of the electric field (see Fig. 1461 and Ruffini jl89| |). It is interesting that the Chandrasekhar-Fermi 
inequality just leads to an extreme Reissner-Nordstrom solution. 

In both these cases the Reissner-Nordstrom geometry appears indeed to be the relevant model for GRB 991216 as 
discussed in the previous sections. We shall return to non sp heric al configuration in forthcoming publications and/or 
when requested by observational evidence (see Ruffini et al. |15,'^ ). 



XXVII. SOME PROPAEDEUTIC ANALYSIS FOR THE DYNAMICAL FORMATION OF THE EMBH 

While the formation in time of the dyadosphere is the fundamental phenomena we are interested in, we can get 
an insight on the issue of gravitational collapse of an electrically charged star core studying in details a simplified 
model, namely a thin shell of charged dust. In De la Cruz & Israel [s^, Israel js^] it is shown that the problem of a 
collapsing charged shell in general relativity can be reduced to a set of ordinary differential equations. We reconsider 
here the following relativistic system: a spherical shell of electrically charged dust which is moving radially in the 
Reissner-Nordstrom background of an already formed nonrotating EMBH of mass Mi and charge Qi, with Qi < Mi. 

The world surface spanned by the shell divides the space-time into two regions: an internal one A4- and an external 
one A^_|-. The line element in Schwarzschild like coordinate is (Cherubini et al. |2^) 

2 _ / -f+dtl + U'dr^ + r^dn^ in M+ ,^ 
'^'^ ~\ - f_de + fl'dr^ + r^dn^ in M- ' ^ ' 

where /+ = 1 — ^ + /_ = 1 — + ^ and i_ and are the Schwarzschild- like time coordinates in 
and respectively. M is the total mass-energy of the system formed by the shell and the EMBH, measured by an 
observer at rest at infinity and Q ^ Qo + Qi is the total charge: sum of the charge Qo of the shell and the charge Qi 
of the internal EMBH. 

Indicating by R the radi us of the shell and by T± its time coordinate, the equations of motion of the shell become 
(Rufiini & Vitaghano [E^) 



^ dr ) A L' 



1 (m - Ml - M - g - ' U (R) , (185) 



M'o \ 2fl 2R R 

dr - Mof±{R) J'^-'l + 2R 2R R )' ^ ' 

where Mq is the rest mass of the shell and r is its proper time. Eas. H185ll86() (together with Eg. (|1 8411 ') completely 

describe a 5-parameter (M, Q, Mi, Qi, Mq) family of solutions of the Einstein-Maxwell equations. Note that 
Eas. H185ll86() imply that 

M - Ml - ^ - > (187) 



holds for R> M + ^ AP - ii Q < M and for R > Mi + ^ M^ - Qj if Q > M. 




Figure 46: Quantitative description of the gravitational collapse to a neutron star and to a black hole of the core of a rotating 
progenitor rotating. The core is estimated to have a mass equal to 3Mq, to have an initial radius r = r© and a rotation 
period of 15 days. Although the initial rotational energy is of the order of 10~^^ of the total energy, the total rotational energy, 
in principle extractable, of the rotating black hole can be as high as of the order of 29%. On the lower- right side the same 
considerations are applied to the case of a neutral star formed by a core oppositely charged from its outermost envelope. The 
core is expected to have a mass of 3Mq, a radius equal to tq and electromagnetic energy Q/M = 0.1. Although the initial 
Coulomb energy is only ~ lO"'^ of the total energy, which is in turn hundred times smaller than the gravitational energy, the 
final Coulomb energy can be as high as 2.5 x 10""^ of the total energy. In both cases, the amplification of the rotational energy 
and of the Coulomb energy, which indeed are the only two extractable forms of energy from a black hole, is due to the process 
of gravitational collapse. 



For astrophysical applications (Ruffini et al. the trajectory of the shell R = R (1+) is obtained as a function 

of the time coordinate relative to the space-time region Ai + . In the following we drop the + index from T^. From 
Eas. p85ll86|) we have 



dT ~ dr dT ~ O 

where 



F^/+(i?,) = l-2f + 2!, (189) 

^ ^ r - (190) 
T^M^. (191) 



Since we are interested in an imploding shell, only the minus sign case in (I188II will be studied. We can give the 
following physical interpretation of F. If M — Mi > AIq, T coincides with the Lorentz 7 factor of the imploding shell 
at infinity; from Ea. H188|l it satisfies 



r = 



yi-(#)R=„ 



> 1. 



(192) 



When AI — Mi < Mq then there is a turning point R* , defined by ^|^_^. = 0. In this case F coincides with the 
"effective potential" at R* : 



T = ^/flW) + Mo 
The solution of the differential equation H188|l is given by 

dT=- 



-1 i-Ml 4_ Qo. 

2R' ^ 2R* 



5iQo 
R' 



< 1. 



n 



dR. 



(193) 



(194) 



The functional form of the integral H194() crucially depends on the degree of the polynomial P (R) = R^ (j}? — F) , 
which is gcncrically two, but in special cases has lower values. We therefore distinguish the following cases: 



1. M — Mo + Mi; Qi = Mi; Q — M: P (R) is equal to 0, we simply have 

R{T) = const. 

2. M = Mo + Mi; M^ - = M^ - Qj; Q ^ M: P (i?) is a constant, we have 



T = const 
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(196) 



3. M = Mo + Mi; M^ - ^ M^ - Ql: P (R) is a first order polynomial and 



T = const + 2R\/n'^ - F 



MqR 

3(M2-Q2-Af2^Q2) 



{Mg+Q^-QlY-9MMo{Mi+Q^-Ql) + 12M^M^+2Q^M^ 
3{M2-Q''-Mf+Qiy 

r^arctanh(;^^^£) 



„2 , 
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(197) 



where r2± = {r±). 
4. M Mo + Mi: P (R) is a second order polynomial and 



T = const 
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Figure 47: Left) Collapse curves in the plane (T, R) for M — 2OAf0 and for different values of the parameter ^. The asymptotic 
behavior is the clear manifestation of general relativistic effects as the horizon of the EMBHf is approached. Right) Electric 
field behaviour at the surface of the shell for M — 20Mq and for different values of the parameter ^. The asymptotic behavior 
is the clear manifestation of general relativistic effects as the horizon of the EMBH is approached. 



Of particular interest is the time varying electric field £u = on the external surface of the shell. In order to 
study the variability of Sr with time it is useful to consider in the tridimensional space of parameters (i?, T, Sn) 

the parametric curve C : (^R = A, T = T{X), Er — . In astrophysical applications (Ruffini et al. jl59l |) we are 
specially interested in the family of solutions such that ^ is when R — oo which implies that F = 1. In Fig. 1471 we 
plot the collapse curves in the plane {T,R) for different values of the parameter ^=-^,0<^<l. The initial data 
(To, i?o) are chosen so that the integration constant in equation H197|l is equal to 0. In all the cases we can follow the 
details of the approach to the horizon which is reached in an infinite Schwarzschild time coordinate. In Fig. 1471 we plot 
the parametric curves C in the space (i?, T, Er) for different values of ^. Again we can follow the exact asymptotic 
behavior of the curves C, £r reaching the asymptotic value ^-. The detailed knowledge of this asymptotic behavior 

is of great relevance for the observational properties of the EMBH formation (see e.g. Ruffini & Vitagliano |l5fil |). 
In the case of a shell falling in a fiat background [Mi = Qi = 0) Ea. H185|l reduces to 

(fr=ife(^^+f (199) 

Introducing the total radial momentum P = Mqu^ = Mq ^ of the shell, we can express the kinetic energy of the shell 
as measured by static observers in as T = —M^u^^t — Mq = ^ P"^ + A/q — Mq. Then from equation H199|l we 



have 



+ % + \/P' + M^ ^ Afo +r- H + (200) 



where we choose the positive root solution due to the constraint (|187|l . Ea. (|200() is the mass formula of the shell, 
which depends on the time-dependent radial coordinate R and kinetic energy T. If M > Q, an EMBH is formed and 
we have 

M = Mo + T+~^ + ^, (201) 



where = T {r+) and r+ = M + ^ NP- — jg the radius of external horizon of the EMBH. We know from the 
Christodoulou-Rufhni EMBH mass formula that 



M = Mi„ + ^, (202) 



so it follows that 



M„. = Mo-^+T+, (203) 



namely that M\„ is the sum of only three contributions: the rest mass Afg, the gravitational potential energy and the 
kinetic energy of the rest mass evaluated at the horizon. M\„ is independent of the electromagnetic energy, a fact 
noticed by Bekenstein (Bekenstein [Tol|'). We have taken one further step here by identifying the independent physical 
contributions to Mi„. 

Next we consider the physical interpretation of the electromagnetic term ^ , which can be obtained by evaluating 
the conserved Killing integral 
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(204) 



where T,f is the space-like hypersurface in described by the equation tj^- = t — const, with dY^ as its surface ele- 
ment vector and where T^™'' = — ^ {Ff^PPpi, -\- \g^vF'"^ Fp„^ is the energy-momentum tensor of the electromagnetic 
field. The quantity in Ea. (|204|l differs from the purely electromagnetic energy 



JT.t " J R 



where n^|^ = f^^^'^^^ is the unit normal to the integration hypersurface and = /+• This is similar to the analogous 



situation for the total energy of a static spherical star of energy density e within a radius R, m{R) — A-k dr r^e. 



which differs from the pure matter energy mp (i?) = An jj^ dr^fg^Pt by the gravitational energy (see Misner, Thorne, 
& Wheeler |l03j ). Therefore the term ^ in the mass formula H200|l is the total energy of the electromagnetic field 
and includes its own gravitational binding energy. This energy is stored throughout the region , extending from R 
to infinity. 

We now turn to the problem of extracting the electromagnetic energy from an EMBH see (see Christodoulou & 
Ruffini [8(il|). We can distinguish between two conceptually physically different processes, depending on whether the 

electric field strength f = is smaller or greater than the critical value £c = ■ Here rUe and e are the mass and 
the charge of the electron. As already mentioned in this paper an electric field £ > £c polarizes the vacuum creating 

electron-positron pairs (see Heisenberg & Euler IzSl). The maximum value = -3- of the electric field around an 

''+ 

EMBH is reached at the horizon. We then have the following: 



1. For < £c the leading energy extraction mechanism consists of a sequence of discrete elementary decay 
processes of a particle into two oppositely charged particles. The condition < Ec implies 



where Ac is the Compton wavelength of the electron. Denardo & RuHini [Sg and Denardo et al. 37] have 
defined as the effective ergosphere the region around an EMBH where the energy extraction processes occur. 
This region extends from the horizon r+ up to a radius 



'Ecre — -72- 



(206) 



The energy extraction occurs in a finite number A'pd of such discrete elementary processes, each one corre- 
sponding to a decrease of the EMBH charge. We have 

iVpD ^ f . (207) 

Since the total extracted energy is (see Eq. H2()2|l l = we obtain for the mean energy per accelerated 
particle {E)^^ = ^ 

(i;) = ^ = 1 i ^ TOeC^ ~ hi^^ rnec^^ (208) 



which gives 



X lo^iey if 1^ < 106 



M, 







One of the crucial aspects of the energy extraction process from an EMBH is its back reaction on the irreducible 
mass expressed in Christodoulou & Ruffini ■ Although the energy extraction processes can occur in the entire 
effective ergosphere defined by Eq. H206() , only the limiting processes occurring on the horizon with zero kinetic 
energy can reach the maximum efficiency while approaching the condition of total reversibility (see Fig. 2 in 
Christodoulou & Ruffini for details). The farther from the horizon that a decay occurs, the more it increases 
the irreducible mass and loses efficiency. Only in the complete reversibility limit (Christodoulou & Ruffini 's^) 
can the energy extraction process from an extreme EMBH reach the upper value of 50% of the total EMBH 
energy. 

2. For £4. > £c the leading extraction process is a collective process based on an electron-positron plasma generated 
by the vacuum polarization, (see Fig.^ as discussed in section ITTll The condition > £c implies 

GM/c^ ( __e_\ ^ o 1 n-6 M 



^i^j -2.10-^<e<l. (210) 

This vacuum polarization process can occur only for an EMBH with mass smaller than 2 • 10^ Mq. The electron- 
positron pairs are now produced in the dyadosphere of the EMBH, (note that the dyadosphere is a subregion 
of the effective ergosphere) whose radius r^s is given in Eq.©. We have Vds ^ '"Ecrg- The number of par- 
ticles created and the total energy stored in dyadosphere are given in Eas. ()10ll2|l respectively and we have 
approximately 

^e^e--(5^)f, (211) 
Edya ^ ^ (212) 

The mean energy per particle produced in the dyadosphere (E)^^ = J^o^" is then 

(^^>d. - i (^) ^ ^ (213) 

which can be also rewritten as 



(^^>d. ^ ^ V MfelO^fcey . (214) 



Such a process of vacuum polarization, occurring not at the horizon but in the extended dyadosphere region 
(r-|_ < r < Tds) around an EMBH, has been observed to reach the maximum efficiency limit of 50% of the total 
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Figure 48: Space-time diagram of the collapse process leading to the formation of the dyadosphere. As the collapsing core 
crosses the dyadosphere radius the pair creation process starts, and the pairs thermalize in a neutral plasma configuration. 
Then also the horizon is crossed and the singularity is formed. 



mass-energy of an extreme EMBH (see e.g. Preparata et al. jl24i| 'l. The conceptual justification of this result 
follows from the present work: the e+e" creation process occurs at the expence of the Coulomb energy given by 
Eq. (|2()4|l and does not affect the irreducible mass given by Eq. H2()3(l . which indeed, as we have proved, does 
not depend of the electromagnetic energy. In this sense, SMii,- = and the transformation is fully reversible. 
This result will be further validated by the study of the dynamical formation of th e dyadosphere, which we have 
obtained using the present work and Cherubini et al. (see Ruffini et al. |l59l |). 

Let us now compare and contrast these two processes. We have 

TEcrg ^ (^) r (215) 
iVdya ^ (^) A^PD, (216) 
(^)dya-te)(^>PD- (217) 

Moreover we see (Eqs. (I209II . (|214|) ') that {E)pjj is in the range of energies of UHECR, while for ^ ~ 0.1 and 
M ~ IOMq, (i?)cis is in the gamma ray range. In other words, the discrete particle decay process involves a small 
number of particles with ultra high energies (~ lO^^'^eT^), while vacuum polarization involves a much larger number 
of particles with lower mean energies (~ lOMeV). 

Having so established and clarified the basic conceptual processes of the energetic of the EMBH, we are now ready 
to approach, using the new analytic solution obtained, the dynamical process of vacuum polarization occurring during 
the formation of an EMBH as qualitatively represented in Fig. 1481 The study of the dyadosphere dynamical formation 
as well as of the electron-positron plasma dynamical evolution will lead to the first possibility of directly observing 
the general relativistic effects approaching the EMBH horizon. 

Before closing we would like to emphasize once more a basic point: all the considerations presented in the description 
of the preceding eras are based on the approximations in the description of the dyadosphere presented in section UTTl 
This treatment is very appropriate in estimating the general dependence of the energy of the P-GRB, the kinetic 
energy of the ABM pulse and consequently the intensity of the afterglow, as well as the overall time structure of the 
GRB and especially the time of the release of the P-GRB in respect to the moment of gravitational collapse and its 
relative intensity with respect to the afterglow. If, however, is addressed the issue of the detailed temporal structure of 
the P-GRB and its detailed spect ral d istribution, the above dynamical considerations on the dyadosphere formation 
are needed (see also Ruffini et al. |l59j V In turn, this detailed analysis is needed if the general relativistic effects close 
to the horizon formation have to be followed. As expressed already in section. IXIII all general relativistic quantum 
field theory effects are encoded in the fine structure of the P-GRB. As emphasized in section |3 the only way to 
differentiate between solutions with same E^ya but different EMBH mass and charge is to observe the P-GRBs in the 
limit B 0, namely, to observe the short GRBs. 



XXVIII. CONTRIBUTION OF THE EMBH MODEL TO THE BLACK HOLE THEORY 



The aim of this section is to point out how the knowledge obtained from the EMBH model is of relevance also for 
the basic theory of black holes and further how very high precision verification of general relativistic effects in the 
very strong field near the formation of the horizon should be expected in the near future. 

We shah first see how Ea. H203|l for M{„, 

Mrr = Mo - ^ + r+ , (218) 

leads to a deeper physical understanding of the role of the gravitational interaction in the maximum energy extraction 
process of an EMBH. This formula can also be of assistance in clarifying some long lasting epistemological issue on 
the role of general relativity, quantum theory and thermodynamics. 

It is well known that if a spherically symmetric mass distribution without any electromagnetic structure undergoes 
free gravitational collapse, its total mass-energy A/ is conserved according to the Birkhoff theorem: the increase in the 
kinetic energy of implosion is balanced by the increase in the gravitational energy of the system. If one considers the 
possibility that part of the kinetic energy of implosion is extracted then the situation is very different: configurations 
of smaller mass-energy and greater density can be attained without violating Birkhoff theorem. 

We illustrate our considerations with two examples: one has found confirmation from astr qphy sical observations, 
the other promises to be of relevance for gamma ray bursts (GRBs) (see Ruffini & Vitagliano llSGj ). Concerning the 
first example, it is well known from the work of Landau _.86] that at the endpoint of thermonuclear evolution, the 
gravitational collapse of a spherically symmetric star can be stopped by the Fermi pressure of the degenerate electron 
gas (white dwarf). A configuration of equilibrium can be found all the way up to the critical number of particles 

TVeHt = 0.775^, (219) 

where the factor 0.775 comes from the coefficient ^'^t^ of the solution of the Lane-Emden equation with polytropic 

index n = 3, and mpi = is the Planck mass, mo is the nucleon mass and /i the average number of electrons per 
nucleon. As the kinetic energy of implosion is carried away by radiation the star settles down to a configuration of 
mass 

M = iVcritmo - U, (220) 

where the gravitational binding energy U can be as high as 5.72 x 10~^A'crit'7io. 

Similarly Gamov (see Gamow & Critchfield [bJ) has shown that a gravitational collapse process to still higher 
densitie s can be stopped by the Fermi pressure of the neutrons (neutron star) and Oppenheimer (Oppenheimer & 
Volkoff 'l06l|) has shown that, if the effects of strong interactions are neglected, a configuration of equilibrium exists 
also in this case all the way up to a critical number of particles 

A^erit = 0.398^, (221) 

where the factor 0.398 comes now from the integration of the Tolman-Oppenheimer- Volkoff equation (see e.g. Har- 
rison et al. Iz^). If the kinetic energy of implosion is again carried away by radiation of photons or neutrinos and 
antineutrinos the final configuration is characterized by the formula H220(l with U < 2.48 x lO^^A'critn^o- These 
considerations and the existence of such large values of the gravitational binding energy have been at the heart of the 
explanation of astrophysical phenomena such as red-giant stars and supcrnovae: the corresponding measurements of 
the masses of neutron stars and white dwarfs have been carried out with unprecedented accuracy in binary systems 
(Gursky & Ruffini 

From a theoretical physics point of view it is still an open question how far such a sequence can go: using causality 
nonviolating interactions, can one find a sequence of braking and energy extraction processes by which the density 
and the gravitational binding energy can increase indefinitely and the mass-energy of the collapsed object be reduced 
at will? This question can also be formulat ed in the mass-formula language of a black hole given in Christodoulou 
& Ruffini f30j (see also Ruffini & Vitagliano |l5(Tj ) : given a collapsing core of nucleons with a given rest mass-energy 
Mq, what is the minimum irreducible mass of the black ho le wh ich is formed? 

Following Cherubini et al. and Ruffini & Vitagliano (l56|, consider a spherical shell of rest mass Mq collapsing 
in a flat space-time. In the neutral case the irreducible mass of the final black hole satisfies the equation (see Ruffini 
& Vitaghano (156j 'l 

Afi„. = M = Mo - ^ + r+, (222) 
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Figure 49; Collapse curves for neutral shells with rest mass Mo starting at rest at selected radii R* computed by using the 
exact solutions given in Cherubini et al. [2^ . A different value of Mirr (and therefore of r+) corresponds to each curve. The 
time parameter is the Schwarzschild time coordinate t and the asymptotic behaviour at the respective horizons is evident. The 
limiting configuration Mirr = (solid line) corresponds to the case in which the shell is trapped, at the very beginning of its 
motion, by the formation of the horizon. 



where M is the total energy of the collapsing shell and T+ the kinetic energy at the horizon r+ . Recall that the area 
S of the horizon is (Christodoulou & Rufhni IsO]) 

S = 47rr^ levrM^.^ (223) 

where — 2Miii- is the horizon radius. The minimum irreducible mass M^^^'^^ is obtained when the kinetic energy 
at the horizon is 0, that is when the entire kinetic energy has been extracted. We then obtain the simple result 

m/™^ = fi. (224) 

We conclude that in the gravitational collapse of a spherical shell of rest mass Mq at rest at infinity (initial energy 
Mi = Mq), an energy up to 50% of Mqc^ can in principle be extracted, by braking processes of the kinetic energy. 
In this limiting case the shell crosses the horizon with T+ = 0. The limit in the extractable kinetic energy can 
further increase if the collapsing shell is endowed with kinetic energy at infinity, since all that kinetic energy is in 
principle extractable. 

In order to illustrate the physical reasons for this result, using the formulas of Cherubini et al. j23, we have 
represented in Fig.02|the world lines of spherical shells of the same rest mass Mq, starting their gravitational collapse 
at rest at selected radii R*. These initial conditions can be implemented by performing suitable braking of the 
collapsing shell and concurrent kinetic energy extraction processes at progressively smaller radii (see also Fig. I50|l . 
The reason for the existence of the minimum (|224|) in the black hole mass is the "self closure" occurring by the 
formation of a horizon in the initial configuration (thick line in Fig. I49|l . 

Is the limit Miir — > ^ actually attainable without violating causality? Let us consider a collapsing shell with 
charge Q. If M > Q an EMBH is formed. As pointed out in RuHini & Vitagliano |l5(Tj | the irreducible mass of the 
final EMBH does not depend on the charge Q. Therefore Eqs. (|222|l and (|224|l still hold in the charged case with 
r+ = M + a/ Af 2 — jj^ Fig.lSOlwe consider the special case in which the shell is initially at rest at infinity, i.e. has 
initial energy Mj = Mq, for three different values of the charge Q. We plot the initial energy Mi, the energy of the 
system when all the kinetic energy of implosion has been extracted as well as the sum of the rest mass energy and 

the gravitational binding energy — of the system (here R is the radius of the shell). In the extreme case Q = Mq, 
the shell is in equilibrium at all radii (see Cherubini et al. 28]) and the kinetic energy is identically zero. In all three 

cases, the sum of the extractable kinetic energy T and the electromagnetic energy ^ reaches 50% of the rest mass 
energy at the horizon, according to Eq. H224(l . 

What is the role of the electromagnetic field here? If we consider the case of a charged shell with Q ~ Mq, 
the electromagnetic repulsion implements the braking process and the extractable energy is entirely stored in the 
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Figure 50: Energetics of a shell such that Mi = Mq, for selected values of the charge. In the first diagram Q = 0; the dashed 
line represents the total energy for a gravitational collapse without any braking process as a function of the radius R of the 
shell; the solid, stepwise line represents a collapse with suitable braking of the kinetc energy of implosion at selected radii; the 
dotted line represents the rest mass energy plus the gravitational binding energy. In the second and third diagram Q/Mq = 0.7, 
Q/Mo = 1 respectively; the dashed and the dotted lines have the same meaning as above; the solid lines represent the total 
energy minus the kinetic energy. The region between the solid line and the dotted line corresponds to the stored electromagnetic 
energy. The region between the dashed line and the solid line corresponds to the kinetic energy of collapse. In all the cases 
the sum of the kinetic energy and the electromagnetic energy at the horizon is 50% of Mq. Both the electromagnetic and the 
kinetic energy are extractable. It is most remarkable that the same underlying process occurs in the three cases: the role of the 
electromagnetic interaction is twofold: a) to reduce the kinetic energy of implosion by the Coulomb repulsion of the shell; b) 
to store such an en ergy in the region around the EMBH. The stored electromagnetic energy is extractable as shown in Ruffini 
& Vitaghano [l5^ . 



electromagnetic field surrounding the EMBH (see Ruffini & Vitagliano .ISfij ) . In Ruffini & Vitagfiano |l5fil | we have 
outhned two different processes of electromagnetic energy extraction. We emphasize here that the extraction of 50% 
of the mass-energy of an EMBH is not specifically linked to the electromagnetic field but depends on three factors: 
a) the increase of the gravitational energy during the collapse, b) the formation of a horizon, c) the reduction of the 
kinetic energy of implosion. Such conditions are naturally met during the formation of an extreme EMBH but are 
more general and can indeed occur in a variety of different situations, e.g. during the formation of a Schwarzschild 
black hole by a suitable extraction of the kinetic energy of implosion (see Fig. 1491 and Fig. I50|l . 

Now consider a test particle of mass m in the gravitational field of an already formed Schwarzschild black hole of 
mass M and go through such a sequence of braking and energy extraction processes. Kaplan (Kaplan 83] ) found for 
the energy E of the particle as a function of the radius r 

E = my^l-^. (225) 

It would appear from this formula that the entire energy of a particle could be extracted in the limit r 2M . 
Such 100% efficiency of energy extraction has often been quoted as evidence for incompatibility between General 
Relativity and the second principle of Thermodynamics (see Bekenstein |ll) and references therein). J. Bekenstein 
and S. Hawking have gone as far as to consider General Relativity not to be a complete theory and to conclude 
that in order to avoid inconsistencies with thermodynamics, the theory should be implemented through a quantum 
description (Bekenstein [lll|. Hawking [t^). Einstein himself often expressed the opposite point of view (see e.g. 



Dyson gg). 

The analytic treatment presented in Cherubini et al. '2^1 can clarify this fundamental issue. It allows to express 
the energy increase i? of a black hole of mass Mi through the accretion of a shell of mass Mq starting its motion at 
rest at a radius R in the following formula which generalizes Eq. H225|) : 



£;eeM-Mi = -||+Mo^1-^, (226) 

where M = Mi + E is clearly the mass-energy of the final black hole. This formula differs from the Kaplan formula 
(|225|l in three respects: a) it takes into account the increase of the horizon area due to the accretion of the shell; b) 
it shows the role of the gravitational self energy of the imploding shell; c) it expresses the combined effects of a) and 
b) in an exact closed formula. 

The minimum value -Emin of E is attained for the minimum value of the radius R = 2M: the horizon of the final 
black hole. This corresponds to the maximum efficiency of the energy extraction. We have 



£;,„i„ ^-^+ Mo^l - ^ - -ilMTTk-) + ^-^0^1 - MT^, (227) 
or solving the quadratic equation and choosing the positive solution for physical reasons 



Eniin = i ( \/M^ + M2 - Ml ] . (228) 
The corresponding efficiency of energy extraction is 




\ 1 + 9^-1- (229) 



„ _ A/o--Emi, 
'/max — Mo 

which is strictly smaller than 100% for any given Mq ^ 0. It is interesting that this analytic formula, in the limit 
Ml <C Mq, properly reproduces the result of equation H224|l . corresponding to an efficiency of 50%. In the opposite 
limit Ml ^ Mq we have 

^max^l-i^. (230) 

Only for AIq 0, Eq. (|229|l corresponds to an efficienc y o f 100% and correctly represents the limiting reversible 
transformations introduced in Christodoulou & Ruffini It seems that the difficulties of reconciling General 

Relativity and Thermodynamics are ascribable not to an incompleteness of General Relativity but to the use of the 
Kaplan formula in a regime in which it is not valid. The generalization of the above results to stationary black holes 
is being considered. 



XXIX. CONCLUSIONS 



The EMBH theory has been here applied for the first time to fit the experimental data of GRB 991216. This 
process has given us the opportunity to rethink the entire GRB process in an unitary description starting from the 
moment of gravitational collapse all the way up to the latest phases of the afterglow. We have identified the three 
fundamental actors of the GRB phenomenon in: 

1- Edya- Having reanalyzed in section Hill the physics of the dyadosphere we have pointed out in Fig. El that the 
same value of E^ya can be obtained from an entire family of (/i, ^) parameters (i.e. Edya is degenerate in (/i, ^)). 
We have then shown in the reexamination of all the GRB eras that all the results depend only on the value of 
Edya and not on the particular value of (/i, <^) (see sections IVIIIIIXIXIIIIXIV)l . The only exception to this occurs 
in the era I (see section |VII|) which is the only one relevant for short GRBs. 

2. B. The crucial role played by the baryonic remnant of the progenitor star in determining the relative intensity 
ratio and the time delay between the P-GRB and the E-APE has been summarized already in the two Figs. lllllTTl 
in the introduction. 

3. ISM. The density Uism of the interstellar medium and its inhomogeneities appears to have a fundamental role 
in the intensity and the temporal substructures of the E-APE and the afterglow. 



The observational data agree with the predictions of the model on: 



1. the intensity ratio, 1.58 x 10 ^, between the P-GRB and the E-APE, which strongly depends on the parameter 
B; 

2. the absolute intensities for both the P-GRB and the E-APE, respectively 7.54 x 10^^ erg and 4.75 x 10^^; 

3. the arrival time of the P-GRB and the peak of the E-APE, respectively 8.41 x 10^^ s and 19.87 s; 

4. the power-law index n of the afterglow, predicted ntheo = —1.6 and observed rtofcs — —1-616 ± 0.067 (see 
sections IXVTTTllXTXll : 

5. the temporal structure of the E-APE and its correlation with the inhomogeneity in the ISM; 

6. the spectral distribution of the X-ray and 7-ray emission. 

Concerning the total energy of GRB 991216, Edya = 4.83 x 10^^ erg is found in the EMBH theory. This value is 
systematically larger than the ones quoted in the current literature by Panaitescu & Kumar ^lllj and by Halpern et 
al. llSl due to the fact that they respectively consider beaming angles of 3° — 4° and 6°. These considerations have 
been shown to be untenable in section IXIXI There is still a difference of ~ 28% between the total energy implied by 
the EMBH theory (4.83 x 10^^ erg) and the value quoted by Halpern {Edya = 6.7 x 10^'^ erg) in the case of spherical 
emission. We trust that this is a consequence of the underlying assumption of the spectral distribution of the radiation 
assumed by Halpern et al. (see e.g. Frail et al. which should be reassessed on the ground of our theoretical 

results (see also RufSni et al. |l53j V 

These results can certainly be considered the success of the EMBH theory. 

Before closing, we like to stress how GRBs, if duly theoretically interpreted, can open a main avenue of inquiring on 
totally new physical and astrophysical regimes. This program is very likely one of the greatest computational efforts 
in physics and astrophysics and cannot be actuated using shortcuts. 

From the point of view of fundamental physics new regimes are explored: 

1. The process of energy extraction from black holes. It is interesting that the analysis of GRBs has promoted a 
new effort in developing new theoretical tools for approaching the dynamical phase of collapse as expressed in 
section. IXXVIII These results have further clarified some basic issue related to the energy extraction process 
from bla ck ho le (see e.g. Ruffini |139| | '). It was already known from the definition of the ergosphere (see Ruffini & 
Wheeler |135| ') that the rotational energy extraction process do occur in an extended region around the horizon 
of a black hole. The fortunate situation that the energy extraction process in GRBs occurs in a condition of 
almost perfect spherical symmetry have allowed us to focus on the second fundamental parameter of black holes, 
namely the electric energy. The spherical symmetry has allowed as well to develop some powerful theoretical 
tools (see section IXXVIIII which have allowed to reach a better understanding of the role of kinetic energy of 
implosion in the process of gravitational collapse, in the storage of electromagnetic energy in the region around 
black holes and to establish as well a new upper limit in the energy extraction process in the gravitational 
collapse up to 50% of the initial rest mass of the system (see section IXXVIII|I . These results are of general 
validity and do transcend the work on the EMBH theory, although they are motivated by these researches. 
Interestingly this work, by giving a new expression for the efficiency of transforming gravitational energy into 
mechanical work (see section lXXVIIIll , has opened up a new opportunity of debating the relation between general 
relativity, thermodynamics and quantum theory, which is certainly one of the most profound and important 
topic of research in the entire realm of fundamental physics. 

2. The quantum and general relativistic effects of matter- antimatter creation near the black hole horizon. It is 
well known that one of the most important topics pursued in the last seventy years in physics has been the 
possibility, postulated by Sauter, Heisemberg, Euler, Schwinger to create matter-antimatter from the vacuum. 
In order to have the first experimental and observational evidence for this phenomenon, three major approaches 
are being followed: 

a) In central collisions of heavy ions near the C oulomb barri er, as first proposed in Ger stein & Y. B. Zel'dovich |^ 
|6ff| (see also Popov & Rozhdestvenskaya |l2l| . Popov |l22j . Zel'dovich & Popov |l92j |'). Efforts in experimentally 
implementing this idea at GSI were made since early 80's. Despite some apparently encouraging result (Schweppe 
et al. I l7l| ). such efforts have failed so far due to the small contact time of the colliding ions (see e.g. Ahmad 
et al. 0, Bar et al. Ganz et al. Heinz et al. [l^li Leinberger et al. [13 )■ Typically the electromagnetic 
energy involved in the collisions of heavy ions with impact parameter 10~^^cm is Ei ~ lOj^erg. 

b) At the focus of an X-ray free electron laser (XFEL) (see Ringwald ^132*1, Roberts et al. |133| and references 
therein). This idea will be possibly testable at DESY, where the XFEL is part of the design of the collider 
TESLA, as well as at SLAG, where the so-called Linac Coherent Light Source (LCLS) has been proposed. The 
electromagnetic energy at the focus of an XFEL is E2 ^ lO^erg concentrated in a region of linear extension 



h - 10-8 cm (Ringwald [Tp l. 

c) Around an electromagnetic black hole (EMBH) (Damour & R.uf Bni 1 3 J. Preparata et al. \l2?l Il2^ 'l. giving 
rise to the observed phenomenon of GRBs (see e.g. Rufhni et al. |144L Il45. 146. ,148] ). The electromagnetic 
energy of an EMBH of mass M ~ IOMq and charge Q~ 0.1M/-\/G is E3 ^ 10^ ^ ergs and it is deposited in a 
region of linear extension ^3 ~ lO^m (Preparata et al. |l24j . Ruffini & Vitagliano |l56j '). 

There is the very distinct possibility that in this race the success will be reached by the observations in relativistic 
astrophysics more than from the high energy experiments on the Earth. This will be certainly a splendid success 
which will be only second to the discovery of Helium first in the stars and then on the Earth! Quite apart from 
the discovery in itself, the detection of vacuum polarization in the astrophysical settings presents distinctively 
new physical phenomena as Rufhni et al. |l60| . The very important topic to be covered in the forthcoming 
months is the study of the dynamical phase of gravitational collapse and to follow the effects of such process of 
vacuum polarization in the dynamical phase. It will be also impo rtan t to follow the development of this process 
all the way to the emission of the P-GRB (Ruffini & Vitagliano |l58j |'). 

3. The physics of ultrarelativisitc shock waves with Lorentz gamma factor 7 > 100. We are expecting much 
progress in this topic from the understanding of th e instantaneous spectrum of GRBs. Some preliminary results 
along this line are presented in Ruffini et al. |l5Clj |. See also section IXXIVI 

From the point of view of astronomy and astrophysics also new regimes are explored: 

1. The occurrence of gravitational collapse to a black hole from a critical mass core of mass M > IOM0, which 
clearly differs from the values of the critical mass encountered in the study of stars "catalyzed at the endpoint 
of thermonuclear evolution" (white dwarfs and neutron stars). 

2. The extremely high efficiency of the spherical collapse to a black hole, where almost 99.99% of the core mass 
collapses leaving negligible remnant. The EMBH theory offers an unprecedented tool in order to map with great 
accuracy all the matter distribution around the newly formed EMBH from the horizon all the way to the ISM. 
This concept was pioneered by Dermer & Mitman who proposed to use GRB sources as "tomographic images 
of the density distributions of the medium surrounding the sources of GRBs" . It is important to emphasize 
that the very precise reading of the matter distribution encoded in the data of the P-GRB, the E-APE and the 
afterglow in GRB 991216 is in marked disagreement with the matt er di stribution postulated by the "coUapsar" 
scenario (see MacFadyen & Woosley Paczyhski |llOj . Woosley |l9n| |). This conclusion is evidenced not only 
by the absence of beaming already mentioned above, but also for the paucity of the baryonic matter encountered 
by the PEM pulse in its way out from the EMBH. There is no evidence for the presence either of a baryonic 
disk component nor of a conspicuous baryonic remnant. We actually have B — 3.0 x 10^'^. Unlike the case of 
formation of a neutron star, the mass of the remnant of the progenitor star is very small indeed. This mass, 
determined by B, is very accurately inferable from the relative intensity and temporal distance between the 
P-GRB and the E-APE (see above). In the present case we have Mb ~ 8.1 x IQ-^Mq. The presence of the 
remnant is also important for guaranteeing the overall charge neutrality of the system formed by the oppositely 
charged collapsing core and the remnant. It has been pointed out in section IXXVII that this condition of charge 
separation between the collapsing core and the remnant occurs only during the relevant part of the gravitational 
collapse process which, we recall, for a IOMq is of the order of 30 seconds. 

3. The necessity of developing a fine tuning in the final phases of thermonuclear evolution of the stars, both for 
the star collapsing to the black hole and the surrounding ones, in order to explain the possible occurrence of 
the "induced gravitational collapse" . 

New regimes are as well encountered from the point of view of nature of GRBs: 

1. The basic structure of GRBs is uniquely composed by a proper-GRB (P-GRB) and the afterglow. The most 
general GRB contains three different components: the P-GRB, the E-APE and the rest of the afterglow. The 
ratio between the P-GRB and the E-APE intensity and their temporal separation is a function of the B parameter 
(see Figs. ITTIITT)) . The best fit is obtained for i? = 3.0 x 10"^ (see section IXVj) . We recall that in the present 
case for B < 2.5 x 10"^ the energy of the P-GRB would be larger than the one of the E-APE and the energy of 
the dyadosphere would be mainly emitted in what have been called the "short bursts" , while for B > 2.5 x 10^^ 
the energy of the E-APE would predominate and the energy of the dyadosphere would be mainly carried by the 
ABM pulse and emitted in the afterglow. 

2. The long bursts are then simply explained as the peak of the afterglow (the E-APE) and their observed time 
variability is explained in terms of inhomogeneities in the interstellar medium (ISM) . The difficulties encountered 
by all theoretical models, through the years, in order to explain the so called "long bursts" are resolved in a 



drastic way (see section |XVI|) . The so called "long bursts" are not bursts at all. They represent just the E-APE 
which was interpreted as a burst only due to the noise threshold in the BATSE observations (see Fig. The 
E-APE is emitted at distances from the EMBH in the range 1.0 x 10^^ ~ 1.0 x IQ-^'' cm, see Tab.|D namely well 
outside the size of the progenitor star and already deep in interstellar space. The fact that the crossing of such 
distance, which is a typical dimension of an interstellar cloud, appears to occur in arrival time in only ~ 100 
seconds is perfectly explained by the rclativistic transformations encoded in the RSTT paradigm corresponding 
to a gamma factor between 100 and 300 (see section Ivl and Tab. IJ). This effect would be interpreted within 
a classical and incorrect astronomical picture by a "superluminal" behaviour propagating at ~ 3.6 x lO'^c (see 
Tab.©. 

3. The short bursts are identified with the P-GRBs and the crucial information on general rclativistic and vacuum 
polarization effects are encoded in their spectra and intensity time variability. In the limit S — + the entire 
dyadosphcrc energy is emitted in the P-GRB. These events represents the "short bursts" class, for which 
the afterglow intensity is smaller than the P-GRB emission and below the actual observational limits (see 
section IXIip . It is interesting that the proposed differentiation between the "short bursts" and "long bursts" 
within the EMBH theory is merely due to the amount of baryonic matter in the remnant, described by the 
B parameter, and totally independent from the process of gravitational collapse which is clearly identical in 
both cases. This explains at once the recently found conclusion thatthe distribution of short and long GRBs 
have essentially the same characteristic peak luminosity (Schmidt |l7Clj '). Also the result expressed in Fig. |H| 
that the average temperature corresponding to the P-GRB emission does increase for decreasing values of the 
B parameter can explain the observed fact that the "short bursts", which are obtained in the limit B ^ 0, are 
systematically harder than "long bursts" (Kouveliotou et al. 

A new class of space missions to acquire information on such extreme new regimes are urgently needed. The 
detailed observations of the yet unexplored region in the range up to 10 seconds in Fig. 1291 and the corresponding 
observations of the "short bursts" by a new class of space missions with higher sensitivity than the BATSE instrument 
appear to be of great importance. Such observations should allow to directly observe for the first time the general 
rclativistic and extreme quantum field theory effects connected to the process of formation of the EMBH. It can be 
of some interest to explore the possibility of observing in these regimes the "gravitationally induced electromagnetic 
radiation" (Johnston et al. and the "electromagnetically induced gravitational radiation" (Johnston et al. Jsl]) 
phenomena as well as to explore the possibility of developing neutrino detectors. This will need further developments 
of the predictions of the EMBH theory in these general rclativistic and ultra-high-energy particle phenomena. 
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